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SOUND PROPAGATION IN AND RADIATION FROM ACOUSTICALLY LINED FLOW DUCTS: 
A COMPARISON OF EXPERIMENT AND THEORY 


by Harry E. Plumblee, Jr., Peter D. Dean, 
George A. Wynne and Robert H. Burrin 
The Lockheed-Georg i a Company 


SUMMARY 


An experimental and theoretical comparison of many of the fundamental de- 
tails of sound propagation in acoustically lined flow ducts comprises the main 
body of this study. 

A description of the theory used for comparison with the experimental re- 
sults is included as the first section of this report. The theory of sound 
propagation in acoustically lined annular ducts is outlined and the method of 
solution is presented. Methods for calculating energy dissipation in lined 
flow ducts are discussed. Then, the theory for determining the complex radia- 
tion impedance of higher order modes of an annulus is outlined. 

A facility at the Lockheed-Georg i a Company Research Laboratory, designed 
for the explicit purpose of testing the predictions of the theory of sound 
propagation in flow ducts, is described. Methods of simulating acoustic duct 
modes and the fundamental characteristics of fan-tone noise are described. 

A detailed measurement program on sound propagation in rigid wall annular 
ducts with and without air flow is discussed. Techniques are described for 
measuring cut-on frequency, modal phase speed and radial and angular mode 
shapes. Measured effects of flow velocity on cut-on frequency and phase speed 
are presented. Comparisons are made with theoretical predictions for all 
effects studied. 

The theory of the two microphone method of impedance measurement is out- 
lined. Measurement error effects are discussed. An experimental program and 
results are discussed which detail the effects of flow on the impedance of 
locally reacting acoustic liner materials. Acoustic impedance for five sets 
of liners is presented under various flow conditions. 

A numerical study on sound propagation in annular ducts with one or both 
sides acoustically lined is presented. Following this, an experimental pro- 
gram which presents the measured values of attenuation, phase speed and mode 
shapes for all five liner sets is discussed in detail. Comparisons are made 
with theory where possible. 


A method for measuring the modal radiation (or termination) impedance 
from a finite-length annular duct is described and experimental results are 
presented and compared with theory for the case without flow. 

Finally, a series of radiation patterns from an annular duct, measured in 
an anechoic room, are presented. This experimental study shows the effect of 
mode order, frequency, jet flow and multiple mode pressure distribution on the 
radiated field patterns. 


I. INTRODUCTION 


Aircraft engine acoustic duct liners have been successfully designed for 
reducing the level of internally generated noise radiated to the community. 
Three large commercial passenger aircraft which incorporated community noise 
specifications into the design are now operational. These aircraft, the 
Lockheed L-1011, the Douglas DC-10 and the Boeing 7^7 have successfully passed 
the FAA certifications in force at the time of their initial deliveries. 

These aircraft are positive evidence that acoustic materials can, with reason- 
able assurance, be used to reduce radiation of internally generated engine 
noise. The question for the future is not, ‘Can internally generated noise be 
effectively reduced to reasonable levels? 1 , but instead, 'how efficiently can 
the acoustic materials be utilized in reducing internally generated noise? 1 . 

The primary goal of this research program was to provide data which would 
imply answers to the second question above. The program was not a research 
study on optimization of acoustic liner materials, as the above stated goal 
could imply. It was instead, primarily an experimental program, designed 
around a unique flow duct facility, to provide data points for correlation 
with theoretical predictions. The rationale was, 'provided the experimental 
data could be used to verify the appropriate theory, then the theory could be 
used for design optimization studies'. This, in simplistic terms, is the 
basic motivation underlying the study to be described in what follows. 

The subject of sound propagation in ducts holds great current interest. 

In the eight year period, 196A - 1972, twenty-six papers on the subject were 
published in the Journal of Sound and Vibration. To show that the interest is 
not waning, fifteen of these were published in the year 1971. Of all the 
papers surveyed, only a few contained experimental data and very few had data 
suitable for detailed correlation with theory. 

There were many aspects of the theory of duct acoustics involving assump- 
tions which needed experimental verification. One aspect which required vali- 
dation involved effects of parameters, such as geometry, flow velocity, tem- 
perature and liner impedance, on the wave properties in the duct. Another 
aspect involved validation of the mathematical representation of the physical 
problem as an eigenvalue problem. Still another aspect of the subject was the 
radiation of sound pressure from the termination of the duct through the jet 
flow to the surrounding ambient medium. In this case, there were no represen- 
tative theoretical results with which to correlate, thus the experiment should 
be used as a guide for future theoretical studies. 
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In particular, a facility was needed to provide experimental data for de- 
termining the effects on sound propagation within a duct of realistic geometry 
(circular or circular annular), liner impedance (both outer and inner walls in 
the case of the annulus), source sound pressure distribution, finite-length, 
mean flow velocity and velocity profile and temperature. Since most of the 
theoretical efforts were concentrating on various eigenvalue approaches, the 
data was required to be comprehensive enough to show effects of the above 
listed parameters on both the complex propagation constant and mode shape. 
Additionally, information on radiated noise from the duct was very sparse, 
especially for correlation with in-duct conditions. Thus, another feature of 
the tests was to provide radiated noise data which showed the effect of some 
of the above parameters and which could be used in future correlations with 
theory. 

The work reported on herein was conducted in the Lockheed-Georg i a Company 
Research Laboratory in a facility specifically designed to satisfy the re- 
quirements of the stated goal. The work included the development of tech- 
niques for operating the facility as well as for obtaining the desired data. 

The work presented involved detailed measurements of the sound transmis- 
sion through and radiation from hard wall and acoustically lined annular ducts 
with and without mean air flow velocity. It was necessary to develop techni- 
ques for measuring the impedance of acoustic materials in the presence of high 
velocity air flow. Also, techniques were developed for measuring the modal 
radiation impedance in the presence of air flow. Additionally, the radiated 
acoustic field for a series of duct modes was measured in an anechoic room. 


II. THEORETICAL SUMMARY 


The theoretical study of excitation, transmission and radiation of sound 
from sources in ducts has elicited considerable interest in the published 
literature during the past few years. At first, there were numerous engineer- 
ing studies which utilized the Morse theory for rectangular ducts (ref. 1). 
Other studies then presented numerical solutions for sound propagation in 
lined circular ducts, such as that of Rice (ref. 2). In an attempt to more 
adequately represent the physical situation usually encountered, Mungur and 
Gladwell (ref. 3) incorporated the ideas of Pri dmo re-Brown (ref. A) and Tack 
and Lambert (ref. 5) in a study of the effects of shear flow on sound trans- 
mission through the duct. In a later study, Mungur and Plumblee (ref. 6) 
studied sound transmission through a lined annular duct carrying a sheared 
mean flow. Several recent studies have presented numerical variations on the 
reference 6 approach, but have incorporated the same basic assumptions. 

All of the studies mentioned were concerned with sound propagation in an 
infinite length duct with uniform axial characteristics. In order to more 
closely represent the actual physical situation encountered in practice, that 
of a fan in a circular or annular duct, it is necessary to take into account 
source characteristics, the effect of finite duct length on the sound field, 
and finally, the radiated sound field from the duct termination. An approach 
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to the complete problem was given by Sofrin and Tyler (ref. 7) and Morfey (ref. 
8 ), based on concepts introduced by Rayleigh (ref. 9). These studies , primar- 
ily analytical but also including some experimental work, were made for hard 
wall ducts. But the concept of considering the source, the duct and the radia- 
tion field as separate acoustic elements was outlined and justified in detail. 

An entirely different approach for calculating the radiated sound field 
for a source distribution within the duct used a mathematical method based on 
the Weiner-Hopf technique for the solution of partial differential equations. 
Examples are outlined in a text by Noble (ref. 10) and problems concerned with 
sound radiation from unflanged circular pipes have been solved by Levine and 
Schwinger (ref. 11) and Carrier (ref. 12). These techniques have been ex- 
tended by Lansing, Zorumski and their co-workers at NASA (ref. 13 , 1^, 15 ) to 
include effects of mean flow through the duct and axial discontinuities in the 
liner impedance, as well as computing exact radiation patterns for representa- 
tive pressure distributions within the duct. In these ca 1 culat ions , the 
radiation field does not contain the very significant effects of refraction 
due to the jet flow, and it is questionable that the convection effects in the 
radiation field calculations are representative. 

The first method, that of analyzing each effect in relative isolation, 
appeared to be most attractive in the early theoretical studies at Lockheed, 
because of the building block approach which could be adopted in preparing 
computer analyses. For this reason, and since the Weiner-Hopf method becomes 
intractable when realistic duct and jet flows are considered, the first 
approach was adopted. 

The following sub-sections outline the theory which is to be evaluated 
for comparison with experimental results in the sections to follow. The theory 
is classified into sub-sections on sound propagation in lined ducts, energy 
flow in lined ducts and modal radiation impedance. This presentation gives 
only a brief outline of the derivations. Amplification of the details can 
be obtained from the relevant reference material. 


11.1, Sound Propagation in Ducts 


The following presentation begins with a summary of the derivation of the 
convected wave equation as presented by Mungur and Plumblee (ref. 6 ). The 
starting point for this analysis, as with most ae roacous t i cs problems, is with 
the Nav i e r-Stokes momentum equations, the mass continuity equation and the 
energy equation. For the purposes cf this study, a cylindrical coordinate 
system is appropriate. The variables are defined in the sketch at the top of 
the following page. 
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Cylindrical coordinate system used for wave equation 


The basic assumption used in deriving the linearized acoustic equations 

a re : 

(1) viscous effects are negligible, 

(2) products of fluctuating quantities are negligible, 

(3) the mean flow is axial and is uniform in that direction, 
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(5) radial gradients of the axial mean flow are permitted. (i.e. ^ 0). 
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If the density, pressure and velocity variables are perturbed and the equations 
resulting from time averaging the fully perturbed equations are subtracted, the 
following linearized equations describing acoustic motion for the specified 
physical situation result. The linearized momentum equations are, 
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where V is the mean axial velocity, v , v z and v Q are the acoustic velocity 
components, p* is the sound pressure and p is the mean density. The linear- 
ized mass continuity equation is, 
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where p 1 is the acoustic density component. 

From this point, the objective is to derive a wave-type equation by elim- 
ination of the acoustic velocity terms. If one differentiates equations (l), 
(2) and (3) with respect to z, r, and e respectively and sums the resulting 
equations, the single equation which follows is 
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If one observes that the term in curly brackets is the divergence in cylindri- 
cal coordinates less a v /r term and corresponds to the bracketed term in the 
continuity equation, a useful substitution can be made. The result of com- 
bining equation (4) with (5) is 
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The previous assumptions concerning neglect of viscous effects in conjunction 
with an assumption of negligible thermal conductivity allows the energy 
equation to reduce to 
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with c being the speed of sound in the fluid. For a particular frequency, w, 
variables are assumed to have the following axial dependence 
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The complex propagation constant, K, is represented as K = 3 - ia and k = w/c. 
From equation (2), using the assumption in equation (8), it is determined 

that 


v r (r>9) ” " ip o ck(l-MK) 

where M ■ ( r)/c. Substitution of equation (7), (8), and (9) into equation 
(6) results in the following 1 convect i ve- ref ract i ve 1 wave equation, as given 
in reference 6, 
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The effect of convection is represented by the terms, k 2 (-2MK + M 2 K 2 )p' and 
the shear refraction is accounted for by the term 


2K 0M 3p 1 
1-MK dr 9r * 


If the boundary conditions on the inner and outer radii are independent of 
angle, 0, then the equation is separable. Also, it is convenient to non- 
dimens i ona 1 ize the radial variable, r, with respect to a convenient length. 
For the case of a cylinder, the radius is appropriate and in the case of an 
annular duct, the duct width is a desirable parameter. The new variable is 


X = r/L. 


The separation variables are p*(X,6) = P(x) 6(0). The resulting angular 
differential equation is 



» 


( 12 ) 


7 



from which it is determined that 


0 = 


sin m0 
cos mQ 


(13) 


and k. = 0 , 1 ,2 , . . . ,m , . . . ,N . 

U 

The resulting radial differential equation is 
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The solution of this equation yields the eigenvalue, K = $ - ia, for a parti- 
cular set of boundary conditions, angular mode number, m, frequency, kL, and 
Mach number profile, M(X). 


Several techniques are available for determining the solution. For the 
idealistic case of plug flow, an exact analytical eigen solution exists. This 
solution has been discussed in detail for a cylinder, with and without flow, 
by Rice (refs, 2, 16) (but with a flow boundary condition which assumes con- 
tinuity of particle velocity rather than particle displacement (see Mungur 
(ref. 17) and Tester (ref. 22)). 

For the case with shear flow, Mungur and Plumblee numerically integrated 
equation (14) and utilized an eigenvalue search procedure in the complex 
eigenvalue plane. The initial numerical program was relatively inefficient 
and in many cases, skipped modes. Mungur later improved this numerical 
routine by incorporating a steepest descent prediction for locating eigen- 
values. 


In a more recent study of equation (14), Wynne and Plumblee (ref. 18) 
arranged the finite-difference equations in the form of an eigenvalue problem. 
The difficulty in this approach resulted from the fact that the eigenvalue 
formulation was non-linear. The form of the eigenvalue equation, to which the 
finite difference was applied, is displayed in equation (15), which is a 
slight rearrangement of equation (14); 
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The eigenvalue, x , is defined, from equation (14), as 


X n = - (kL) 2 (1 - K 2 ). 
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The finite-difference equation which results from a simple partitioning scheme 
i s 
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The impedance boundary conditions are applied at station and X . In the 
finite-difference scheme utilized, all variables were defined at 5:he central 
point of the finite-difference increment. The program was written for varia- 
ble increments so that the step size could be optimized in regions of rapid 
change in the variables. The derivative is 
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The form of the eigenvalue matrix, from equation (17) is 


{[A] - A n [|]}P = 0 


( 21 ) 


with the A matrix being tri-diagonal as observed in equation (17). 

In order to determine the complex eigenvalues in the presence of a 
sheared flow, M.(X.), an iterative procedure was established. First the no- 
flow e t gen va 1 ue" 1 is^ introduced into the [A] matrix as an initial guess and the 
eigenvalue problem is solved for the eigenvalue of smallest magnitude. This 
first solution incorporates some flow effects. This new eigenvalue or a 
suitable modification thereto by the Newton-Ra 1 phson scheme, is then substi- 
tuted into [A] matrix once again, after convergence tests have been performed. 
The iteration is continued until a converged value is obtained for the lowest 
eigenvalue. This procedure is repeated until the desired number of solutions 
is obtained. 

Needless to say, convergence problems were encountered at many stages in 
the solution process. One of the more rudimentary techniques used for coun- 
teracting the slow convergence was a gradual introduction of the flow magni- 
tude and shear gradients. The region where most convergence problems occurred 
was in the vicinity of a branch point. Most of these were cured, but there 
are still certain conditions where the numerical program allows the solution 
process to oscillate between two eigenvalues, across the branch line. When 
this occurs, the numerical solution is very slow and sometimes fails. 

The solutions which result from the procedure described above also con- 
tain the mode shapes appropriate to the set of eigenvalues. It is noted that 
the concept of modes is very useful from the standpoint of a mathematical 
examination of sound transmission through ducts, however, the problem of phys- 
ical concern is that of propagation of a source pressure or acoustic velocity 
distribution. Thus, while individual modes are, in many instances, physically 
meaningless they are required as a means of representing the behavior of a 
complex sound pressure distribution as it propagates through the duct. 

Based on the above, the numerical analysis was written such that propaga- 
tion of an arbitrary sound pressure distribution could be studied. In order 
to simulate the complex radial pressure distribution, a least square fit to 
that pressure distribution was accomplished, making use of the set of radial 
modes for a particular angular mode. The error function for the least square 
representation is 
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t h 

where P (r.) is the measured pressure distribution for the nr angular mode at 
the i t ^ ,m raclial station and <J> ( r t ) is the complex mode shape. The mode co- 
efficients are determined by m mi n imiz ing the total error magnitude with respect 
to the coefficients. Thus, if 
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where * denotes the complex conjugate and where r j has been included as the 
appropriate weighting function for the cylindrical coordinate system, then 
the conditions 
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gives a set of simultaneous equations from which A™ can be determined. A 
similar set of conditions could be used to match an acoustic velocity source 
distribution with the complex mode shapes. 

This forms the basic analysis for sound propagation in circular or annu- 
lar lined flow ducts. A number of similar analyses utilizing different 
numerical approaches have been done within the past three years. These have 
been different in the solution technique but have generally utilized the basic 
assumptions leading to equation (14). Thus, further discussion of the similar 
reference material is not felt to be necessary. 


11.2. Energy Dissipation in Acoustically Lined Flow Ducts 

The results from the previous analysis yield eigenvalues which define the 
attenuation rate of individual modes. However, as will be demonstrated in a 
subsequent numerical study, the modal attenuation rate can be very misleading 
when trying to predict the performance of a liner. The overall power attenua- 
tion can differ from the attenuation of the least damped mode for at least two 
reasons. The first results from reflections from the termination and the 
second is a result of interference between multiple propagating modes. Other 
causes of variations between energy attenuation and attenuation of the least 
damped mode are due to source sound pressure distribution, source impedance 
and refraction due to inlet or exhaust flows. Some of these have been intro- 
duced in the numerical procedures but termination and source impedance and 
inlet and exhaust flow effects were not included. 

The subject of acoustic energy in non-uniform flows has been discussed 
extensively by Morfey (refs. 19 and 20), Mohring (ref. 21) and Tester (ref. 
22). Cantrell and Hart (ref. 23) derived an expression for acoustic energy 
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which, Morfey (ref. 19) showed, is conserved in non-uniform but irrotational 
mean flows. Morfey also included in his studies a description of energy pro- 
duction terms due to interactions in rotational (sheared) mean flows. Tester 
used the work of Morfey and that of Mohring to examine certain "strange" modes 
(modes which have phase velocities in the direction opposite to the direction 
of decay), and showed that the Cantrell and Hart energy flow appears to be 
amplifying in the uniform flow region, whereas the Mohring energy flow in the 
thin shear layer is in a d i recti on oppos i te to that in the uniform flow region 
and compensates for the apparent energy "amplification." Tester also showed 
that the energy of well cut-on modes could be determined with the Cantrell and 
Hart energy expression. 

Although the work of Mohring results in a definition of energy which has 
no production terms in a parallel sheared flow, it is not valid for a match 
condition between a duct and the radiated field. The work would require 
reformulation to represent the energy in free-field with either an irrota- 
tional or a rotational flow. However, as long as the flow is irrotational, 
for instance an ideal inlet flow, Cantrell and Hart‘s energy will be conserved. 

In the presentation which follows, the Cantrell and Hart energy expres- 
sions as given by Morfey are used for determining power level attenuation. 

The results are expected to be valid for modes which are cut-on, but could be 
questionable for non-cut-on modes. 

The form of Morfey 1 s energy is 
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For the axial flow of the present analysis, the axial intensity is 
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where * denotes the complex conjugate. It is necessary to derive an expression 
for the axial acoustic velocity, v . This is accomplished by substituting 
equation (9) and (8) into equation (1), with the result for the n^.mode being, 
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After substitution of equation (27) into (26), the total power at a particular 
cross section can be found with the following integral relationship, 
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where X- is the inner radius and X^ is the outer radius. The ratio of the 
powers at two different cross-sections is taken to be the parameter which most 
closely represents the attenuation performance of a particular liner design in 
the radiated field from a lined duct if effects such as flow refraction, 
diffraction and reflections at the end are neglected. 


11,3, Modal Reflection Coefficients for Finite-Length Ducts 


A very elegant analysis for determining the modal reflection coefficient 
in an arbitrarily shaped duct of finite length was performed by Doak (ref. 

24). The modal reflection coefficients, in this instance, can be used to 
completely specify the sound field within the duct and in the radiated field, 
based upon a specific source sound pressure or acoustic velocity distribution. 
Thus, for practical purposes the specific modal impedance coefficient, re- 
sulting from the modal reflection coefficient is equivalent to the radiation 
impedance for a piston of the same cross-section as that of the duct termina- 
tion in a rigid baffle and with a velocity equal to the modal velocity 
distribution over the duct cross-section. 

The analysis by Doak, which is a more general study than that to be sum- 
marized in the following subsection, leads to a derivation of the modal 
reflection coefficients for an arbitrarily shaped f i n i te- 1 ength duct with 
radiation from the termination at either end. The analysis is only applica- 
ble to orthogonal duct modes, but as is pointed out in ref. 24 any set of 
non-o rt hogonal modes can be o rthogona 1 i zed . Thus, in principle at least, the 
analysis can be applied to ducts of arbitrary geometry where the modes would 
not normally be orthogonal. 

The method, as described in the following, was applied to a rigid-wall 
duct of annular cross section by Haddle (ref. 15). Calculations were per- 
formed by numerical evaluation of a t r ip le- integra 1 over the surface of the 
annulus. The results for the real part of the impedance were compared with 
calculations by Morfey (ref. 26) and shown to be identical, although the 
methods of analysis in the two cases were considerably different. 

11.3.1. Method of analysis 


It is shown in ref. 24 that the sound pressure, P“(X.;aj), inside a 
duct can be represented as (in the absence of mean flow, V = J 0). 
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where iJj satisfies the equation, 
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and the boundary condition, 
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on the duct boundaries (the (-) superscript refers to conditions within the 

duct). These boundary conditions permit determination of the wave numbers, 

k: , from which the axial wave number k can be determined as 
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The axial acoustic particle velocity from equation (1) for the case of no 
mean flow is, after spectral transformation. 
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Substitution of equation (29) yields the following relationship, 
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where the complex phase velocity, c , is related to the wave number, k , by 
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At this point, the modal reflection coefficient is defined as the amplitude of 
the reflected wave to the amplitude of the incident wave, (the possibility of 
reflection in several modes with only one mode incident occurs and in some 
cases is very important; see Lansing, et al., ref. 14 ) 


R . = b JH . ( 36 ) 
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The analysis which follows ignores the modal reflection coupling coef- 
ficients and considers only the self-modal reflection coefficients. Thus, for 
the mn^-h mode, it can be shown that at the end z = 0, the inner solution can 
be used to define a modal impedance of the following form: 
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Doak then proceeded to determine an acoustic solution appropriate to the 
region external to the duct, i.e., a solution which satisfies the boundary 
conditions in the exterior region, but utilizing the set of orthonormal char- 
acteristic functions of the duct termination cross-section. The pressure at 
an arbitrary point in the acoustic field of a piston in a rigid baffle is 
formulated in terms of the Green's function for a unit point source on the 
surface of a plane, rigid baffle of infinite extent. 


The general Green's function solution, 


( X j * y. 


i s 


P6 ( yy = l « mn (X 3 ;x|) ^ mn (Xi,x 2 ) ^ mn (x-,x-) (38) 

J J mn 

where the (+) superscript refers to locations in the outer region, and the 
coef f i cients a are , 
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The Green's function appropriate to the case of interest is 
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P*(Xi,X 2 ,Jt,Xi,X 2 ,A) ~ *^P ( 
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Thus, the acoustic pressure in terms of the Green's function and the acoustic 
velocity on the outer surface of the duct can be shown to be (Haddle, ref. 25) 


iklX-X', , 

+ . . v r re 1 1 + 

p = -!P o ck l [ I 2ti x _x> v m 
mn v S-' 


dX ' dX ' 
mn 12 


(41) 


with the nomenclature defined in the sketch below. This step is equivalent to 
that discussed in Morse and Ingard (ref. 27) for the circular piston. 
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The solution for the reflection coefficient, R , can be accomplished with 
application of the condition of continuity of sound pressure and particle 
velocity across the interface between the duct and the radiated field, i.e.; 
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Use of the orthogonality relationships for the duct modes is necessary to 
complete the evaluation. Also, the modes are normalized such that 
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where T = E (k a). The E functions were defined by Sofrin and Tyler (ref. 
^ mn mn v mn 

Upon performing the required manipulations, the modal termination impe- 
dance can be formulated as 
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In this expression, N mn is the 
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The nomenclature for the annulus is defined in the sketch below. 
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It was possible to reduce the integral expression of equation (kk) to a triple 
integral by changing variables such that 
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and performing the indicated integration with respect to 0‘ . The resulting 
integral, which was evaluated numerically is 
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where r = r (a, a 1 ,<j>) . 

Haddle's numerical results were compared with Morfey's (ref. 26) calcula- 
tions for the radiation resistance for hub/tip ratios of 0.0, 0.3 and 0.5 and 
were found to be in agreement. Some of these results will be shown in a later 
section for the experimental configuration discussed in this report. 


Mi. EXPERIMENTAL GOALS AND FACILITY DESCRIPTION AND OPERATION 


As stated briefly in the introduction, and further indicated in the 
theoretical synopsis, the primary goal of this program was to .provide experi- 
mental data in sufficient quantity and detail and with sufficient accuracy to 
verify the major predictions of the theory for sound propagation in a lined 
aVinular duct with mean' flow, - In the implementation of 1 this go a 1 , a f aci 1 i ty 
bliT 1 tat' the Toe" khee d - G e o rg i a Company, with Lockheed fixed asset and IR and D 
funding, was used exclusively. The following subsections include a detailed 
explanation of the program goals and an explicit description of the facility 
and the techniques developed for operating the facility. 


I 11.1. Major Objectives 


There were several major areas for detailed research, as well as subsid- 
iary research tasks necessary for calibration of the facility and complete 
definition of the test conditions, such as determination of wall impedance and 
radiation impedance. A basic aspect of the primary goal was to be able to 
measure the complete sound field, at a discrete frequency, within a finite- 
length acoustically lined annular duct sustaining a mean flow velocity, for 
comparison with theoretical predictions. However, the approach to implemen- 
tation of this objective was considerably more basic than simply launching 
into a detailed measurement program. 

It was shown in the previous section that solutions to the convected wave 
equation, which describes the duct sound field, are in the form of eigenvalues 
(the complex modal propagation constant which defines attenuation rate and 
phase speed) and eigenvectors (the complex mode shape or sound pressure dis- 
tribution in the duct corresponding to a particular eigenvalue). It is a 
basic premise that these fundamental solutions or functions can be super- 
imposed to represent any arbitrary pressure distribution at any point in the 
duct. Validation of these basic solutions provides proof sufficient for pre- 
diction of the complete sound field. 

The experimental program was sub-divided into four major areas for 
research, with subsidiary research investigations into mode simulation, wall 
impedance measurements and radiation impedance measurements. The four basic 
investigations were measurement of sound pressure wave propagation 
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characteristics in (1) a rigid wall annular duct containing an ambient medium, 
(2) a rigid wall annular duct with an axial mean flow velocity, (3) an acous- 
tically lined annular duct containing an ambient medium, and (4) an acousti- 
cally lined annular duct with an axial mean flow, A brief description of each 
of the listed investigations, as planned, follows, 

111.1,1. Rigid wall duct without flow 

This task, in addition to providing basic information for comparison 
with theory, was to serve as a calibration for the facility and as a test bed 
for developing operational and measurement techniques. First, methods for 
simulating (or exciting) angular standing modes (acoustically identical to 
spinning modes) were to be developed. Then a number of angular modes were to 
be excited individually and the mode shapes, phase speed, axial standing wave 
ratio and radiated field were to be determined. Excitation frequencies were 
to encompass the first three radial modes. 


I I 1 . 1 . 2 . Rigid wall duct with mean flow velocity 

The survey of mode characteristics initiated above was to continue 
in the presence of axial mean flow velocity. As designed, the facility could 
provide a sustained velocity of 200 fps with blow-down capability of 450 fps 
for short runs (on the order of one minute). The useful test range with flow 
was to be determined. 


I 1 I . 1 . 3. Finite wall impedance without mean flow 

A minimum of two sets of acoustic duct liners were to be constructed 
for measurement of attenuation, phase speed and mode shape for a frequency range 
encompassing the first three radial modes. The tests were to include detailed 
studies of the characteristics of individual modes and of more complex sound 
pressure patterns at a discrete frequency. 

111.1,4. Finite wall impedance with mean flow velocity 

This series of tests was to be similar to the one above, but would 
include mean flow velocity up to the maximum possible as determined in the 
earlier series of hard wall tests. 


III. 2. Facility Description 


The development of experimental techniques for evaluating acoustic liner 
design has followed a somewhat different course from that of the theoretical 
developments described in Section II. Experimental facilities have been 
primarily designed to give the average acoustic energy attenuation of 
Helmholtz type acoustic liner materials, developed over the past decade for 
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use in fan jet engine ducting. These facilities, termed flow duct facilities, 
have been built at the major airframe companies and extensive material evalu- 
ation tests have been conducted. Descriptions of some of these tests are 
given in NASA contractor summary reports (ref. 28, 29) and the proceedings of 
NASA Langley contractors conferences on Noise Alleviation of Large Subsonic 
Jet Aircraft (ref. 30, 31). Design considerations for the flow duct type of 
facility were discussed in detail by Melling and Doak (ref. 32). Melling also 
discusses the limitations and interpretations of test results from the double- 
reverberant chamber flow duct facilities. 

The major points of concern in using data from this type of facility for 
detailed confirmation of the analytical predictions are: 

o Sound pressure distribution at the duct entrance is uncontrolled 
but is generally diffuse. 

o The sound pressure in the source chamber is broadband noise in 
most cases (narrow band or discrete frequency peaks may exist). 

o The receiving chamber is reverberant. 

As a result of these facility characteristics, it is very unlikely that 
many of the details of the theory can be verified using flow duct facility 
data . 


As a result of the above considerations, a new type of facility was de- 
signed for verification of many of the theoretical trends and predictions. In 
particular, the facility was designed to overcome the above objections, thus 
permitting experimental determination of (1) the individual mode cut-on 
frequencies and shapes, (2) modal propagation constant, (3) duct termination 
impedance, and (A) the radiated sound field from the duct termination for 
single modes and complex pressure distributions. In addition, the facility 
was designed to permit evaluation of effects of mean flow and finite wall 
impedance on the above fundamental acoustic properties of the duct. 

The facility to be described was designed and constructed at the Lockheed 
Georgia Company's Research Laboratory in late 19^9 and early 1970. Initial 
testing techniques were developed under Lockheed IR and D funding during 1970. 

A basic plan view of the complete facility is shown in figure 1 and a 
side elevation of the flow duct is shown in figure 2. Figure 3 contains a 
sketch of the complete facility and a photograph of the source section. 
Basically, the facility consists of (with the numbered items corresponding to 
the circled numbers in figures 1 and 2): (l) a high-pressure air supply, 

(2) a control valve, (3) a wide-angle expansion, (A) a muffler section, 

(5) another wide-angle expansion, (6) an acoustically lined plenum, (7) a flow 
straightening and turbu lence- reduci ng honeycomb section, (8) a contraction, 

(9) a final muffler, (10) the acoustic source section, (11) the test section, 
(12) an anechoic room, and (13) an exhaust and muffler. A brief description 
of each of these components is given below. 
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(1) The air supply system (which is also used by the Lockheed wind 
tunnels) continuously delivers 20 lb/sec. of cleaned, dry air com- 
pressed to 300 psi. If required, a low-pressure pump can be added 
to the system for delivery of 30 lb/sec. at 100 psi. The air 
system also has blow-down capability of 12,000 cubic feet of air 
at 300 lb/in 2 . 

(2) The valve is a 6-inch throttling ball valve with a pressure servo- 
controller. It is designed to be most effective on pressure control 
for speed ranges of 50 to 450 ft/sec in the test section. 

(3) The wide-angle conical expansion connects the 6-inch valve to a ' 
20-inch muffler section. 

(4) The muffler section is 13 feet long and is lined with 2 inches of 
polyurethane foam, retained with a perforated metal sheet. The 
muffler also contains a full-length, similarly treated, cruciform 
sect i on . 

(5) The wide-angle expansion section has strategically placed screens to 
effect a smooth flow transition from the 20-inch supply pipe to the 
36 -inch plenum. 

(6) The plenum is made from a 36 -inch pipe and is lined with polyure- 
thane foam retained with perforated metal. The plenum contains a 
10-3/4 inch pipe which forms an annular flow region. 

(7) The delivery end of the plenum contains an 8-inch long section of a 
1/8- inch cell honeycomb for reducing the turbulence level in the 
supp 1 y air. 

(8) The contraction section connects the 36-inch plenum to a 19"inch 
inner diameter muffler section. 

(9) The muffler section is 4 feet long with a lining similar to that in 
the other acoustically lined sections. 

(10) The noise source section also shown in figure 3 is equipped with 16 
loudspeaker driver units spaced at 22.5° intervals and connected to 
the outer wall of the source section. The ends of the driver tubes, 
at the juncture with the source section, are covered with 20-Ray 1 
felted metal to minimize acoustic flow resonance and to prevent 
turbulence generation. The driver units are individually controlled, 
both with regard to amplitude and phase, and the useful range ex- 
tends from 300 to 6000 Hz. Electrical power input to the driver 
units is up to 100 watts peak. The source section has a joint, 
pressure-sealed with Teflon, immediately downstream of the driver 
tubes. The downstream section is rotatable with respect to the 
driver units, permitting acoustic pressure or velocity to be mapped 
azimuthally, using transducers fixed in the downstream source 
section or in the test section, which is fixed to the rotatable part 
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of the source section. Radial traversing of transducers is also 
achieved. Two configurations of the microphone traversing units 
are shown in figures k and 5. 

Figure k is a photograph of the basic radial traversing unit which 
was designed for automatic traversing and recording of microphones, 
hot-wire anemometer probes, pitot tubes and thermocouples. It con- 
tains a variable speed D.C. motor which drives the probe through a 
belt and a linear wire-wound potentiometer for generation of a D.C, 
voltage proportional to frequency. Figure 5 is a modi fi cat ion of 
the basic traversing mechanism for obtaining SPL measurements in 
the region upstream of the lined test section. It ut i 1 i zed two of 
the traversing mechanisms, one mounted at the duct termination and 
another mounted in the test section. A one-half inch microphone 
with a long probe tube attached was fixed to the ends of the radial 
traversing mechanisms. Axial locations were obtained by way of 
manual positioning of this device. An automatic axial traversing 
microphone would have been desirable for some of the tests, however, 
construction and installation was not practical considering cost and 
schedule commitments. 

(11) The test section consists of a basic steel shell three feet in 
length with a 23 inch inner diameter. The center pipe has a 10-3A 
inch outer diameter. Liner materials are placed on the inner and 
outer walls of this section as desired. In all the tests reported 
herein, the i.d. of the outer liner remained at 19 inches (this 
provides for application of material up to 2 inches deep on the 
walls). The o.d. of the inner liner was 10-3/^ inches for all 
tests. Provision was made for inserting transducers with radial 
drive capability, as discussed, on a 1-inch axial spacing. Azimuthal 
mapping is accomplished by rotating the test section with respect to 
the driver units. 

(12) The anechoic room illustrated in figure 3, was converted from a re- 
verberant room by adding 6 inches of polyurethane foam to all 
interior surfaces. The room is anechoic for frequencies above 500 
Hz. 


(13) An exhaust port is provided in the wall of the anechoic room as 

shown in the photograph of figures 6 and 7. A 1/10 scale flow model 
of the room was constructed and tested to determine an optimum ex- 
haust port configuration for minimizing the flow circulation in the 
anechoic room. It was finally determined that the collector shown, 
with a 20° angle from a normal to the wall, resulted in minimum 
circulation within the major acoustic measurement volume in the 
room. 

Before the tests were completed, five sets of acoustic liners were tested. 
Two of these liner sets are shown In figures 8-11. The liner set illustrated 
in figures 8 and 9 was constructed of a felted metal acous t i c . mate ri a 1 over a 
1/2 inch plastic honeycomb with 3/8 inch cells. Two sets of inner and outer 
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Figure 9 . 20 Rayl Outer Liner Installed in Test Duct 
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Figure 10. Welded Fibrous Metal Outer Liner 



Figure 11 


Welded Fibrous Metal Outer Liner Illustrating 
Holes for Probe Microphone 


liners of this design were constructed, utilizing a nominal 20 Rayl acoustic 
material on one set and a nominal AO Rayl acoustic material on the other. 

The liner illustrated in figures 10 and 11 consisted of a single unit for the 
outer wall of the test section and used a welded fibrous type acoustic materi- 
al over a 1-1/2 inch honeycomb with a 1/8 inch diameter aluminum cellular core. 


I I 1.3. Operational Details 


III. 3.1. Instrumentation and Controls 

The controls are divided into three sections, viz, air supply, microphone 
position, and acoustic excitation and are illustrated with a block diagram in 
Figure 12. 

The air supply is automatically controlled by a 6-inch throttling ball 
valve and servo controller, with the velocity head monitored by a pitot probe 
upstream from the plenum, located in an area of low turbulence. 

The duct probe microphone position is controlled angularly with respect 
to the drivers by a motor which rotates the duct test section, the angle being 
monitored by a cable deflection box. This signal is input to the X-Y plotter 
and computer interface via suitable signal conditioners. The radial probe is 
motor driven and its position is monitored by an interconnected potentiometer, 
which is also connected to the plotter and computer interface. The boom 
microphone's angular travel is also detected by a deflection box. All micro- 
phone position controllers have limit switches to eliminate the possibility of 
excess travel and consequent equipment damage. 

The frequency control circuit for acoustic excitation utilizes a voltage 
controlled oscillator coupled to an analog computer which provides for auto- 
matic frequency compensation for temperature changes sensed by a thermocouple 
in the duct. The frequency controller has been designed to hold wavenumber, 
k = 2 tt f/c, constant, where c is the speed of sound in the duct and is, 
therefore, proportional to square root of duct temperature. 

The data acquisition system has two functions; the first is to provide a 
visual instantaneous record on an X-Y plotter; the second is to store data in 
digitized form on the disc storage of a UNIVAC Al8 computer. This data can 
then be manipulated as required, after being transferred to magnetic digital 
tape. The computer storage is controlled by a teletype and MAC 16 computer 
in the laboratory. A block diagram of the complete recording systems can be 
seen in figure 13* 

The probe microphone in the duct (set normal to the flow) is a standard 
2mm diameter B&K probe (suitably damped and modified to have sensing ports in 
the side instead of the end, as shown in figure’ 6A) , coupled to a standard 
1/2" B&K microphone. The reference microphone is a 1/A inch diameter B&K 
microphone mounted in the rotating section at an axial position such that it 
is projecting upstream and in the plane of the drivers. This microphone is 
used to set up the individual drivers for amplitude and phase such that their 
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combination will produce the desired circumferential standing mode (as will be 
described later). The boom microphone is a 1/2 inch B&K with a 2mm diameter 
"static pressure" probe attached. This type of probe minimizes the pressure 
fluctuations due to probe ' aerodynami c interactions. All the microphone out- 
puts are filtered through a 2 Hz tracking filter, which eliminates the broad- 
band components of flow turbulence. The resultant signal to noise ratio of 
the excitation peaks is better than 60 dB for zero flow tests, reducing the 
order of 10 to 20 dB at 200 fps mean flow velocity. However, unless one 
reverted to signal averaging or signal recovery techniques, the random fluc- 
tuations introduced by the flow noise made exact comparisons with theory 
difficult. Thus, all acoustic comparison tests were restricted to 150 fps 
f low vel oci ty . 


Ml. 3 . 2. Facility dimensions and characteristics 

The test section geometry and dimensions are illustrated in figure 
lA which also includes pertinent dimensions for the various acoustic liners 
which were tested. The test section is three feet long, the center pipe or 
liner is 10-3/A inches in diameter and the outer rigid pipe or liner has a 19 
inch inner diameter. This results in an annulus of 4-1/8 inches width and a 
hub-tip ratio of 0.566. 

Air flow velocity can be maintained at 200 fps but all tests reported in 
this study were limited to 150 fps. Circumferential velocity profiles 
measured in the test section are shown in figure 15. This shows approximately 
10% variations in mean flow velocity due to non-uniformities in the surface of 
the final muffler stage (figure 2, muffler (9)). 

The radial velocity profiles shown in figure 16 were reconstructed from 
the velocity measurements as described in the following discussion. 

Measurements of velocity cannot be made very close to the duct walls, 
primarily because probe errors become large in the high velocity gradient 
region. As a result, useful measurements usually cease some distance from the 
surface and a layer of large shear remains which cannot be easily explored 
exper i mental 1 y. 

In order to provide reliable estimates of the velocity gradients in the 
wall layer, calculations were made using the universal law of the wall (ref. 
38) for two-dimensional flow. 

The wall layer is idealized into a viscous sublayer, in which viscosity 
is assumed to be the only source of shear stresses, and a ful ly-turbulent 
layer in which viscosity is neglected; no blending region is admitted between 
them. 

In the viscous sublayer, the velocity, u, is proportional to the 
distance, y, from the surface: 


(47) 
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and where x is the wall shear stress, p the density of the fluid, and v its 
kinematic viscosity. 

In the fully turbulent region, the velocity varies logarithmically with 

y : 


U + = jj- (1ny + + C) , 


(50) 


where k is von Karman's constant, equal to about 0 . 4 1 8 , and C is a second 
empirical constant, equal to about 2.31 for smooth walls (ref. 34). 

The "edge" of the viscous sublayer is given by the intersection of equa- 
tions 47 and 50, and corresponds to a value of y + , over an appreciable frac- 
tion (.1 or .2, perhaps) of the duct width. Its outer boundary of validity 
is ill-defined, but it certainly extends out into the region where reliable 
probe measurements can be made. 

Both equations 47 and 50 contain x as a parameter, and its value is 
unknown at the outset. The procedure adopted for finding it was based on the 
assumption that several reliable data points lay in the region of validity of 
equation 50. It was then only necessary to choose a value of x w which ensured 
that equation 50 merged smoothly into the measured velocity profile. This 
method is equivalent to the use of a Clauses plot for determining wall shear 
stress (ref. 35). The values of x w obtained were of the right order of magni- 
tude for a duct of the appropriate Reynolds number. 

Once x w is known, the complete vel oci ty. p rof Me, in the region of 
interest, can be calculated. The velocity gradient can then be obtained by 
differentiation. In the viscous sublayer: 


3u _ T w 
3y p"u * 


and in the logarithmic region: 
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The velocity gradients calculated by this procedure are shown in figure 17. 


III. 3 . 3 . Sound pressure spatial control: preliminary study 

A primary requirement for the facility was the ability to set up a 
particular spatial pattern of sound pressure in the operating frequency range 
to simulate the sound pressure distribution in an angular mode. In the design 
study of the facility, it was decided that radial control of sound pressure in 
the source section would be impractical to achieve, thus the design focused on 
angular control in the source section. 

A detailed study on source characteristics was conducted. Two source 
distribution schemes were considered in detail. One was a crude method for 
exciting one higher order radial mode as illustrated In figure 18a. This 
scheme would result in maximum sound pressure level in the test section, but 
would pose difficult problems in spatial control. However, construction dif- 
ficulty and lack of confidence in the final control over the sound pressure 
eliminated this technique. It was finally determined that for the purpose of 
this facility the scheme shown in figure 1 8 b would be most practical, both 
from a construction and a control standpoint. This scheme assumed 16 loud- 
speaker driver units with independent level and phase controls for each unit. 

A simplified example, as derived by Doak (ref. 36 ) during a preliminary 
study of the facility design, can be used to illustrate the expected response 
in a hard wall duct from a phased array of transducers mounted in the duct 
walls, using rectangular geometry as being representative of a thin annulus. 
The purpose of the example is to show the expected enhancement of one mode 
relative to all other cut-on modes at a particular frequency w. 

For simplicity it is assumed that all transducers have the same dimen- 
sions (square, of edge length d) , all have their centers on the intersection 
of the x = o and z = o planes, and all are driven at the single frequency, 00. 

Let be the complex velocity amplitude of the transducer (see 
figure 19), with its center at (o,y£,0). Then the pressure field in the 
region z>d/2 is given by the following equation, as derived by Doak. 
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Suppose, for example, that it was desired to enhance the response in the 
(0,4) mode. One might do this by using five (but in effect, only four) trans- 
ducers, placed at (0,0,0), (0,b/4,0), (0,b/2,0), (Q,3b/4,0), and (0,b,0), as 
shown in figure 19, with respective amp 1 i tudes (v q /2 , -V Q , V Q , -V q &\/ o / 2). 
Strictly speaking, of course, a finite-sized transducer cannot have its center 
at (0,0,0), etc., in a rectangular duct. However, if the rectangular duct is 
considered as an analog of an annular duct, this difficulty does not arise. 

In other words, the idealized example being considered here is really that of 
a duct which is infinite in the y-direction, excited by an infinite number of 
transducers with centers at (0, ±£b/4,0). 

For this array of transducers £ V cos (mvy /b) is zero for the modes 
(0,0), (0,1), (0,2) and (0,3), but £ for the mode (0,4), it is 4V Q , The sum 
is similarly zero for all modes (m,n) for which n is less than 4. Hence only 
those modes (m,n) for which n is equal to or greater than four can be excited 
by this transducer array, whatever the frequency. 

In general, for this array, 
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For all odd values of n the sum is zero. For even values of n (n = 2m), 
the sum can be written as 
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If m is odd the sum is again zero. For m even, m=2k, the sum becomes 



2 Vo | 1 - cos (krr) j 


( 55 ) 


This is zero for k=0 , 2, 4, 6..., and equal to 4 V q when k=1, 3, 5.... 

Therefore, for this array of transducers, and independently of frequency, 
the amplitude of excitation of all (m,n) modes is zero except for the modes 
(m, 4k) where k=1 , 3, 5» .... For these modes the array excitation factor is 


/ n7T V o \ 

l V £ COS (~b — ) = ^ Vo * (56) 

The implication of this result is clearly, that as far as the (0,n) modes 
are concerned, a phased arrangement of transducers can provide excitation of 
only, say, the (0,j) mode, over a frequency range from the cut-on frequency 
for that mode,o> 0 »j» up to a frequency of an octave or more higher than the 
cut-on frequency, the range depending on the mode number concerned and also 
perhaps on the array geometry. In a particular case, as considered above, 
the frequency range can easily be determined. 

In practice, however, another factor is more likely to limit the range 
over which an isolated mode can be excited. If all (m,n) modes are considered, 
as well as just the (0,n) modes, then among the modes that can be excited by 
the above array, the (1,4) mode may well have a cut-on frequency lower than 
that of the (0,12) mode. In a rectangular duct of cross-section shape repre- 
sentative of an annular duct, one might have, say, a breadth equal to 1/8 the 
height (i .e. , a = b/8). Then the ratios of the cut-on frequencies concerned, 
namely (0,4) (1,4) and (0,12) are 1:2:3. Thus, the frequency range over which 
the (0,4) mode could be excited in isolation in this case is limited to the 
octave beginning with its cut-on frequency. However, since for this cross- 
section shape the cut-on frequency ratio (relative to that of the (0,4) mode) 
for the (2,4) mode is even larger than for the (0,12) mode, it follows that in 
the frequency range from the (0,4) cut-on frequency to the (0,12) cut-on fre- 
quency (a range of a one-to-three ratio in frequency) only the (0,4) mode will 
be excited over the first octave of the range, and only the (0,4) and (1,4) 
modes will be excited over the remainder of the range. Furthermore, relative 
amplitudes and phases of excitation of the (0,4) and (1,4) modes over the whole 
range can readily be estimated from the appropriate theoretical expressions. 

To suppress the (1,4) mode would require twice as many transducers, the 
additional transducers being placed on the wall- at x = a. 

It should be mentioned that the enhancement-suppression properties of the 
transducer array considered above are to a great extent insensitive to the 
sizes and shapes of. the transducers, and of the velocity distributions over 
their faces. Also, it can be shown that variation of the transducer velocity 
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over the face could occur without substantially altering the mode response, as 
long as the average normal velocity times the transducer face area remained 
approximately constant. In such a case, of course, the "acoustic center" of 
the transducer might not correspond to its geometric center. 

This brief consideration of duct excitation by an array of transducers 
shows that practical arrays can be designed and operated to give predictable, 
isolated generation of selected duct modes, for a relatively large range of 
modes in a given duct. The major limitation of this idealized modal excita- 
tion, in practice, is the variability of the acoustic velocity distribution of 
the individual transducers, over the degree of interaction isolation which 
can be maintained between individual units in the phased array and the 
local geometric variations in the duct at the mounting positions of the 
transducers. Each of these factors will cause variations from the ideal- 
ized situation discussed above with the primary effect being excitation of 
modes other than the desired excitation pattern. 


I I I . 3 . 4 . Sound pressure spatial control: experimental results 

After initial operating experience with the facility, it was deter- 
mined that interaction effects between individual transducers in the phased 
array were negligible, but that each transducer produced its own individual 
"signature" in the duct. The method of characterizing the output of each 
transducer is illustrated in figure 20. Each driver was excited individually 
and the microphone in the source plane, at the driver exit plane, was used to 
record the circumferential distribution of sound pressure and phase in that 
plane around the circumference of the source section as shown by the dashed 
line in figure 20. The method of recording the driver output was standard- 
ized. The microphone ci rcumferent i al traverse was always initiated at the 
same angular position, i.e., Q^O. The resulting sound pressure distribution 
is illustrated in figure 21 for driver number eight which is located at 157.5°. 
This driver characteristic curve was made at 876 Hertz (at 71°F), which is 
slightly above the angular mode cut-on frequency. 

The most obvious freature of the SPL distribution in the source plane as 
shown in figure 21 is the near field pressure behavior at the transducer exit 
plane near the duct wall. This near field is characterized by a peak in both 
the amplitude and phase. In this near field dominated region, the wave is 
propagating outward, thus the phase, change should be linear as the microphone 
traverses the circumferential contour. If the source is considered as a point 
source, then in the free field of that source the spherical wave generated 
will propagate outward, decaying at 6 dB per doubling of distance with a phase 
change of k per unit distance traversed. In this particular case, the phase 
change for the propagating wave is calculated to be 3.8 degrees phase/degree 
test section angle. The experimental behavior. for the point source phase 
change and amplitude decay are compared with the theoretical point source be- 
havior in figure 21. This illustration is particularly outstanding for 
demonstrating the point source influence. 




Figure 20. 


Illustration of Recording Technique for Driver "Signature 
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At some point around the duct circumference, the propagating wave merges 
with the modal structure in the duct which was excited by that source. In 
this particular example, which is just above the calculated (3,0) mode cut-on 
frequency, the primary response is in the (3,0) mode as can be seen by 
examining the sound pressure amplitude and phase response in the duct. The 
shape of the amplitude should be that of log (cos (30)) and the phase, for a 
standing wave, should have 180 degree changes at the node lines. This pattern 
is not strictly attained, as the phase shows the characteristic of a propa- 
gating wave in some regions (i.e. a linear increase with distance). Figure 22 
is included to show how closely matched the individual driver units are. 

Figure 22 shows that the characteristic for driver 2 is exceptionally similar 
to that of driver 8 shown in figure 21, with only very minor differences in 
amplitude and phase shape. Other drivers did have greater deviations, proba- 
bly due to duct dimensional differences as much as anything. Differences of 
the order of 1 to 2 decibels and 5 to 15° phase were noted in the character- 
istic shapes at a particular frequency at points away from the peak amplitude. 
While these differences away from the peak sound pressure under the driver may 
appear minor, they can give rise to significant errors when attempting a modal 
simulation. 

In the discussion in the preceding section, it was shown how, in the 
ideal situation, all modes other than the one desired, could be excluded by 
simple adjustment of the amplitude and phase of the piston type transducer as 
appropriate to the mode amplitude at the location of that transducer. For the 
example, which utilized four transducers (effectively), the amplitudes were 
Vo/2, -Vo, Vo, -Vo, Vo/2, where the end points were effectively 1/2 trans- 
ducers. If the transducers were identical, then the amplitude of each driver 
unit, directly under the exit of that unit, could simply be computed from the 
modal amplitude and phase desired at that location. An example of this is 
shown in figures 23-25. In figure 23, the amplitude and phase of each driver 
is shown, with phase indicated by a plus or minus sign. Figure 24 shows the 
response in the source section at 1000 Hertz with the idealized setup and in- 
dicates the location of each driver. The near field sound pressure and phase 
of the individual drivers has a predominant influence on the spatial distri- 
bution. However, these higher order spatial variations would tend to excite 
the (8,0) mode and this mode does not propagate at 1000 Hertz, since its cut- 
on frequency is of the order of 2200 Hertz. Figure 25 shows the wave as it 
propagated into the test section. Although the waveform and phase are very 
similar to the required log (cos0) and ±180 degrees shapes, there is signifi- 
cant deviation. Although not immediately obvious here, much of the deviation 
is caused by a spinning mode component of the same mode form. Other attempts 
at mode set-up were more effective, utilizing this simplified approach. 
However, in general, it was determined that the simple scheme was inadequate 
for satisfactory mode simulation over a broad frequency range. 

In order to obtain a more accurate modal simulation a more rigorous 
approach was used. In this method, all the details of sound pressure in the 
source section were incorporated. The more exact, mathematical, method for 
mode simulation is described below. 
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Figure 21. Characteristic Function for Driver Number 8 at 876 Hertz 



Figure 22. Characteristic Function for Driver Number 2 at 876 Hertz 
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(1 ,0) MODE SET-UP 


111 .7 


O.OdB 


in. 7 


SOUND PRESSURE 
AMPLITUDE 



Figure 23. (1,0) Mode Set-up for Idealized Driver Output 








First, the individual driver response, as shown in figure 21, was charac 
terized as an individual complex function. This function was digitally 
recorded during an angular sweep at 128 equally spaced locations around the 
circumference of the duct. Both dB amplitude and phase in degrees were re- 
corded. This function was then converted to a complex pressure for each of 
the driver units. With this information available, the following expression 
can be written to show the simulation of the desired angular function, f(0), 
by the complex driver response functions, (f>(0), 


16 

n=i Vn (6) = f(6) + E(6) (57) 


where A n can be represented as the complex amplitude of each transducer 
directly in front of the transducer outlet. A least squares approach for 
simulating the function, f (0) , is ideally suited to this requirement. The 
function, f(0) is transferred to the left hand side of the above equation and 
the resulting error, E(0), is minimized with respect to the complex amplitude, 
A , after squaring. The operations are 

16 

l Vn (0) ‘ f ( 0 > = E(e). (58) 

n= 1 


The modulus of the squared error function is 



16 


16 

G(e) = E ( 0 ) E*(8) = 

l A n * n (e)-f(e) 

. n= 1 


y a -f 5 '((|)) 

Ln= 1 J 


The unknown coefficients, A , are determined by the minimization indicated 
be low: n 


3G 3G 

3Re(A n ) °> a|m(A n ) °‘ 


( 60 ) 


These operations were performed with a digital computer analysis and the 
resulting amplitude and phase at the exit plane of each driver was numerically 
determined. 

The calculated sound pressure and phase for each driver unit for the case 
illustrated in figures 23"25 is shown in figure 26 , These are the driver out- 
puts determined as necessary to produce a wave of 120 decibels amplitude. 
Because of the differences in individual driver response, the amplitudes and 
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CALCULATED AMPLIFIER AND PHASE POT SETTINGS 
FOR THE FOLLOWING TEST CONDITIONS 


REQUIRED PRESSURE MODE - COS (2*PI* 1 *THETA) 
16 DRIVER RESPONSE FUNCTIONS USED IN FIT 
128 DATA POINTS PER DRIVER RESPONSE FUNCTION 
TEST FREQUENCY = 1000. HERTZ 


DRIVER 

DRIVER 

DRIVER 

NUMBER 

OUTPUT 

PHASE 


SPL-DB 

DEGREES 

1 

122.1 

302.4 

2 

122.4 

308.1 

3 

120.4 

308.3 

4 

114.8 

306.8 

5 

87.6 

214.1 

6 

115.0 

132.0 

7 

119.9 

131.8 

8 

122.6 

128.7 

9 

123.4 

130.4 

10 

122.2 

126.5 

11 

120.0 

129.3 

12 

114.6 

131.2 

13 

78.4 

1 7.7 

14 

115.2 

313.6 

15 

120.9 

311.7 

16 

121.3 

310.8 


Figure 26. Computer Output Listing Dr i ver Ampl i tude and Phase for "Least-Square" 
(1,0) Mode Set-up 
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phases computed in this fashion deviate significantly from those computed from 
the ideal model. The test section response for the more exact set-up is shown 
in figure 27 and should be compared with figure 25 to show the increase in 
fidelity. 

Finally, to show the results of the least square minimization procedure 
in predicting the sound pressure distribution at the source location, figure 
28 displays the measured and calculated sound pressure and phase distribution 
after determination of the complex amplitude function, A . The plot is repre- 
sented by 

16 

P( 0 ) = l A n <J> n ( 6l ) 

n- 1 


and is shown in dB amplitude and phase in degrees. The solid dots are the 
computed amplitude, with open triangles being the computed phase. The compar- 
ison with the measured angular distribution of these functions is quite close, 
as it was in all cases tested, indicating an insignificant interaction between 
transducers when all are operating. 

Thus, in summary for this part of the investigation, the objectives of 
the full-scale investigation have been outlined, the facility was described 
and the means for simulating angular modes in the annular duct facility was 
described and justified in detail. 


IV. SOUND PROPAGATION IN RIGID WALL ANNULAR DUCTS 


The primary objective of this series of tests was to provide basic infor- 
mation on cut-on frequency, modal phase speed, mode shapes, frequency response 
of the duct, and radiation impedance of and radiation characteristics from the 
unflanged termination. The first three items are described in this section 
and the latter two are described in sections VII and VIII. Effects of flow 
through the duct were also determined and are included in this discussion. 


The first sub-section will briefly present a description of the 
ments conducted, followed by a summary of analytical predictions for 
wall annular duct concerning cut-on frequency, phase speed, and mode 
Then the discussion will present experimental observations and where 
these will be compared with theory. 


exper i- 
the hard 
shapes . 
poss ib le, 


^9 


SOUND PRESSURE LEVEL - DECIBELS^ S0UN0 PRESSURE level in DECIBELS (re. 0.0002 dynes/cm 




Figure 28. Comparison of Calculated SPL and Phase D i st r ibu t i on at Source Plane for "Least-Square" 
(3,0) Mode Set-up at 876 Hertz 
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IV. 1. Description of Experiments 


The fundamental parameters in the modal theory of sound propagation in 
rigid wall tubes or ducts are the eigenvalues and corresponding eigenvectors. 
The eigenvalue, in the cases of infinite or finite length ducts, is either 
real or imaginary, with the frequency parameter at the transition point be- 
tween real and imaginary being the cut-on frequency. The first item selected 
for experimental evaluation was a determination of the set of cut-on 
frequencies for the annular duct. Two basic techniques were utilized in loca- 
ting the cut-on frequency. One involved experimental determination of the 
frequency of peak response amplitude in the region of cut-on frequency, for a 
single driver simulating a point source excitation. The other involved visual 
examination of the circumferential or radial sound pressure distribution in 
the test section, for a single driver excitation, to determine which frequency 
resulted in- a spatial response closest to the ideal modal pattern. 

Next, a technique for determining the phase speed, c^, was developed. The 
phase speed is inversely proportional to 0, i.e. c^ = c/6. In an infinite 
length duct (one in which there is no reflected wave), or a duct with a rea- 
sonably! high axial attenuation rate, determination of the phase speed is a 
relatively simple matter. It simply involves experimental determination of 
the phase difference between two points with the same radial and angular 
coordinates but with axial separation. .However, with both forward and back- 
ward traveling waves existing simultaneously, as occurs in the hard wall 
finite length duct, the amplitude of each wave must be determined before the 
phase speed can be measured. 

Finally, both angular and radial mode shape measurements were accom- 
plished for several modes at several frequencies. Although radial control of 
the sound pressure distribution within the duct was not possible, it was found 
that radial modes could be excited in relative isolation near the cut-on 
frequency of that mode. 

The following sections will contain more detailed descriptions of the 
experimental procedures developed during this phase of the testing. 


IV. 1.1. Determination of cut-on frequency 

In the early analytical treatments on sound propagation in annular 
ducts by Sofrin and Tyler (ref. 7) and Morfey (ref. 8), it was assumed that 
the source was a harmonic, distributed, pressure source. It was then shown 
that with a pressure source, the duct response resembles a unit step function 
in the region of cut-on frequency. Morfey went one step further and ran a 
simple experiment wherein he excited a inch diameter duct with a siren 
producing a ( 1 , 0 ) mode pattern. The measured response is shown in figure 29. 

Doak (ref. 37) conducted an analysis of a duct excited by a point veloc- 
ity source and computed modal responses. An example of the space integrated 
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Figure 29 * Cut-off Demonstration: Variation of SPL with Frequency 

(from Morfey, ref. 8 ) 
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mean square pressure response and the acoustic power for a source located at 
the center of one of the duct walls is given in figure 30. 

Since the infinite duct is effectively a finite duct with a pc termina- 
tion at the ends, an escape route would appear to exist for any resonant energy 
that might be built up in the duct, ft is perhaps surprising that the mode 
amplitude shown in figure 30 goes to infinity at the cut-on frequency, and is 
very large in the neighborhood of this frequency. Such infinite amplitudes 
would, of course, be expected in a completely enclosed, hard wall room at the 
normal mode resonant frequencies (see, e.g. Morse (ref. 38), Chapter VIII). 
Although this resonance effect in the duct mode case is perhaps less familiar, 
it nevertheless has been noted previously (e.g., Morfey (ref. 8), in the same 
reference in which figure 29 was presented, commented on the response to a 
point velocity source). 

A physical explanation of the behavior is given in reference 39 (see 
Chapter 2). In brief, it is that, in the neighborhood of cut-on, the mode 
wave motion is almost (completely so, at cut-on) a pure two-dimensional 
standing wave, just like a two-dimensional standing wave in a hard walled 
room, in fact. Therefore as no energy is required to maintain motion of such 
a standing wave in a lossless medium, the amplitude tends to infinity when the 
mode motion is forced by a volume velocity source at the cut-on frequency. At 
the same time, however, the impedance presented by the mode to the source also 
becomes infinite, so that to maintain, in practice, any finite volume velocity 
amplitude, 0. o , of the source, would require application of an infinite force at 
the cut-on frequency. Thus,' the back reaction of the mode on the source can 
be large at frequencies near the mode cut-on frequency. 

Therefore, in practice, the amplitude reached by a mode at the cut-on 
frequency is determined by two factors: not only the rather obvious factor of 

the wall absorption (neglected here), but also the effective internal imped- 
ance of the source mechanism. 

Morfey has shown theoretically that, for a constant pressure amplitude 

source (the other extreme from the constant volume velocity source used here), 

the mode frequency response factors corresponding to those shown in figure 30 

approach the asymptotic values of unity from below, as u mn /arK), without the 

resonant peaks at w /tn a 1 , as shown in figure 29. 

mn 

In practice, therefore, without specification of the source type, one can 
expect mode response behavior near cut-on to lie somewhere between the two ex- 
tremes: either a violent overshoot for sources behaving nearly like constant 

volume velocity sources as shown in figure 30* or asymptotic behavior for 
sources acting more like constant pressure sources, as illustrated in figure 
29. 

An example of the measured response in the annular duct facility due to a 
single driver excitation on the wall is shown in figure 31 (ref. AO). Based 
on the preceding discussion, it can be assumed that the source is a reasonable 
approximation to the point velocity source. In recording this response, the 
source pressure amplitude was held constant with frequency. It is obvious 
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DIMENSIONLESS POWER OR PRESSURE RATIO 


0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 
DIMENSIONLESS FREQUENCY PARAMETER - kR 


Figure 30. Dimensionless Acoustic Power and Space Integrated Mean 
Square Pressure for Source at Position x = a/2, y = b 
(see figure 19) (from Doak, ref. 37) 


FREQUENCY RATIO - f/f. 


Figure 31. Variation of Annular Duct Acoustic Response in the Vicinity 
of the First Angular Mode (1,0) Cufon Frequency 
(from Plumb lee and Dean, ref, 40) 







that, as expected, the response amplitude is not infinite. Damping mechanisms, 
not accounted for in the model, as discussed, limit the response amplitude and 
increase the frequency of maximum response. 

The analogy of the modal response to a point velocity source is that of a 
damped accelerometer with a base harmonic displacement. The physical analogy 
is shown in figure 32(a) and the response is illustrated in figure 32(b). From 
this analogy, two useful data analysis tools can be extracted. First, the 
ratio of the peak amplitude of response to the final "s teady-state" response 
amp 1 i tude i s 



( 62 ) 


On examination of figure 31, it is seen that the peak pressure amplitude is 
approximately 9 dB above the "steady-state" pressure amplitude at higher fre- 
quency (for constant pressure vs frequency excitation at the source). This 
measurable quantity can be used to determine an "effective" damping ratio, y, 
for the mode from the above relationship, Wi th a 9 dB ratio, X /X q = 2.82 
and y = 0 . 1 8 . The frequency at which peak response occurs is 


ft = 



(63) 


For a value of Y = 0 . 1 8 , a frequency ratio, ft, is determined to be 1.035. That 
is, the frequency of peak response, based on an experimentally determined 
damping ratio, is 3.5% greater than the calculated undamped natural frequency. 
As will be shown later, this is the exact ratio of experimentally determined 
peak frequency to calculated cut-on frequency, as shown by the dashed line in 
figure 31.- 

Another useful technique for identifying mode cut-on frequency was also 
adopted for this study. This method involved examination of the circumfer- 
ential sound pressure response in the test section, for a single driver 
excitation, at various frequencies above and below cut-on frequency. in 
figure 33, angular SPL response patterns in the test section at the outer wall 
are shown for three frequencies; one below peak response, one at peak response 
and one above peak response. In determining the cut-on frequency of a mode, 
once again it is assumed that the response of a mode at cut-on will signifi- 
cantly exceed the response of all other modes for a point velocity source 
excitation. In figure 33, this assumption is confirmed. The peak response 
is virtually all (2,0) mode while at 5 Hz below peak response the mode contri- 
bution is virtually all due to the (0,0) and (1,0) modes. At 5 Hz above 
cut-on the mode contribution is primarily (2,0) but the response amplitude 
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is significantly reduced from that at 589 Hz. Using. this criterion to deter- 
mine mode cut-on frequency, a series of circumferential SPL distributions was 
made for small increments in frequency and then a judgment was made, based on 
a visual balance of peak response amplitude and distortion of the mode shape, 
to determine the cut-on frequency. 


Both of the above techniques were used extensively during this program as 
a means of determining cut-on frequency. The latter technique, as described, 
involved more judgment on the part of the experimenter than the first method; 
however, it does give a more positive identification to a response peak. 


IV. 1.2. Method for experimentally determining phase speed 

The analytical form for the axial propagating wave is given as 


- i k3z 
e 

where 3 is the wavelength constant. Both the axial wavelength and the propa- 
gation speed are inversely proportional to 3> i«e. ^ m>n = ^/3 m>n and C(f>m,n = 
c/3 m n where m,n refers to the mode number and X is tfie free space wavelength 
for £ particular frequency, f, at a speed of sound, c. Thus, in the case of 
the finite length duct, wherein a standing wave pattern exists in the axial 
direction (assuming reflection occurs at the open end), the most straight 
forward technique for measuring 3 m n is to accurately measure the frequency, 
speed of sound, and the axial standing wave pattern for a mode in relative . 
isolation. This method was suggested by Morfey (ref. 8). 

The primary problem in implementing this technique was the lack of an 
axial traversing microphone. Thus, measurements were made at discrete points 
in the direction of the axial coordinate, by inserting a microphone in the 
microphone probe holes shown in figure 14, and recording an angular sweep of 
sound pressure distribution such as is shown in figure 27. These were 
recorded at each of 33 axial positions and the peak amplitude was determined 
from the angular SPL plots. 

A plot of the resulting axial standing wave amplitude and phase is shown 
in figure 3^ for the (0,0), (1 ,0) , (2,0), and (3,0) modes, A least square 
curve fit routine, using the assumed forward and backward wave system was used 
to determine the best representation of the composite wave. This is shown by 
the curves plotted through the data points. Examining the phase plots, it is 
observed that as the reflection coefficient becomes very small, as in the case 
of the (1,0) mode, the phase increases almost linearly with axial distance; 
but as the reflection coefficient approaches unity as for the (3,0) mode, the 
phase variation approaches the ± 1 80 degree square wave function characteristic 
of the pure standing wave. The phase change with distance can be used to 
determine 3, only when no reflected (or standing) wave component exists, as 
mentioned earl ier. 
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AXIAL DISTANCE - INCHES 

Figure 34, Measured Axial Variation of SPL and Phase for (0,0), (1,0), (2,0) and 
(3,0) Modes in Hard Wall Test Duct at 1000 Hertz 
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The propagation or wavelength constant 3 mn . ' s determined by measuring 
the distance between minima, and the relationship, 


= x/e 

mn mn 


i 

is used to determine 3 .. 

mn 


IV. 1.3. Mode shape measurement 

The technique for measuring a mode shape was relatively simple, after 
it was determined that the mode was excited in isolation. Techniques for 
assuring this were described under cut-on frequency determination (IV.l.l), 

At the frequency of peak response for a particular mode, the driver units 
were set up as described in the section on modal simulation and the angular 
and radial variations of sound pressure and phase were recorded at the 
appropriate stations in the test section. 


IV. 2. Analytical Predictions 


The first analysis for prediction of the sound field in a hard wall 
annular duct is that by Tyler and Sofrin (ref. 7)* Their analysis resulted in 
a tabulation of the characteristic numbers for the E functions, which are the 
radial mode shapes. The characteristic numbers, which are a function of the 
hub-tip ratio, are directly related to the cut-on frequencies. 

The numerical analysis, described in a following section on sound propaga 
tion in lined annular ducts degenerates to the Sofrin-Tyler analysis when 
presented with rigid wall boundary conditions and no mean flow. The results 
are solutions to the following set of equations: 


E (k 1 ) = j' (k* ) + 0*; a V (k 1 ) = 0 

m N m my m 

j 

E (ok 1 ) = j' (ok 1 ) + oJfV (ok 1 ) = 0 
m m my m 


(64) 

(65) 


where primes denote differentiation with respect to the argument, a is the 
hub-tip ratio, k‘ is the characteristic number and J and Y are Bessel functions 
of the first and second kind. The computer solutions also yield the propaga- 
tion constant, 3, which is a function of k'. 
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IV. 2.1. Mode shape predictions 


The mode shape in the hard wall annular duct is a combination of a 
sinusoidal function in the angular direction and the E function in the radial 
direction. Because of complications in the angular wave form which can be 
introduced by 'noise' in the set-up procedure, a single mode might be excited 
in isolation, but the resulting waveform might appear to be very different 
from the theoretical shape. Occasionally, a mode set-up will visually appear 
to contain a distortion in shape, evidenced by a low peak to minimum ratio. 

This is illustrated in figure 35. However, upon analysis of the waveform 

using the technique described in section VII on RADIATION IMPEDANCE, it is 
observed that this type of behavior is usually the result of combining a 

spinning mode and standing mode of the same mode number. An analysis which 

reduces the combined wave to its spinning and standing mode components was 
programmed on a digital computer. The results of this for the waveform shown 
in figure 35 are displayed in figure 36. It shows that the spinning wave is 
approximately 7 decibels below the standing wave, A somewhat trivial observa- 
tion is that the minimum level of the combined wave is equal to the spinning 
mode component, if no other modes are involved. 

A generalized complex Fourier computer analysis was written to determine 
the Fourier components of the angular sound pressure distribution. An example 
of the computer printout is shown in figure 37. The first line identifies the 
axial position in the test duct and the second line gives the necessary infor- 
mation for identifying the data run. Next, the results of the Fourier 
analysis are listed, giving the amplitude and phase angle for the cosine and 
sine components. Finally, these cosine and sine components are combined to 
give the amplitude and phase of the combined wave, the amplitude and phase of 
the standing and spinning wave components and In the last column, the angular 
location in the test section corresponding to the first minimum of the com- 
bined wave. In this particular example, the set-up at the source plane corres- 
ponded to an equal amplitude in the first four modes [(0,0), (1,0), (2,0), and 

(3.0) ]. The analysis shows that this first condition still exists with ±2 dB 
at station 36, which is the entrance to the test section. Furthermore, the 

(4.0) mode is 25 dB below the plane mode component, which Indicates the degree 
of isolation achieved in the mode simulation procedure. 


IV. 3. Comparison of Experiment with Theory 
IV. 3.1. Cut-on frequency 

Using the methods described in the previous section, the non- 
dimensional cut-on frequency, kR 0 , was measured for nine modes and is indica- 
ted by asterisks in table I. The computed eigenvalues are shown in the table 
for the (0,0) mode through the (6,2) mode. In all cases, the experimental 
value is higher than the computed value, with the percentage ranging from 3.6 
percent to 0.7 percent. The decrease in error is consistent with increasing 
mode number, indicating a correlated effect such as ascribed to damping in 
the simple one-dimensional model presented in the previous section. 
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SOUND PRESSURE LEVEL - DECIBELS 


TEST SECTION ANGLE - DEGREES 

Figure 35. Typical Angular Mode "Distortion 




J - 


STANDING WAVE SPL AND PHASE 


SPINNING WAVE SPL AND PHASE 


-COMBINED WAVE SPL AND PHASE 


TEST SECTION ANGLE - DEGREES 


Figure 36 . Components of Stationary and Spinning Modes Contributing to Angular Mode "Distort! 
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HOLE NUMBER 36 MODE SET (4 MODES, EQUAL AMPLITUDE) 


FREQ 

1000.0 HZ 

TEMP 72.0 DEG F 

PHASE 

REF 0 DEG 


COS AMP 

COS PHASE 

SIN AMP 

SIN PHASE 

MODE NO. 

DB 

DEG 

DB 

DEG 

0 

109.6 

207.1 



1 

107.3 

133.3 

109.6 

315.9 

2 

97.5 

193.2 

113.1 

173.8 

3 

113.2 

143.2 

106.9 

129.1 

4 

84.0 

330.4 

77.3 

146.6 



COMBINED 

WAVE 

SPINNING 

WAVE 

STANDI NG 

WAVE 



AMPLITUDE 

PHASE 

AMPLITUDE 

PHASE 

AMPLITUDE 

PHASE 1 

PHASE SHIFT 

NUMBER 

DB 

DEG 

DB 

DEG 

DB 

DEG 

DEG 

1 

111.64 

314.97 

78.43 

82.53 

111.45 

314.97 

127.56 

2 

113.20 

174.28 

87.84 

255.32 

112.72 

174.28 

81.05 

3 

114.07 

140.57 

93.68 

165.83 

113.20 

140.57 

25.26 

4 

84.85 

149.77 

52.94 

125.02 

85.07 

149.77 

155.25 


Figure 37* Computer Output Illustrating Fourier Spinning/Stationary Mode Analysis 
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TABLE I 


COMPARISON OF THEORETICAL AND EXPERIMENTAL HARD-WALLED 
ANNULAR DUCT EIGENVALUES (a = 0.566) 


X 




cut-on frequency, R Q - duct outer radius, c - speed of sound 


Radi a 1 
Modf» No. 

Circumferential Mode No., m 

n 


0 

1 

2 

3 

4 

5 

6 

0 


0.0 

m 

*2.5976 

2.5716 

*3.8550 

3.8268 

*5.0895 

5.0492 

*6.2783 

6.2351 

7.3860 

1 


7.3219 

7.4509 

*7.9022 

7.8282 

8.4284 

9.2173 

*10.2948 

10.1587 

*11.2913 
11 .2167 

2 

1 *1.5155 

14.5776 


15.0672 

*15.6355 

15.4855 

16.0108 

16.6358 


*Expe r i menta 1 ly determined value 


The computed values of 'equivalent' s i ngle-degree-of-f reedom damping 
ranged from 8% to 1 8% for the various modes. The computed ratio of pressure 
at damped resonance to the 'steady-state 1 pressure, well above cut-on, varied 
from 9 dB to 16 dB for the modes measured above. This ratio was not measured 
for all the modes and no consistent comparison such as that given for figure 
31 has been made. This comparison could only be made if the mode is excited 
in isolation. When a point source frequency response curve was used as a 
means of looking at cut-on behavior, usually more than one mode was excited. 


IV. 3.2. Axial propagation constant 

In order to determine the propagation constant, 3, a number of SPL 
and phase measurements were made as described in Section IV. 1.2 and illustra- 
ted in figure 34. The curves drawn through the data points in figure 34 are 
the result of a least square curve fit through the data points at the theo- 
retical value of 3. In some cases, it appeared that a longer or shorter wave- 
length, corresponding to different values of 3. (or more basic, different 
values of cut-on frequency), would provide a closer fit to the measured data. 
As a means of determining the accuracy of the curve fit, and the 'goodness' of 
the theoretical propagation constants, the 'standard error of the estimate' . 
(the square root of the 'sample variance 1 ) was computed. This computation is 
described in more detail in the section on RADIATION IMPEDANCE. 
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In attempting to determine the best value of 3 to fit the measured data, 
the curve-fit routine was used repetitively and incremental values of 3 were 
used. The error was examined for each computation and the 3 value which pro- 
duced the minimum error was taken as the optimum value. 

This procedure seemed inherently prone to the whims of experimental 
error, since one bad data point could result in a minimized value of error at 
a cut-on frequency (or wavelength) significantly removed from the actual value. 
After examining the results, it was found that the theoretical value of propa- 
gation constant gave the ’best' fit to the measured data points in most cases. 
The conclusion from this was that within minor error limits, the propagation 
constant is as computed, within experimental tolerance, and thus, the theory 
agrees exceptionally well with experiment. 


IV. 3.3. Mode shapes 

Sound pressure distributions in both the angular and radial coordi- 
nate directions were measured for several thousand combinations of the 
controlling variables, which were: (l) frequency, (2) transducer mode simula- 

tion, (3) axial position, and (4) radial position for angular plots and 
angular position for radial plots. In most cases, the data were plotted on an 
X-Y recorder and also recorded on digital tape. 

Several modes were studied extensively, in order to determine axial stand- 
ing wave ratios, propagation speed, mode purity and comparison with theory. It 
was difficult to determine how much of this data could usefully be presented. 
Finally, it was decided to present only one angular and/or one radial plot for 
each mode studied, exactly as that plot was recorded and to show the calcu- 
lated mode shape in comparison with that measured data. Mode shapes are 
presented for the (0,0), (1,0), (2,0), (3.0), (*4,0), (5,0), (2,1), (5,1), and 
(6,2) modes in figures 38 through 52. These modes were recorded at various 
frequencies as noted. 

Figure 38 illustrates the angular sound pressure distribution for a plane 
mode simulation at 1000 Hertz. There is a variation in the angular direction 
of 2 decibels and 15 degrees phase angle, indicating the presence of other 
modes at very low level. At the frequency of this test, three other higher 
order modes, i.e. the (1,0), (2,0), and (3,0) were also cut on. Thus, this 
plot demonstrates the ability to excite one mode in exclusion of other propa- 
gating modes. Figure 39 demonstrates the radial behavior of the plane mode 
and shows no discernible variation from the theoretical shape. 

Figures 40 through 47 demonstrate the purity of both angular and radial 
mode shapes for higher order angular modes in combination with the lowest 
order radial mode. The angular variation of the (1,0) mode is shown in figure 
40 , with the computed amplitude and phase shape represented by solid circles 
and triangles. The measured phase shape is relatively true to the theoretical 
shape except in the region of very low amplitude. This deviation is not sur- 
prising, since at this low amplitude, other spuriously excited modes and 
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SOUND PRESSURE LEVEL IN DECIBELS (re. 0.0002 DYNES/CM 2 ) 



0 20 1,0 60 80 100 120 140 160 180 200 220 24o 260 280 300 320 

TEST SECTION ANGLE IN DEGREES 

Figure k6. Comparison of Experimental and Theoretical SPL and Phase for Angular Hard Wall 
Shape at 1428 Hertz 
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Figure 47. Comparison of Experimental and Theoretical SPL and Phase for Radial Hard Wall (5,0) Mode Shape 
at 1428 Hertz 
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spinning mode components not evident in the amplitude plot can seriously 
affect the phase meter output. 

Figures 41 and 42 show the comparison between theory and experiment for 
the (2,0) angular and radial mode shapes. Note that the predicted 1 dB lower 
amplitude at the inner wall as compared to the outer wall for the radial shape 
is reproduced by the experiments with great fidelity. Comparisons with simi- 
lar accuracy are shown in the cases reproduced in Figures 43 through 47. 

Three higher order radial modes were studied. Results from these tests 
are illustrated In Figures 48 through 52. Figure 49 is, of course, similar to 
Figure 4 1 , however, it is excited at a frequency very near the (2,1) mode cut- 
on frequency. The radial behavior is shown in Figure 48. The comparison is 
quite good with the exception of a zero shift. This is, in all likelihood, 
due to a small content of (2,0) mode, which when summed with the (2,1) mode, 
would result in the observed zero shift. A similar observation can be made 
about the comparison between the (5 , 1 ) theoretical and experimental radial 
mode shapes. Finally, a radial mode with two minima, the (6,2) mode, is shown 
in Figure 52. Here, the comparison is exceptional, especially considering 
that this test was conducted at 3778 Hertz. 

in summary, it can be concluded that duct modes in a rigid wall annulus 
are quite accurately predicted by the theory, both in shape and cut-on 
frequency. 


IV. 3.4. Frequency response In hard wall duct 

It was shown earlier in this section that the cut-on frequency could 
be determined from the frequency of peak response and the ratio of peak to 
steady state response amplitude for a single mode, such as shown in figure 31- 
In order to obtain a curve such as shown in figure 31 » it was necessary to 
excite the mode of interest in isolation. This is very difficult to do, since 
the mode set-up is dependent on frequency, remembering that the drivers have 
different characteristic functions (illustrated by figure 21) as frequency 
varies. Thus, to conduct a frequency sweep of more than 30-50 Hertz range for 
a single mode excitation would be very tedious, since the mode set-up proce- 
dure would have to be repeated every 30~50 Hertz. 

However, another frequency response curve, which produces informative 
results, is obtained by driving the duct with a single transducer, simulating 
a point velocity source on the wall, while maintaining the sound pressure out- 
put of that transducer at a constant level with increasing frequency. At low 
frequencies, where only a few modes are cut on (and propagating) the cut-on 
characteristics of each mode can be examined. Figure 53 is an example of the 
hard wall duct frequency response from 280 Hertz to 3000 Hertz. The mode cut- 
on frequencies are indicated for those which are clearly identifiable. 

The peaks following the primary cut-on peak are a result of interference 
between multiple modes, due to the differing wavelengths (caused by differing 
propagation speeds). This is more clearly illustrated in figure 54 for 
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gure 51 • Compar 1 son of Experimental and Theoretical SPL and Phase for Radial Hard Wall (5»1) Mode Shape 
at 2335 Hertz * 


PHASE ANGLE 









SOUND PRESSURE LEVEL 



Figure 52. Comparison of Experimental and Theoretical SPL and Phase for Radial Hard Wall (6,2) Mode Shape 
at 3778 Hertz 
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another position in the duct. Only the first four modes are involved in this 
illustration. At the cut-on frequency of the (l ,0) mode the wavelength 
constant, 3 -|q is zero, thus the axial wavelength is infinite. As frequency 
increases, this wavelength decreases very rapidly with the wavelength of the 
plane mode remaining inversely proportional to frequency. The wavelength for 
the higher order modes, A mn , is equal to 


mn 


f vl " 




( 66 ) 


This relationship is illustrated in figure 55 for the plane and (1,0) modes, 
showing the axial wavelength differences near cut-on. 

A very simple theory is hypothesized to explain the response peaks 
following the cut-on response maximum value. Assume the duct extends to in- 
finity in both directions. Then the phase difference between two modes at 
z, excited by a point source at z = 0, is, assuming that the two modes are 
in phase at the source, 


<£ = (k i - k 2 )z. 


(67) 


It is seen that when this phase difference Is 0, 2 tt, Att , ...., reinforcement 

between the modes occurs and when the difference is tt, 3tt, 5 mr cancell a 

tion results. Thus, in the infinite hard wall duct, assuming no attenuation, 
the relationship, 


nTi = (kj - k 2 )z, 


( 68 ) 


will result in either a maximum or a minimum if n is even or odd. If k]_ is 
assumed to represent the plane mode and k 2 a higher order mode, then 


nir 


27rf 

c 


( 1 - 6 ), 


(69) 


recognizing that 3 = 




2 

, whe re 


f 

n 


Solving this relationship for the frequency 
of n, the result is 


s the cut-on frequency, 
f, which corresponds to a value 


f 



nc 


nc 

K 7 * 


(70) 
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AXIAL WAVELENGTH, X - FEET 
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For the (1,0) mode, which has a cut-on frequency of 293.4 Hertz at a temper 
ature of 66.9°F (corresponding to a speed of sound of 1124.9 fps), the fre- 
quencies for constructive and destructive interference at a point z = 5.30 
feet from the center of the source are given in the following table. 


TABLE I I 

FREQUENCIES FOR CONSTRUCTIVE AND DESTRUCTIVE 
INTERFERENCE 


n 


0 


f 


00 


(Hertz) 


1 

2 


458.8 

309.0 


3 

4 

5 

6 

7 

8 


294.4 corresponds 

313.8 corresponds 

346.3 

385.9 
429.2 
475.0 


to pos i t i ve 3 


to negat i ve 0 


Solutions in the table show that as n is increased, the frequencies first de- 
crease and then increase. The values of n corresponding to increasing values 
of frequency represent waves traveling in the direction opposite to the plane 
wave, but still in phase with the plane wave at z = 0. These solutions are 
meaningless in the infinite duct and misleading in the finite length duct 
unless the end had a rigid termination at a distance from the measuring point 
equal to the distance from the measuring point to the source (i.e., an image 
source). Thus, only three values have any meaning in this comparison. A 
minimum is predicted at 458.8 Hz, and it appears that this is reasonably 
reproduced by the experiment. The maximum predicted at 309 Hz is also in 
reasonable comparison with experiment as well as the minimum at 294 Hz. 

However, it is obvious that several more maxima and minima exist in the actual 
situation. The only possible explanation for these is end reflection, result- 
ing in waves propagating in both directions. An analysis similar to the one 
Just presented could similarly be performed, accounting for the waves reflected 
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from the termination and the phase shift introduced by the radiation impedance. 
No doubt with this added complication, frequencies of all the other maxima and 
minima could be predicted. 

I V . 3 . 5 • Flow effects on cut-on frequency and mode shape 

It has been shown that the duct mode cut-on frequency is a function 
of both the space averaged value of. mean flow through the duct as well as the 
wall shear layer (Mungur and Glddwell, ref. 3) • But for the range of velocity, 
boundary layer thickness and frequency encountered in the tests described, 
it would also be expected that the plug flow velocity effect derived by Morfey 
(ref. 20) would apply equal ly' wel 1 . Morfey showed that the mode cut-on fre- 
quencies are shifted by the flow as follows, 



v/i 


-M 2 k 


R 

mn o 


(71) 


The introduction of flow into the test facility usually resulted in a signi- 
ficant drop in temperature, of the order of 10 to 20 degrees Farenheit. This, 
of course, resulted in a very significant change in cut-on frequency. However, 
after temperature was stabilized at the test velocity, a new cut-on frequency 
was determined by one of the methods described in section IV. 1.1. After the 
test, the results were converted to dimensionless frequency, kR Q . 

As described in section IV. 1.1, the location of the exact cut-on fre- 
quency relies, to a great extent, on the judgment of the experimenter. During 
this series of tests, a different investigator was conducting the cut-on fre- 
quency determinations. Because of this, the values of cut-on frequency at no- 
flow were higher than determined in the previous no-flow cut-on frequency tests 
and reported earlier in this section. The cut-on frequency test results from 
both series of tests are shown in figure 56, for the (1 ,0) , (2,0) , (3,0), 

(5,0), (2,1), and (5,1) modes,” The lower curve for each mode Is a plot of the 
above cut-on relationship, based on the theoretical value of no- flow cut-on 
frequency. This shape is then transferred for comparison with measured data. 

The trends displayed by the measured data certainly confirms the rela- 
tionship presented by Morfey. While there is some scatter in the data, it 
must be realized that this scatter about the theoretical shape is less than 
0,5% for all data recorded and except for 2 data points (the no-flow data for 
the (2,1) and (5,1) mod es) , t he scatter is less than 0.2% and generally 
•follows the predicted ^1-M 2 shape. 

Theoretically there is no effect of flow on mode shape in the hard wall 
duct- if the plug flow type of analysis, such as Morfey used, is chosen. There 
is, however, the possibility of a significant effect due to the wall shear 
layer. The numerical analysis described in the next section was used to com- 
pute mode shapes for the experimental flow velocity profiles. The modification 
of the mode shape due to the very thin shear layer for this facility was 
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computed to be insignificant. Examples of experimentally determined mode 
shapes for the (3,0) angular mode and the (5,1) radial mode at various veloci- 
ties are shown in figures 57 and 58. The degradation of shape with increasing 
velocity is due to flow noise, either present in the duct or else generated by 
the microphone probe. The data, as presented, were filtered through a 2 Hertz 
narrow band filter. The shapes in figures 57 and 58 are the plots recorded at 
the experimentally determined value of cut-on frequency. The radial shapes in 
figure 57 confirm that the effect of flow on the shape is insignificant, if 
the flow noise effects are ignored. 


V. WALL IMPEDANCE MEASUREMENTS 


The prime requi remen t , as defined in section VI for the accurate predic- 
tion of acoustic attenuation in a lined duct is an accurate determination of 
the boundary conditions (in this case, the acoustic impedance) at the wall. 

The following method of wall acoustic impedance measurement, commonly 
termed the “two microphone" method has been developed using as a rationale the 
fact that impedance values measured in the duct environment would be more 
representative than values determined from any other technique. That is to 
say, other types of impedance measurements (e.g., from an impedance tube test) 
have been made in the past and found to be inadequate for the duct application, 
particularly when mean flow or peculiar duct geometries are used in the atten- 
uation tests. 

The two microphone method was first suggested by Sivian (ref. 42), but 
for some reason was not considered practical until recent times when Philips 
(ref. 44) first attempted to use the method (although in his case the position 
of the measurement locations has raised some questions on the validity of his 
results). Soon after this, this method was used by Binek. (ref. 45) with some 
success, using a turbulent boundary layer as the source pressure excitation. 

In this report the theoretical justifications (together wi th the implied 
limitations) follows that given by Tester and Dean (ref. 46) and will only 
be briefly summarized here. The method is based on the fact that a simple 
relationship between acoustic pressure and particle velocity exists. It is 
assumed that the particle velocities on either side of the cavity surface 
layer (considered to be infinitely thin) are identical.* The condition of 
application for this method is that the liner is locally reacting, i.e. that 
the wavelength is large compared to the cavity cross dimensions and that there 
is no transmission through the lining in a direction parallel to the liner 
face (rigid cavity walls). Only two measured quantities are required for the 
two microphone method: (1) the acoustic pressure at any point in the cavity, 

(2) the acoustic pressure above the cavity surface and its phase relationship 
with the cavity pressure. 


*This discussion particularly concerns a resonant cavity type of liner. 
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SOUND PRESSURE 10 dB 



(a) Angular SPL Distribution 
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(b) Angular Phase Distribution 

Figure 57- Effect of Flow Velocity on Angular Hard Wall (3,0) Mode Shape 
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(a) Radial SPL Distribution 



Figure 58 . Effect of Flow Velocity on Radial Hard Wall ( 5 , 1 ) Mode Shape 
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V.l. Basic Theory of Local Impedance Measurement “ 
Two Microphone Method 


Consider a lining as in Figure 59 with a y-axis origin at the back wall. 
The liner consists of a porous surface and an air space backing with trans- 
verse dimensions less than a half-wavelength; i.e., wave motion is in the y 
di rect ion on ly . 


RESISTIVE 



Figure 59. Classical Resonant Cavity Liner Element 


The acoustic pressure in the cavity is given by, 

P = 2 P B e Jll,t cos(ky), 

where P D is the cavity back wall acoustic pressure* but 

D 

P |*jr = (the linearized momentum equation). 

This results in the following expression for particle velocity, u, 

u = 2 sin ky # 

PC o 


(72) 


(73) 


( 7*0 
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The impedance, Z, is the ratio of pressure to particle velocity, 


1 = ' J cot(ky) ' 


with Z being the Impedance at any cavity point y behind the resistive facing, 
normal i zed to pc. 

The normalized acoustic impedance at the front surface of the resistive 
facing is 


P F -jP F e J * 

Z = — = -r cosec (k 1 ) , 

U P B 


( 76 ) 


where <j) Is phase angle between Pp and P^. 

In terms of sound pressure level measurements, the impedance is 


Z = -j 10^ SPL F" SPL B^ /20 e-* <t> cosec (kl ) . 


( 77 ) 


The resistive and reactive components are 


R = 10 


(SPL C - SPLj/20 sinj 


B 


in (kl) 


w _ _(SPU - SPLJ/20 cos <J> 

X - 10 P B sTnTkTT 


( 78 ) 

( 79 ) 


V.2. Measurement Error Effects 


V.2.1. General considerations 

There are two major disadvantages associated with the two microphone 
method. First, the measurement position of the microphone in front of the 
liner must be close to the liner surface; that is, at a distance small com- 
pared with the wavelength (otherwise the amplitude and phase measurements will 
be in error, depending of course on the properties of the duct acoustics 
field). But it must not be so close that an unrepresentative pressure is 
measured; for example, if the liner is a simple perforate, then close to the 
surface the pressure adjacent to a rigid part of the surface will differ from 
that adjacent to a hole. Thus, the effective measuring point must be in the 
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“far field" of the liner but still at a distance small compared with the wave- 
length. Thus, it is even possible that this proper measurement point may not 
exist. 

The second disadvantage concerns a similar difficulty} that is, if this 
distance from the measurement point to the surface is comparable with the ^ 
boundary layer thickness, then the basic assumptions of the. 'two microphone' 
method are invalidated and the impedance estimate will, strict ly. speaking be 
incorrect. In addition, the measurement of an acoustic signal within the 
boundary layer becomes difficult in terms of separating the acoustic signal 
from the other pressure components of the turbulence in the boundary layer. It 
should also be noted that the two microphone method can strictly be applied, 
only to 'simple lumped systems' adhering to the laws of homogeneous isotropic 
medi a. 


V.2.2. Error effects on two-microphone method evaluations 

For the expected values of resistance, R, with a nominal ^0 rayl (felt- 
metal) liner with a 1/2 " cavity depth, the changes in sound pressure level 
between the front and back walls of the cavity are very small, particularly at 
the low end of the frequency range, thus considerable care in measurement is 
required. The necessary increase in accuracy was accomplished by subtracting 
the two signals (filtered through a 2 Hz filter) with an analog computer. and 
then amplifying the difference signal. The result was recorded from a digital 
voltmeter to an accuracy of ± 0.03 dB. 

The phase differences between each probe microphone are small. in the low 
frequency region and thus the resistance (proportional to sin A<j>) is very sen- 
sitive to phase inaccuracies (the meter used has, for all practical purposes, 
an accuracy to 0.6 degrees). From the theory presented earlier the other. two 
major variables are temperature (which controls the speed of sound) and micro- 
phone position* from the back wall. Temperature is capable of being measured 
to an accuracy of 0.2 degrees F and microphone position (measured by a 
precision potentiometer) to better than 0.005 inch. 

Figure 60 shows the percentage errors of resistance and reactance versus 
frequency for a 10% variation in each of the above variables for the *t0 rayl 
Brunsmet liner with a 1/2" cavity depth. The sign interactions for these 
changes are given below in Table III from which it can be seen that positive 
changes in any error quantity lead to positive resistance errors and negative 
reactance errors. 


*i.e., this equals the cavity depth for the two microphone evaluation. 
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Figure 60. Effect of 10% Variation in SPL, Phase Angie, Temperature and Cavity 
Depth on a Nominal 40 Rayl Liner with 1/2" Cavity Depth 
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TABLE I I I 


SIGN CHANGES OF ERRORS IN IMPEDANCE 



Sign of Error 

Sign Change 

of Impedance 

Error Quantity 

Quantity 

R 

X 

Sound Pressure Level (d B) 

+ 

+ 

- 


- 

- 

+ 

Phase Angle (Degrees) 

+ 

+ 

- 


- 

- 

+ 

Temperature (°F) 

+ 

+ 

- 


- 

- 

+ 

Cavity Depth 

+ 

+ 

- 


— 

— 

+ 


The effect of a 10% phase variation can be seen to cause a change of 
almost 10% in resistance for all frequencies up to about 2 kHz with the error 
reducing to about 7i% at A kHz. The effect of this change on reactance, on 
the other hand, is very small at low frequencies, rising to about 2% at 2 kHz 
and then sharply rising to 10% at A kHz. 

A 10% change in the sound pressure level difference results in a slowly 
increasing frequency characteristic for both resistance and reactance, pro- 
ducing in each case changes of only 2i% at A kHz. The combined effect of 10% 
deviations in both SPL and phase differences was also estimated and can be 
seen to reduce the variations in resistance and to i ncrease the reactance 
changes compared to the 10% phase only changes. 

The variation of temperature affects only the speed of sound. A 10% 
change In degrees Fahrenheit would not be expected to have much effect. This 
can be seen to be less than 1% with an increase in both resistance and reac- 
tance noted. 

On the other hand, a 10% variation in cavity depth produces a peaked re- 
sistance characte rist i c (of about at 1000 Hz) and a decreasing reactance 
characteristic with frequency (from 10% at 300 Hz to 6% at A kHz). Of the 
four variables it appears that cavity depth ( or microphone position) is a 
most important parameter being capable of affecting both the resistance and 
reactance roughly in proportion to its own variability. The conclusion is 
therefore, that precise cavity construction is necessary to produce the 
desired impedance values or conversely, adjacent cavities of different depths 
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will cause proportionate changes in both resistance and reactance, which could 
be an effective means of constructing a broad-band liner. 


V.3. Experimental Configurations 


Two distinct two-microphone measurement configurations were used, depend- 
ing on the two different cavity depths, viz 1/2 inch and 1-1/2 inches, respec- 
tively. For the 1/2" depth cavity liners figures 61 and 62 show the layout of 
the probes in the test duct. The microphone voltages were filtered through a 
2 Hz filter, converted to logarithmic DC voltages (as described earlier in the 
facility description) and then input to the analog computer circuit of figure 
63. This produced an amplified difference signal which was measured by a 
digital voltmeter (or by a probability density analysis on a real-time corre- 
lator, if fluctuations caused by low signal to noise ratios existed when high 
air flows were present). The phase angle was measured between each microphone 
with a digital phase meter, using the filtered amplified microphone outputs 
di rect ly . 

Calibration of the microphone and amplifier amplitude and phase charac- 
teristics were achieved by placing the two microphones at a point where they 
would be measuring the same acoustic field. For convenience, this point was 
chosen as the back wall of the liner cavity. (The liners can be considered to 
be locally reacting, thus the acoustic field is normal to the back wall of the 
cavi ty . ) 

For each chosen frequency the procedure is as follows. First, the front 
wall microphone is set such that the probe tip is at the same point as the 
fixed back wall microphone. Then, the amplitude and phase calibrations are 
carried out. Next, the front wall microphone is moved until it is just flush 
with the inner (duct) surface of the liner. The amplitude and phase differ- 
ences are again measured and the differences between these values and the 
calibration values are used to evaluate the wall impedance. 

This method was used on both the 20 Rayl and 40 Rayl liners at 6 differ- 
ent hole positions, i.e. hole numbers 16 , 18, 22, 24, 28, 30 (the hole number 
represents the distance in inches from the exit plane of the duct liner). 

Three values of mean flow velocity were tested, i.e., zero, 50 ft/sec, 150 ft/ 
sec. 


For the 1-1/2" cavity depth liner a slight modification of the fixed two 
microphone technique described above was used. One microphone was traversed 
from the back wall to the front wall and out Into the duct and sound pressure 
amplitude and phase were recorded directly with an X-Y recorder. The wall 
pressure amplitude and phase values (phase was measured relative to an oscil- 
lator) were taken directly from this traverse by reading from the plots at the 
appropriate location. In order to eliminate any probe end effects, the probe 
tube was modified, as shown In figure 64 such that the liner surface hole was 
always filled with the probe tube. The sound pressure was detected by two 
small holes drilled in the probe tube at a right angle to the flow. 
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Figure 61 . Schemat i c and Photograph of Two- M i crophone Test Arrangement 



Figure 62. Overall View of Two-Microphone Test Installation 







Figure 63. Analog Computer Circuit for Accurate Measurement of Two-Microphone 
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Figure Sk. Travers i ng "One-Microphone" Modi fi cation of Two-Microphone Impedance 
Measur i ng Method 









This procedure was used for this liner as the honeycomb core (necessary 
to satisfy the local reactivity assumption) had a transverse dimension of 1/8 
inch (little more than twice the probe diameter), making the installation of a 
second microphone in the back wall of the same cavity exceeding difficult, 
without communicating with adjacent cavities. Figures 65 and 66 illustrate 
some typical pressure and phase traverses at different frequencies. The theo- 
retical (or "classical") behavior of the sound pressure in the cavity can be 
demonstrated by a comparison with a plot of 20 log (cos kl) superimposed on 
the curves (with an arbitrary vertical displacement for clarity). In all 
cases up to about 0.1 inch from the liner surface the theoretical and measured 
curves are parallel thus implying that the "lumped" impedance concept of the 
liner surface is valid with quite well defined edges and that the liner is 
essentially locally reactive. In figure 66 the 180 degree phase change within 
the cavity, expected in a standing wave situation, is shown particularly well, 
the pressure minimum being outside the dynamic range of the recorder thus 
producing a "flat" on the pressure traverse wave at this point. 

The determination of the back wall and front wall points which define the 
values of pressure and phase to be used require careful consideration. 

The back wall situation is relatively simple: as the probe tube is moved 

in, the measurement holes are progressively exposed, thus the pressure level 
very rapidly rises, then, after complete exposure, the cavity level decreases 
with distance, initially at a very slow rate (in fact at 1 = zero the gradient 
of pressure, proportional to sin(kl), is zero). Thus, the back wall pressure 
can be estimated from the peak level which will occur approximately within one 
half of a measurement hole diameter from the geometric back wall point (i.e. 
within 0.010 to 0,020 inches). As an additional check the phase changes are 
also very rapid until the measurement holes are fully exposed after which 
point the phase will level off (unless some cavity dissipation is present due 
to wall flexure or bond leakage). The front wall situation is slightly more 
difficult. The obvious starting point is the geometric front surface, 
however, immediately prior to that point, the pressure and phase traverses 
through the resistive layer will almost invariably show a high rate of change, 
thus, after the measurement holes are exposed on the duct side of the surface, 
a pronounced "knee" occurs followed by a finite pressure and phase gradient. 
The point of measurement is thus taken to be that termination point of the 
"knee" just before the commencement of the relatively shallow pressure and 
phase gradients on the duct side of the surface. Reference to figures 65 and 
66 show this point to be approximately 0.010 to 0.020 inches away from the 
surface. 

In the cases when mean flow is present this point is subject to a little 
more uncertainty (see figure 67) due to the fact that the boundary layer 
produces some distortion and turbulence induced fluctuations of pressure and 
phase on the duct side of the surface. However, the "knee" point is still 
quite obvious and can be used to estimate accurately the surface pressure 
level and phase, the level differences for each flow case being quite marked 
(see figure 67 for the zero flow and 50 ft/sec. comparison at 1000 Hz and a 
3000 Hz, 50 ft/sec. case for comparison with the zero flow case shown on 
figure 66). This "knee" estimation problem becomes more difficult in those 
cases where small standing wave pressure gradients occur near the surface, 
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Figure 65 . Effect of Frequency on SPL and Phase Angle Traverse Measurements for 
Welded Fibrous Metal Liner with 1 1/2 Inch Cavity Depth (300 - 2000 Hz) 
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i.e. at low frequencies where the wavelength is more than ten times the cavity 
depth or at higher frequencies where the wavelength is approximately equal to 
cavity depth and the standing wave pressure minimum is near the middle of the 
cavity. In these cases the pressure and phase charges through the liner 
surface are small and the knee points are difficult to determine accurately. 

It should be emphasized that in these cases 2 Hz filtering is necessary in 
order to reduce the turbulent flow distortions. 

Figure 68 is included as an example of the traverse variations encoun- 
tered. The first case exhibited the problem of drilling holes in the liner 
for the probe. The traverse shown is typical of a torn cavity side wall. The 
remaining traverses show evidence (as variations from hole position to hole 
position of pressure level and phase) of surface non-uniformities and/or seal- 
ing imperfections between the honeycomb edges and the surface resistive layer. 

Some variations between hole positions will always occur due to experi- 
mental error so the individual hole measurement must be averaged to produce a 
representative answer. However, it should be emphasized that an "average 11 
impedance will only represent the "average" properties of the liner if the 
variability of the measured wall impedance is of small scale and spatially 
random in character. Any mean axial trends or variations should be looked 
for, particularly as these could produce marked changes in the modal propaga- 
tion constants associated with the acoustic field. The acoustic properties of 
commercially available materials commonly exhibit large scale 'patchiness 1 
with resistance variations possibly even greater than ± 25 %, thus any quoted 
value should be treated cautiously. 

In some high frequency tests of the 1/2" cavity depth liners at high 
flows, for holes far downstream, the signal to noise ratio was so low that 
peculiar results (originating from turbulent pressures) appeared; these were 
ignored in the averaging process. In one case at hole 16 of the 20 Rayl 
liner, all the values were in considerable disagreement. This test was re- 
peated and the same answers resulted, the implication being that, at this 
point, there exists a gross manufacturing defect in the surface resistive 
layer. Figure 69 shows examples of imperfections such as large scale voids 
and surface warping. Results from these regions were not included in the 
averaging process since they were not considered as representative of the 
important liner region further upstream. 


V.^. Results of Measurements 
V.4.1. General comments 

Five different linings were tested, three of the resonant cavity type 
(locally reacting) and two of a non-local ly reacting porous homogeneous 
plastic foam. Two of the resonant cavity linings were constructed to have 1/2 
inch cavity depth with felt-metal resistive surface facings of, nominally, 20 
Rayls and ^0 Rayls flow resistance, respectively. The third sample of this 
type had a cavity depth of H inches with a surface resistance layer of welded 
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random fibers. The plastic foam liners were both 2 inches in thickness, from 
two different companies, one of them being a controlled manufacture acoustical 
foam. 

The two microphone method of impedance evaluation was used for all the 
resonant cavity type liners, while the standard impedance tube method was used 
for the homogeneous plastic foam liners (this method was also used as a check 
on the zero flow values of the resonant cavity type liners). For the 20 Rayl 
and 40 Rayl liners with 1/2 inch cavity depth, both the impedance results 
(from the two microphone method) expressed as a mean over all measurement 
holes, and the standard deviation, can be seen in figures 70-75. Two examples 
of different hole impedance curves (representing a 'bad' measurement case and 
a 'good 1 case, respectively) are also presented (figures 70 and 77) and it can 
be seen from these that local variations of material properties are present, 
particularly in the resistance. This is to be expected, considering the sen- 
sitivity exhibited by this parameter in the error analysis. The averages over 
all holes for these two liners are presented in figures 78 and 79 as the zero 
flow case. 

The reactance can be seen to be far more stable, probably due to the uni- 
formity of the cavity depth between the various test positions. The reactance 
consists of a negative stiffness (the cot(kl) term) plus a positive mass term 
associated with the liner surface. This latter term is very small for the felted 
metal materials (shown by the difference between the measured total reactance 
and cot(kl) curves) and is not consistent between holes due to non- 
uniformities of the resistive surface material. 

In addition to the in-duct measurements, h flat samples of the lining re- 
sistive surface layer bonded to plastic honeycomb were evaluated in an imped- 
ance tube. These measured values should be equivalent to the "in-duct" 

(figures 78, 79) measurements at zero flow and from figures 80 and 81, good 
agreement can be seen (considering the variations possible in the surface 
layers) in all the samples tested. However, in the welded fiber metal 
sample sheet used on the 1-1/2" facing, the derived impedance, as shown in 
figures 82 and 83, appears to be considerably larger (of the order of 0.2 
or greater) and frequency dependent, with an increase at the low and high 
end of the frequency range, especially for the mass reactive terms. 

Considerable variability of resistance was also experienced, as the small 
cell size of the honeycomb used was of the same order of magnitude as the 
surface density striations of this material. Thus some honeycomb cavities, at 
the points of measurement, were covered by a high flow resistive facing and 
others by a low flow resistive facing. 


V. 4. 2 . Effect of material property variations 

In many of the impedance tests, a standard deviation of the impedance 
difference between various measurement positions was computed. In all cases 
the standard deviation observed for both resistance and reactance decreased 
with increasing frequency. These deviations incorporate both material 
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Figure 82. Comparison of Impedance Tube and Two-Microphone Impedance 

Measurements for Liner 1 (Welded Fibrous Metal - 1 1/2 Inch 
Cavity Depth) Without Flow 


Figure 83 . Two-Microphone Wall Impedance Measurement for Liner 1 (Welded Fibrous 
Metal - 1 1/2 Inch Cavity Depth) at 50 fps Mean Flow Velocity 
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property variations and measurement errors, which from the earlier evaluations 
are most likely to be caused by phase errors and/or cavity depth errors for 
the resistance deviations, and cavity depth errors for reactance deviations. 
However, these will be quite small (of order of 5 % at low frequencies) in 
comparison to the observed low frequency deviations of resistance of 20-25% in 
the 20 Rayl liner zero flow case, and 20-30% in the ^0 Rayl liner case. Thus, 
it appears reasonable to ascribe these deviations to gross changes of the 
porous surface material properties. 

The effect of a flow velocity of 50 feet per second over the surface is 
to increase the resistance deviations in the low frequency region, particularly 
with the 20 Rayl liner where a 5% or more increase can be seen by comparing 
figures 71 and 72. In contrast, the reactance deviations are little affected, 
thus implying, from the earlier analysis, that phase changes are the cause of 
this increase. It is also interesting to note that when the mean flow veloci- 
ty is increased to 150 feet per second the low frequency deviations are 
increased but not in proportion to the flow velocity. However, the magnitude 
of the resistance increased as well, thus the percentage deviations remained 
constant or were even reduced. 

The bO Rayl felted-metal liner resistance deviations remained largely 
unaffected with flow velocity, although the increase in mean values produced a 
reduction in percentage error. This implies (from the assumption that bulk 
flow properties in the duct are not affected by the liner change) that the 
variations observed with the 20 Rayl liner could be flow interaction with the 
materials, this interaction having varying degrees of effect depending on the 
local variations near the test holes. Again, with this liner as with the 20 
Rayl liner the reactance variations remained essentially independent of flow 
ve loci ty . 

The welded fiber metal liner was not as extensively tested as the 20 Rayl 
and bO Rayl felted metal liners; however, its mode of construction produced 
striations of varying flow resistance, which can be seen as a cyclic variation 
in resistance, as shown in the multiple hole test data at 1000 Hz in Figure 8 A. 

The impedance of the plastic foam linings was determined from samples cut 
from a large sheet and tested in an impedance tube by the standing wave method. 
The results for the two materials are given in Figures 85-86. As can be seen 
from Table IV, one of the two foams exhibited considerable variations between 
supposedly identical test samples at 1000 Hz. This is unusual for "acoustic" 
foams and as such was considered undesirable test material. The other foam 
did not exhibit such variations. 1 

The impedance of the controlled manufacture acoustic foam is shown in 
figure 86* This material, which has the same peaked characteristic in 
impedance, did not exhibit the impedance variability of the foam previously 
di scussed . 
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igure 84. Effect of Axial Measurement Position on Wall Impedance for Liner 1 
(Welded Fibrous Metal - 1 1/2 Inch Cavity Depth) Without Flow.f = 1000 Hz 











TABLE IV 


IMPEDANCE VARIATIONS OF SIX SAMPLES OF 
2" THICK POLYURETHANE FOAM AT 1000 Hz 


SAMPLE 

RESISTANCE R. 

REACTANCE X 

1 

.755 

.056 

2 

.983 

.168 

3 

.712 

,0A2 

A 

.809 

. 034 

5 

.863 

.086 

6 

0 

00 

• 

-.078 


V.A.3. Mean values and related flow effects 


Figures 78 and 79 summarized the effects of flow on the 20 Rayl and AO 
Rayl liners. A marked frequency effect exists in the resistance values in 
both cases; this being a preferential increase in the low frequency region. 

The 20 Rayl case appears remarkable at 300 Hz as this increase is almost 100% 
and then tapers off to 20% at 1000 Hz. It should be remembered, however, that 
percentage increases in this case are perhaps not particularly meaningful. On 
considering the A0 Rayl case (a material of similar construction), almost the 
same magnitude increases (with flow at various frequencies) can be seen, 
implying that the increase in flow induced resistance can be considered as an 
additional "lumped" element, dependent only on the mean flow. 

The effect of flow on reactance appears small and of a positive nature 
and is often difficult to estimate as it is of the order of the natural 
material reactance variations. However, it is possible to conclude that the 
mass reactance of the liner increases with velocity, this increase then de- 
creasing in magnitude with increasing frequency. The cavity reactance term 
(cot(kl)) is also shown, and as expected dominates the total reactance except, 
of course, near the resonance frequency. It should be noted that these obser- 
vations apply only to the 1/2" liners where all frequencies of interest are 
be low i resonance and thus the cavity reactance term, by definition, is negative 
in sign. 

The welded fiber metal liner was not studied as comprehensively as the 
other two resonant cavity types, thus work concentrated only on a few frequen- 
cies and flows up to 50 feet per second. However, from a comparison of figures 
82 and 83 , the resistance increases little (if any) and the reactance curves 
shift slightly to the left (implying a higher mass term). The drop in reson- 
ance frequency (where reactance is zero), compared to the 1/2" liners, is due 
to the increased cavity depth of 1-1/2" for this liner. 
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V.5. Conclusions 


From the theoretical discussion and the results presented earlier, the 
following comments can be made. First, the in-situ impedance measurement 
technique is a viable one, the success of its implementation being governed by 
an appreciation of the effects of errors, and a regard for the necessary equip- 
ment for minimizing these errors. 

Secondly, the importance of uniformity in the construction of the acoustic 
surface liners has been made clear. The variability of the properties of the 
surface layers is currently too large for preci se classification of the wall 
impedance, thus implying extreme difficulty in making precise attenuation 
predictions. Large sample area tests (such as in an impedance tube) serve only 
to "smear out" the local variations which influence any grazing acoustic field 
application, such as duct linings. However, for some practical situations the 
"average" wall impedance, as presented, could be "adequate" for design 
purposes. 

Finally, the effects of mean flow over the surface of the liners tested 
appears to increase the acoustic resistance (with a low frequency bias) and to 
reduce the reactance slightly. The quantifying of these effects is a prereq- 
uisite for the successful prediction of the attenuation of the sound field in 
a duct. 


VI. SOUND PROPAGATION IN LINED ANNULAR DUCTS 


The primary objective of this series of experiments was to provide a suf- 
ficient amount of data to validate the theoretical predictions for sound pro- 
pagation and attenuation in an acoustically lined annular duct sustaining a 
radially non-uniform mean flow. This involves validation of both the eigen- 
values (which consist of the propagation constant, 3 » and the attenuation 
constant, a) and the eigenvectors (or natural mode shapes). In addition, other 
factors which influence the eigenvalues and eigenvectors such as the wall 
impedance boundary conditions, the termination radiation impedance, the sound 
pressure distribution at the source plane, and the mean flow and temperature 
distributions, must be experimentally determined in order to completely define 
the experiment conditions. 

The experiments conducted in the facility with sound attenuating liners 
provide the nucleus for comparison of annular duct theoretical predictions 
with measured data. Experiments were conducted with five different sets of 
liners. The parameters varied for each liner configuration were: (1) angular 

pressure distribution, (2) frequency, and (3) flow velocity. 

The pressure distribution studies are divided into two regimes; the first 
being below the first higher order radial mode (0,1) cut-on frequency and the 
second being equal to or greater than that frequency. At frequencies below 
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the first higher order radial model cut-on in the hard wall source section 
(1620 Hertz), the radial pressure distribution at the entrance to the lined 
duct is determined by the angular pressure distribution, if it is assumed that 
higher order mode reflections at the interface between the lined and unlined 
sections are negligible. From experimental observations, this assumption 
seems to be reasonable as long as no higher order radial modes will propagate. 
Thus, because of the control which can be exercised over the angular pressure 
distribution, it is possible to drive, in relative isolation in the source 
section, the plane and higher order angular modes, up to and including the 
(5,0) mode, without the presence of higher order radial mode response. 


VI. 1. Description of Experiments 


The experiments conducted with the five liners progressed from a rela- 
tively simple low frequency test to a reasonably complex, multi-mode, high 
frequency test. The degree of success achieved was almost directly propor- 
tional to the complexity of the test. The most simple and probably most 
rewarding test in terms of useful data was that run using a felted metal over 
a 1 i 1 1 deep honeycomb core. A series of tests was conducted at 1000 Hertz for 
both a plane (0,0) mode excitation and a (3,0) mode excitation. The tests 
were conducted with and without mean flow. The most complex and least produc- 
tive tests were those conducted with a felted metal liner over a £" deep 
honeycomb core at high frequencies. 


VI. 1.1. Liner description 

The five liners used In the test series and illustrated previously in 
figure 14 are briefly described below: 

(1) Liner set 1 was a single liner for the outer wall consisting of one 
layer of nominal 40 Rayl felted metal over a 1±" deep honeycomb 
with 1/8" cell diameter. The honeycomb was closed at the back with 
a .050" thick aluminum plate. The honeycomb was bonded to the 
aluminum plate and was tightly squeezed around the acoustic inner 
surface, thereby eliminating any change in porosity due to the 
adhesive. The Inner wall consisted of a i" thick steel pipe. 

(2) Liner set 2 was a 20 Rayl felted metal over a £" core. This con- 
figuration was used to make up both an Inner and outer liner for the 
annulus. The rigid closure for the honeycomb was a 0.050" thick 
layer of cured res in- impregnated fiberglass cloth. 

(3) Liner set 3 was the same as (2) above with a 40 Rayl felted metal 
(not the same manufacturer as for the acoustic material of (1) 
above) . 
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(5) Liner set 5 consisted of the outer wall covered with a 2" thick 
commercially available acoustic foam with closely controlled pore 
size, with the inner wall surface a i“ thick steel pipe. 


VI. 1.2. Data requi rements 

The type of data recorded varied from test to test, depending on the 
motivation for a particular test. In all cases, the circumferential sound 
pressure distribution was recorded for the outer wail position in the plane of 
the drivers. This was used to compare with the computed circumferential dis- 
tribution of sound pressure at the driver plane, for a particular mode simu- 
lation, to determine the validity of the test set-up and to determine if the 
presence of the liner had introduced any distortion in the driver plane sound 
pressure distribution. 

The sound pressure distribution in the plane at the intersection between 
the lined and unlined section was recorded in as much detail as possible. The 
SPL at this plane was chosen as the source distribution. Angular SPL distri- 
butions were recorded at nine equally spaced radial locations from the outer 
wall to inner wall. Also, continuous radial SPL distributions were recorded 
at several angular locations. Further, for several of the tests, these 
measurements were made at four axial positions. 

In tests where one mode was responding in isolation, it was possible to 
measure the axial decay and propagation speed of that mode by measuring the 
peak amplitude at the inner and outer walls for a number of axial stations. 

In addition, temperature, humidity, and barometric pressure were recorded 
for most of the tests conducted. No particular effect was related to changes 
in humidity or barometric pressure, but very significant changes in mode 
response amplitude, cut-on frequency and modal attenuation were observed for 
very slight changes in temperature within the duct. In fact, as noted 
previously, when conducting tests with flow, it was found to be necessary to 
introduce a temperature compensating circuit which maintains a constant value 
of wave number (k = 2Trf/c). 


VI. 1.3. Phenomena to be verified 

The most important phenomena to be considered in the lined duct tests 
are listed and briefly discussed below. 

(l) Natural acoustic pressure mode shapes. 
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(a) Angular modes. Provided the liner has no variation of proper- 
ties in the angular direction, and if no angular non-uniformity 
exists in the mean flow, then the angular mode shape should 
remain in a sinusoid variation. This, however, required experi- 
mental verification. 

(b) The radial mode shape can be critically dependent on both the 
value of wall impedance and the flow velocity profile. Radial 
mode predictions must be verified for as many combinations of 
wall impedance and flow as possible. 

(2) Complex propagation constant. In rigid wall ducts mode propagation 
is very distinctly defined. The propagation constant is either real 
(propagating) or imaginary (exponentially decaying). The presence of 
an attenuating liner results in a complex propagation constant, the 
real part normally corresponding to the rigid-wall propagation con- 
stant with the imaginary part representing the sound attenuation in 
the duct. The presence of mean flow further modifies the propagation 
constant. Both phenomena are subjects of study in the experiments 
for the 1 ined annular duct. 

(3) Radial acoustic velocity mode shape. Analytically, this was given as 
a solution to the radial momentum equation. The radial acoustic 
velocity is a function of both the acoustic pressure mode shape, the 
propagation constant and the mean velocity and mean velocity profile 
(or shear). The intent was to make radial measurements with a hot- 
wire probe and compare with the predictions. However, difficulties 
in processing the hot-wire signals precluded any useful observations. 

(4) Effect of termination impedance on wave field in duct. If the duct 
radiation impedance is not well matched, significant reflections 
occur at the termination. The effect of this on the energy attenua- 
tion rate and the standing wave in the duct was investigated. 

(5) Effect of mean flow on wave field in duct. Flow affects the wave 
field by convection through the duct and by refraction due to the 
flow profile. For certain conditions, either or both effects can 
be very significant. Also, they can either add or cancel according 
to the direction of propagation with respect to the flow direction. 

(6) Effect of locally reacting vs non-locally reacting liners. For the 
honeycomb type liners, there is no appreciable sound propagation 
through the liner, parallel to the duct axis, and use of the locally 
reacting impedance assumption appears to be valid. However, for 
certain of the foam and fiberglass liners, propagation through the 
liner material occurs, with significant alteration in the overall 
attenuation performance of the liner. Tests and comparisons with 
theoretical predictions illustrate this point. 
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VI. 2. Analytical Predictions 


The theoretical study of sound propagation in acoustically lined ducts is 
a very complex subject. An experimental program of the type described in this 
report would be impossible without guidance from theoretical predictions. 
Therefore, before a detailed description and presentation of the experimental 
results, a series of numerical studies which display the predicted behavior of 
the parameters of interest and which were used in designing the experimental 
study, is felt to be useful. 

The numerical studies are presented in graphical form to aid understand- 
ing of the effects. The graphs will show contours of equal attenuation and 
phase speed for individual modes and combinations of modes. Attenuation vs. 
frequency behavior for one of the liners is shown. Radial sound pressure 
distribution as a function of distance along the duct axis is studied, showing 
effects of variation and sensitivity to wall impedance on that radial pressure 
distribution. Also, the effect of wall impedance on the axial power level 
attenuation is studied. 

A computer program, based on the theory presented in the analytical 
review section of this report was prepared within the Lockheed-Georg i a Company 
under I R&D funding. Before entering a detailed discussion of the numerical 
studies, a brief description of the computer program is useful. 


VI. 2.1. Features of computer analysis for sound propagation in ducts 

Basically, the computer program provides mode shapes and eigenvalues 
of the convected wave equation (equation (14)) for sound propagation in ducts 
of rectangular, circular, or circular annular geometry with axi symmet r i c , 
radially non-uniform mean flow. A brief discussion of each feature of the 
program input and solution features follows. 

(1) Geometry. As stated above, the geometry may be either rectangular, 
circular, or circular-annular, 

(2) Wall impedance. The wall impedance may be arbitrarily specified. 
For a rectangular duct, different values of wall impedance can be 
specified for opposite walls. The other pair of walls is assumed 
to be rigid. For the annular duct, the inner and outer walls may 
have different impedance. 

(3) Circumferential or cross-mode order. Specification of the circum- 
ferential mode (“spinning mode") order for circular ducts or the 
cross-mode for rectangular ducts is arbitrary. The presence of 
rigid radial splitters can be included by specifying non-integer 
values for the circular geometry mode order. 

(4) Radial mode order. This input is arbitrary. Any number of modes 
can be calculated, up to the number of finite difference stations, 
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however, accuracy requirements would keep the number of modes at a 
factor of 5-10 below the number of computational stations. 

(5) Flow velocity and velocity profile. One of three options is 
ava i 1 ab le : 

(a) Experimental values may be input at the finite difference 
computational stations; 

(b) Mean velocity and boundary layer thickness can be specified, 
using a ha 1 f-s inewave to simulate the boundary layer profile; 

(c) A number of experimental data points can be input and the 
velocity and shear profiles close to the wall can be computed 
using a modified version of the universal law of the wall. 

(6) Frequency. This may be arbitrarily specified and can be input as 
real frequency or non-dimensional wave-number parameter. 

(7) Source sound pressure distribution. This has two options. One 
assumes a uniform sound pressure in the radial direction. Only 
radial modes of a particular angular mode order may be considered at 
any one time, A two-dimensional curve fit was not built into this 
program. The other option allows manual input of sound pressure 
distribution in the radial direction at each of the computational 
stations. A complex least-square curve fit routine is used in a 
Fourier type analysis to determine the mode amplitudes at the source 
location in order to simulate a specific source sound pressure 

d is trlbut ion . 

The computer program also includes effects of finite length, but in the 
present studies no end reflection is included. Thus, the results given are for 
an infinite length duct. 


VI, 2. 2. Numerical studies 


VI. 2. 2,1, Effects of wail impedance without flow . The parameter of primary 
usefulness in a design study is the total attenuation expected per unit length 
of duct. However, this cannot be specified unless the sound pressure distri- 
bution, and thus the amplitudes of the various modes, is specified at the duct 
entrance. Analytical studies in the past have resorted to various techniques 
for determining this total attenuation per unit duct length based on assumed 
modal distributions or initial pressure distribution. Feder and Dean (ref. 

48) experimented with three assumptions for the modal distribution. These were 
(1) to assume that all the energy propagates in the least attenuated mode, (2) 
to assume that all the propagating modes initially (at the liner entrance) have 
equal energy, and (3) to assume that all the propagating modes initially have 
equal pressure amplitude. Another assumption introduced by Rice (ref. 49 was a 
sound pressure distribution at the liner entrance which conformed to the hard 
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wall plane (0,0) mode or else was specified from measurement. For a single 
frequency of excitation, this assumption is more realistic, especially for 
comparison with experiments where the initial pressure distribution can be 
measured. 

In the discussion immediately following, it was assumed that the sound 
pressure distribution was confined to a particular angular mode set and that 
the sound pressure in the radial direction at the entrance to the lined 
section was plane. For purposes of comparing with experimental data from the 
facility, this condition is very closely approximated at frequencies below 
1620 Hertz. The numerical studies in this sub-section will concentrate on 
the attenuation in an annular duct with hub-tip ratio of 0.566 and an unlined 
inner cylinder. It is quite difficult to display the relationship between the 
complex impedance and the complex eigenvalue without resorting to a contour 
map. Two types of these are used in practice. Horse (ref. 1) used the eigen- 
value as the axis and ordinate of the complex plane plot and the complex im- 
pedance as the contour parameter. It seems more convenient, when studying a 
particular mode or a sum of modes to display the contours of constant eigen- 
value on an impedance plane plot. Thus, effects of variation of impedance can 
eas i ly be studied. 

In the presentation to follow in this section, several modes and combi- 
nations of modes will be displayed in an impedance plane plot with attenuation 
and phase speed the plotting parameters. These plots will be used to show the 
effect of, and sensitivity to, wall impedance. Because of the multiplicity of 
combinations which could occur if inner and outer wall impedance varied inde- 
pendently, the inner wall is either considered as rigid or else its impedance 
is set equal to that of the outer wall. 

Figure 87 illustrates the attenuation predicted for the 3rd angular mode 
set with the assumption of a plane radial sound pressure at the liner en- 
trance. This computation, which displays the power level attenuation at one 
duct width downstream is at a frequency of 1000 Hertz (kL = 1.915). Several 
interesting facts can be gleaned from this contour map. First, examination of 
the data which resulted in the attenuation contours revealed that the contri- 
bution of each mode, necessary to simulate the plane radial distribution at 
z = 0, varies for each value of impedance on the complex plane plot. This 
results from the fact that mode shapes are highly dependent on the wall impe- 
dance boundary conditions. In this particular study, only the outer wall 
contains an acoustic material, with the inner wall specified as rigid. Second, 
the result of this study shows that for the particular annular geometric con- 
figuration, the optimum impedance is approximately Z/pc = 0.05 - 0.10? with a 
power level attenuation greater than 30 ‘decibels at one duct width downstream. 

An interesting study, which illustrates the characteristics of the modes 
contributing to the contours of Figure 87 , is presented in Figures 88-91. 

Fi gure 88 is an impedance plane contour plot of the complex eigenvalue. The 
attenuation contours are Attn(dB) ~ 8.68 kLa and the contours normal to the 
attenuation are the real part of the propagation constant, 3. The details of 
this curve provide data which can be directly compared with measurements from 
the facility. At frequencies below 1620 Hertz, the hard wall (3,0) mode can 
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Frequency - 1000 Hz; kL = 1.915; Hub-Tip Ratio = O .568 
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Figure 88- Isocontours of (3,0) Modal Attenuation and Propagation Constant, 3, 
in Annular Duct 
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89- Variation of (3,0) Modal Amplitude Along Trace of 8 dB 
Contour in Figure 88 


Figure 90 . Variation of (3,0) Modal Amplitude Along Trace of 16 dB Contour 
in Figure 88 





be simulated in the V long source section. This hard wall (3,0) mode, after 
entering the test section, which is nine duct widths long, rapidly converts to 
the (3,0) lined duct mode for most values of impedance on the contour chart. 
Thus, except for impedance values below the branch line, down to reactances of 
approximately ”0.6, the (3»0) lined duct mode can be set up in relative isola 
tion at about two duct lengths from the entrance, since all the higher order 
softwall modes decay much more rapidly than the zero order mode. In the region 
excluded, other modes may have less attenuation than the (3,0) mode. In all 
but the excluded region, it should be possible to measure both the attenuation 
rate ■snd the propagation speed of the (3,0) softwall mode. 

The phase velocity of propagating modes is c^ = c/3, where c is speed of 
sound and 3 is the propagation constant. It is quite intriguing in examining 
figure 88, to note that at resistances less than 0.73, the propagation con- 
stant exceeds 1.0 for most reactance values below the branch line. The 
physical interpretation of a mode with a propagation constant greater than 1.0 
is that the mode is “slow. 11 That is, the phase speed, c^, is less than the 
speed of sound. It is not known if these modes have been experimentally 
detected previously, but the possibility of their existance was noted in a 
conversation with Mungur in 1969. One physically significant point associated 
with the slow propagation speed is that the axial wavelength is reduced at a 
particular impressed frequency (i.e., A mn = c/f3 ). 

It was mentioned earlier that while the contours of figure 88 display 
the attenuation and axial propagation characteristics for the (3,0) mode, the 
radial shape of the (3,0) mode changes drastically for various positions in 
the impedance plane. Figures 89 and 90 illustrate the changing character of 
the mode shape as two of the constant level contours are traversed in a clock- 
wise direction. Mode shape amplitudes are shown in figure 89 for various 
positions in the impedance plane, corresponding to intersections of the 8 dB 
contour and reactance lines beginning at 1.1 and progressing in a clockwise 
direction in steps of 0.1 magnitude in reactance. The mode shapes are normal- 
ized to 120 decibels at the outer wall. For large positive values of reactance, 
sound pressure is greater on the inner hard wall, than on the outer soft wall. 
However, for the low positive and all negative values of reactance, the pres- 
sure is lower on the hard wall, than on the lined wall. At first glance, this 
is adverse to physical concept since one would expect the lined wall to have a 
lower sound pressure for all modes. However, this adverse skewing of the mode 
shape is due to the effect of the reactance. It is observed that the influence 
of the hard wall on forcing a zero gradient at the wall extends quite a long 
way toward the center of the channel. In fact, in some impedance regions 
finite sound pressure exists only near the wall. The mode is very similar to 
a surface wave which decays exponentially away from the wall. Consistent with 
this, the value of sound pressure on the hard wall in the region of small 
resistance and negative reactance appears to approach zero in the limit of 
zero resistance. It is noted, however, that this would violate the boundary 
condi t i ons . 

The set of mode shapes in a clockwise direction around the 16 dB contour 
are very similar, except for a slight dip in some of the shapes. This dip 
sometimes makes mode identification very difficult, especially in the region 
of the b ranch line. 



A characteristic of this first order mode, which is very easy to measure, 
is the SPL and phase difference between the inner and outer walls. These 
parameters are also fairly easy to plot as contours, somewhat similar to the 
plots of eigenvalue. In figure 91 » the differences between outer wall and 
inner wall SPL and phase are plotted as contours of constant level. At first 
glance, the usefulness of a plot of this type appears to be minimal, but as 
will be described later, it provides a check on the effective value of wall 
impedance also obtained from other experimental data. 

A considerable volume of published results deal specifically with the 
least attenuated mode. Also, as mentioned, much of the current and previous 
liner design work is based on the least attenuated mode results. Thus, it is 
of interest to study this mode both analytically and experimentally and also 
to compare the (0,0) mode eigenvalue with that of the (3,0) mode previously 
discussed. First, the (0,0) mode is not always the least attenuated mode. 

For certain reactance values immediately below the branch line, the (0,1), or 
in some cases the (0,2) mode will be least attenuated. 

In the comparison of the (0,0) (figure 92) and (3,0) mode (figure 88) 
contours of equal attenuation, it is immediately obvious that at low resis- 
tance, for positive values of reactance, the contours are normal to the branch 
line for the (0 ,0 ) mode and normal to the zero resistance axis for the (3,0) 
mode. The most useful observation is that the (0,0) mode attenuation is 
generally less than that for the (3,0) mode for a particular value of im- 
pedance. Probably the most significant difference between the impedance plane 
eigenvalue contours is the location of the branch point. The lower order mode 
branch point lies closer to the zero reactance axis. Also, for any value of 
impedance, the phase speed is much slower. 

The propagation of a plane wave of zero angular order is shown in figure 
93. If this is compared with figure 87, two distinct differences can be 
observed. First, the attenuation for most values of impedance is significantly 
lower. The maximum attenuation for the plane mode is of the order of 15 dB at 
the optimum, while the 3rd order, plane radial shape has a maximum attenuation 
of over 30 dB. Second, the optimum value of impedance appears to be approxi- 
mately, Z = 0.15 "0.3i, compared to 0.05 " 0. 1 i for the 3rd order plane mode 
of figure 87. 

Two sets of the test liners have both inner and outer attenuating walls. 
Figures 94 and 95 show the contours of constant attenuation for the (2,0) and 
(2,1) modes at 1800 Hz over a range of impedance encompassing the test values 
for those liners. For this particular impedance-f requency range, both modes 
are propagating with predicted attenuations from 0.5 to 9 decibels. 


VI. 2. 2. 2. Effects of frequency without flow . Another parametric study which 
is of interest^ is' given in figure 96 and" shows the functional behavior of 
the 20 Rayl configuration with frequency with both inner and outer walls 
lined. This chart illustrates the attenuation-frequency cha racter i s t i cs of 
the first three radial modes, for a sin(2G) angular mode variation, up to a 
frequency of 4000 Hertz. Average values of measured impedance were used in 
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Figure 94 . Isocontours of (2,0) Modal Attenuat ion i n Annular Duct 
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Figure 95 . Isocontours of (2,1) Modal Attenuation in Annular Duct 
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computing this frequency response. Similar modal attenuation curves are 
shown in figures 97-99 for the plane (zeroth), fifth and sixth angular 
modes. 

The most interesting result from this study concerns the least attenuated 
mode. If the criterion chosen was that the least attenuated mode dominates 
the liner attenuation, then the maximum attenuation at any frequency for this 
liner would be 6 dB per duct width, with the peak at 2900 Hz, where the wall 
impedance is approximately 0.A - l,2i. However, if a radial source pressure 
distribution is assumed, which is independent of frequency , a completely 
different trend results at frequencies above the first radial mode cut-on 
frequency. This demonstrates the necessity, once again, for more detailed 
knowledge of the in-duct sound pressure distribution due to the noise source. 

In this numerical study a plane radial shape was introduced for all 
frequencies. Coefficients of each mode were determined to give the best 
approximation to the plane wave at the liner entrance. This procedure was 
repeated at each frequency and the resulting sound power attenuation was 
computed at three axial stations (i.e. 1, 2, and 3 duct widths downstream). 
Similar results are shown for the other angular mode sets in figures 96-99 
with the curves being identified in the legend. It is very interesting to 
note how the power attenuation at one duct width adhers to the attenuation of 
the least attenuated mode up to approximately 3000 Hertz, but then becomes 
considerably greater than the attenuation of that mode (which incidently has 
switched from the zero order radial mode to the 2nd order radial mode). It is 
in this region (i.e. the region near liner resonance) that the picture of 
modal attenuation becomes very confusing and it is considered an absolute 
requirement that an initial sound pressure spatial distribution must be given. 

Another interesting feature of the power attenuation curves for an initial 
plane wave is the tendency of the frequency of peak attenuation to shift as a 
function of axial distance. This is a result of the modes with high attenua- 
tion losing effectiveness at distances further downstream from the source 
p 1 ane . 


Vl - 2 * 2 - 3 : Effects of flow. In previous numerical studies, flow effects on 

attenuation rate and propagation speed have been found to be very signifi- 
cant. ^ In addition, mean flow indirectly affects liner performance by modi- 
fication of liner impedance. These effects on liner impedance were discussed 
in Section V. 


A study which duplicates the data of figure 88, but with 50 fps mean flow 
velocity through the annulus, is shown in figure 100. This chart is a con- 
stant level contour plot in the impedance plane for the (3,0) mode at 1000 
Hertz. The mean flow profile is that of a quarter-sine wave in the 2 % thick 
boundary layer { 2 % of annulus width) and uniform flow elsewhere. General 
observations are that for negative reactance values the attenuation is reduced 
by approximately 1 decibel and that the branch line has shifted to more posi- 
tive values of reactance. 
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Figures 101-103 show the attenuation contours for an initial radial plane 
distribution (derived by Fourier summing 10 radial modes at the source plane) 
for the third angular mode at 50 fps , with 7 % boundary layer, for three dif- 
ferent downstream positions. Figure 101 at one duct width downstream is 
comparable to figure 87 for no-flow. An overlay of the two. plots would show 
that the position of the optimum impedance with flow has shifted upward and 
toward greater positive values of wall resistance. At the same time the 
maximum attenuation appears to have been reduced. Figures 102 and 103 ate 
similar to figure 101 at two and three duct widths downstream. The contours 
in figure 102 are plotted in multiples of two times the values of attenuation 
used for contour plotting increments in figure 101. A multiple of three was 
used in plotting the contours at three duct widths downstream. Thus, when 
comparing the attenuation contours at the three positions, the deviation from 
a linear attenuation rate can be ascertained. It is apparent that at distances 
further downstream, the optimum value of impedance for maximum attenuation 
shifts to a higher value of resistance. This effect is a result of the fact 
that modes with greatest attenuation are, for all p ract i cal pu rposes , com- 
pletely decayed at 3 duct widths downstream. Thus, the optimum shifts from, 
the lower resistance, associated with higher order modes, to the higher resis- 
tance optimum value, characteristic of the least attenuated mode. 

Finally, in order to provide further comparison data for the test results, 
contours for the (3,0) mode complex propagation constant are p lotted in figure 
1 04 and the 3rd order plane mode are shown in figure 105 for a flow velocity of 
150 fps. The predominant effect of flow on the (3,0) mode complex eigenvalue 
is a reduction in propagation constant, 3, as would be expected for the down- 
stream propagating mode and as can be seen by comparing figures 88 and.lOA, 

The data necessary to fill in the area in. the vicinity of the branch 1 ine.was 
not computed because of the extra computer run time required for that region. 

In examining the plane mode attenuation for the 3rd angular set, the 150 fps 
data of figure 105 should be compared with figures 87 and 101. Although the 
area of the optimum impedance was not bracketed, it appears that the value of 
optimum impedance has increased in resistance and the reactance has assumed a 
more positive value, when compared to the no- flow case. 


VI. 2. 2. 4. Sound pressure distribution across the duct . An example of the 
effect of wall impedance on the pressure distribution in the duct is shown 
in figures 106-109. These curves at the four axial stations, for each 
impedance value, originate as plane waves at z/L = 0 and ultimately revert 

to the shape of the least attenuated mode. In figure 106, the shape 
remains relatively plane at all axial positions, however, the attenuation 
per duct width is not linear. A glance at figure 110 reveals the behavior 
of this non-linear attenuation vs axial distance, and a study of the propa- 
gation constants reveals that mode "beating" is the reason for the non- 
linear behavior. The case studied in figure 106 has two modes with almost 
equal attenuation but with significantly different propagation speeds as shown 
in the tabulation at the top of figure 110. This particular situation results 
in the two modes slowly moving in and out of phase as the pattern propagates 
down the duct. Figures 107-109 show variations of the downstream pressure 
distribution for different values of impedance in the region near the optimum 
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outer wall reactance 


Frequency - 1000 Hz; kL « 1.915; Hub-Tip Ratio = 0.568 
Hard Inner Wall; Acoustically Lined Outer Wall; 

Velocity Profile i Sine Wave; Boundary Layer Thickness 24 


Frequency - 1000 Hz; kL *= 1.915; Hub-Tip Ratio = 0.568 
Hard Inner Wall; Acoustically Lined Outer Wall; 

Velocity Profile i Sine Wave; Boundary Layer Thickness 24 




Figure 100. Isocontours of (3,0) Hodal Attenuation In Annular Duct with 50 fps 
Sheared Flow 
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Figure 101. Isocontours of Sound Power Attenuation One Diameter Downstream In 
Annular Duct for Plane Radial - Third Angular Mode Set with 50 fps 
Sheared Flow 


Frequency - 1000 Hz; kL “ 1.915; Hub-Tip Ratio = 0.568 
Hard Inner Wall; Acoustically Lined Outer Wall; 

Velocity Profile i Sine Wave; Boundary Layer Thickness 24 



Figure 10Z. (socontours of Sound. Power Attenuation Two Diameters Downstream in 
Annular Duct for Plane Radial - Third Angular Mode Set with 50 fps 
Sheared Flow 


Frequency - 1000 Hz; kL = 1-915; Hub-Tip Ratio - 0.568 
Hard Inner Wall; Acoustically Lined Outer Wall; 

Velocity Profile i Sine Wave; Boundary Layer Thickness 24 



Figure 103 . I socontours of Sound Power Attenuation Three Diameters Downstream In 
Annular Duct for Plane Radial - Third Angular Mode Set with 50 fps 
Sheared Flow 
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Frequency - 1000 Hz; kL — 1.915; Hub-Tip Ratio — 0.568 
Hard Inner Wall; Acoustically Lined Outer Wall 
Velocity Profile £ Sine Wave; Boundary Layer Thickness 2 % 



Figure 10** . Isocontours of (3,0) Modal Attenuation and Propagation Constant, B, 
in Annular Duct with 150 fps Sheared Flow 


Frequency - 1000 Hz; kL = 1.915; Hub-Tip Ratio * 0.568 
Hard Inner Wall; Acoustically Lined Outer Wall; 

Velocity Profile £ Sine Wave; Boundary Layer Thickness 2 % 



Figure 105 • I socontours Sound Power Attenuation One Diameter Downstream in 
Annular Duct for Plane Radial - Th i rd Angular Mode Set with 150 fps 
Sheared Flow 
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impedance. Note in figure 110 how the performance ranking of various liners 
at one duct width downstream changes for distances further downstream. 

There are many interesting variations on this type of presentation; 
however, this serves to illustrate the nature of energy attenuation of a 
particular pressure distribution and how that parameter compares with modal 
attenuation. As an example of a ridiculous comparison, observe the modal 
attenuation for a wall impedance Z = 0.6 - 0.2i. The attenuation of the least 
attenuated mode i s 19.*+ dB and the next mode has 19.8 dB attenuation, with the 
other modes having 50 dB or greater attenuation. At first glance, it would be 
expected that the minimum attenuation for any type of sound pressure wave 
propagating down the duct would be at least 19.4dB. However, the calculated 
energy attenuation which is proportional to 

R q 

1/2 f pv*RdR + Mach no terms, (8°) 

i 

is only 12.5 decibels at one duct length and at three duct lengths becomes 45 
decibels, still 13 decibels below the attenuation of the least attenuated mode. 
The cause for this significant disparity is due almost entirely to mode phasing, 

VI. 3 Comparison of Theory with Experiment 


The previous sub-section on analytical predictions has served to illus- 
trate the type of results expected from the experiments. Also, in the intro- 
duction to this section, the major phenomena for experimental consideration 
were listed. These were (l) propagation constant (attenuation and phase speed), 
(2) sound pressure mode shapes, (3) radial acoustic velocity mode shapes, (4) 
effect of end impedance, (5) effect of mean flow, (6) effect of wall impedance, 
and (7) effect of locally reacting vs non-local ly reacting wall impedance. 

In this section, a detailed set of measurements are presented, followed 
by the relevant comparisons of those measurements with theoretical predictions. 
The fundamental parameters derived from basic theory are the modal eigenvalue 
and eigenvector. Thus, the complete section will be devoted to various com- 
parisons with either the eigenvalue or eigenvector, singularly or in combina- 
tion, together with studies of the effects of flow and wall impedance on those 
pa rameters . 

A considerable portion of the work which follows is devoted to comparisons 
in lined ducts without flow. These tests were multi-purpose. First, it was 
necessary to develop the necessary test techniques and data analysis methods 
prior to introducing the additional complication of mean flow, and second, it 
was necessary to verify the theory for the no-flow case prior to testing with 
fl ow. 

Thus, the following major sub-sections are devoted to no-flow eigenvalue 
comparisons of theory and experiment, no-flow eigenvector comparisons of theory 
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and experiment, and the effect of flow on both the eigenvalue and eigen- 
vector. 


VI. 3 . 1 . Modal attenuation and propagation constant ~ no flow 

In the introduction to this section, it was impl ied that at frequen- 
cies below the first radial mode cut-on it would be relatively easy to measure 
the mode attenuation rates and propagation constants, but that at higher 
frequencies, where more than one radial mode is propagating, determination of 
the characteristics of individual modes would require a more sophisticated 
experimental technique. Because of the different requirements of the experi- 
mental procedures, results of tests from two different frequency regimes will 
be described separately. 

VI. 3. 1.1. Single lower-order radial mode . Liner 1 was designed to have maxi- 
mum attenuation in the frequency range below the first higher order, hard-wall, 
radial mode cut-on frequency. Thus, with this liner, a series of angular 
modes could be studied in isolation in both the hard wall source section and 
in the lined test section. 

In comparisons of experiment and theory, the numerical predictions 
utilized the measured values for impedance as discussed and presented in 
Section V. 

Liner 1 was very useful in two ways. Because of its i mpedance- f req uency 
characteristics, much useful data for comparison of theory and experiment for 
the lowest order radial mode was developed from this liner. It was also used 
for developing data analysis techniques for measurement of higher order radial 
mode characteristics. This will be discussed in the sub-section immediately 
fol lowi ng. 

A useful series of basic tests was conducted with the liner 1 configura- 
tion. A 3rd order angular mode (sin 39) was set up in the source section. 
Radial sound pressure surveys were conducted at an angle corresponding to a 
peak value of the angular mode at a large number of axial stations. An 
interesting parameter is the axial variation of the sound pressure amplitude 
and phase. Figure 111 shows SPL and phase at the inner and outer wall for 
three frequencies plotted versus axial distance down the duct. Examining the 
sound pressure amplitude decay characteristic for 88 A Hertz, it is observed 
that the SPL at the inner wall falls more rapidly than does that at the outer 
wall. The radial pressure pattern at the liner entrance has a 2.5 dB differ- 
ence between the inner and outer wall, which is approximately the level 
difference expected for the hard wall ( 3 , 0 ) mode excited at the source plane. 

As this wave enters the lined section, the natural modes abruptly assume a 
new radial shape. Thus, some of the energy transfers into the least attenu- 
ated radial mode and some transfers into higher order radial modes. In the 
first 4 " - 6 “ , the effect of these higher order modes can be detected by 
observing that the attenuation rate on the outer wall is non-linear, 
signifying that the modes are decaying at different rates. In figure 111(c) 
for 1 4 30 Hertz, it is found that the outer wall now has a lower sound level 
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than the inner wall. This is primarily attributed to the change in the 
shape of the least attenuated mode, which will be shown later. It is empha- 
sized that this data is acquired in the first 12 M of the duct liner in order 
to avoid any waves which may be reflected from the open end. As the attenu 
ation at these test frequencies is of the order of 1.5 dB per inch, the 
attenuation from station 24 to the open end would be 36 dB, with a further 
36 dB for the reflected wave traveling back to station 2k . Thus, at 24" any 
reflected wave from the open end would be at least 70 dB below the incident 
wave . 

If one attempts to determine the linear attenuation rate for these fre- 
quencies, it is wise to use the slope between 32" and 2k ", since the higher 
order modes have, for all practical purposes, decayed and little chance for a 
reflected wave exists. The attenuation rates of the inner wall, the outer 
wall, and an average between these two is given in table V. * 


TABLE V 

COMPLEX PROPAGATION CONSTANT FOR kO RAYL, 1 .5” DEEP LINER FOR (3,0) MODE 


Frequency 
He rtz 

INNER WALL 

OUTER WALL 

AVERAGE 

At tenuat ion 
dB/Duct 
Wi d th 

Wave 

Length 

Constant 

At tenuat i on 
dB/Duct 
Wi dth 

Wave 
Leng th 
Constant 

Attenua ti on 
dB/Duct 
Wi dth 

Wave 

Length 

Constant 

884 

6.5 

.665 

5.6 

.644 

6.0 

.654 

1000 

6.8 

.705 

5.5 

.733 

6.2 

.719 

1430 

§.3 

.710 

6.3 

U~l 

00 

VJ"1 

6.3 

.648 


The computed attenuations and wavelength constants are compared with the 
average measured values in table VI, using the impedance values derived from 
impedance tube measurements as the boundary condition for the computer program 
which solves the convected wave equation in the duct. 

Two values are given for the computed attenuation at each frequency. The 
first is the modal attenuation and the second is the power levc 1 a ttenuat i on 
associated with a wave which was radially plane at the test section entrance, 
but which decayed to the shape of the (3 ,0) mode , as the higher order lined duct 
modes were attenuated. The computed wavelength constants are for the same 
conditions. The comparison between measured and calculated values of attenu- 
ation is excellent. It would be difficult to achieve a more accurate measure 
ment or comparison. The wavelength (or propagation) constant, 3, comparison 
is not nearly as good, but on the whole this comparison is also quite satis- 
factory. 
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TABLE VI 


COMPARISON OF CALCULATED AND MEASURED COMPLEX PROPAGATION CONSTANT 

(3,0) MODE 


F requen cy 
He rtz 

ATTENUATION 

WAVELENGTH CONSTANT, 3 

COMPUTED 

Avg. 
Meas . 

COMPUTED 


Moda 1 

Power Level 
For Plane Wave 

Modal 

P 1 ane 
Wave 

884 

6.2 

6 . 1 

6.0 

.556 

.571 

.654 

1000 

6 . 0 

6.0 

6.2 

.679 

.698 

.719 

1430 

5.8 

6.0 

6.3 

.682 

.698 

.648 


The usual indicator of a liner's performance is its frequency response 
(or attenuation) curve. Figure 112 is a plot of the power level attenuation 
at three axial stations for the frequency range 400-4000 Hertz for an assumed 
radial plane wave at the duct entrance. The attenuation, as stated in a 
previous discussion, is not a linear function of axial distance because of 
the different exponential decay rates and phase speeds of the various modes 
contributing to a particular pressure distribution. Thus, the average atten- 
uation at each of these three axial stations is determined from figures 
111(a) - (c) at stations 32, 28, and 24. These data are plotted in figure 
112 in a direct comparison with the computed attenuation, and reflects the 
excellent comparison made in table VI above, using impedance tube data for the 
calculation. However, if anything, the comparison in figure 112 is better, in 
that it shows the non-linear effect due to radial mode combinations, necessary 
to match the test section entrance pressure distribution. It is unfortunate 
that a test was not run at 1200 Hz to show the weak peak in the liner attenua- 
tion and again at a higher frequency to show the reduced efficiency of the 
liner; however, the three data points shown are sufficient to give a tentative 
verification of the theoretical prediction method. 

However, one matter of concern which must be discussed before accepting 
the above results at face value is that of wall impedance spatial variation. 

A study of the impedance measurements at 1000 Hertz reveals that, by the two 
microphone method^the mean acoustic impedance is Z = 1.25 - 0.93 i with a 
standard deviation of 0.1 on resistance and 0.05 on reactance. If this value 
of impedance is entered on the curve of figure 87, the attenuation and propa- 
gation constants derived are 5.0 to 5.8 dB for a plane wave at one duct width 
downstream and from figure 88, a modal attenuation of 5.8 to 6.8 dB with 
3 - .76 to .80. The plane wave attenuation compares reasonably with the 
measured value of 5.1 dB from figure 111(b) and the modal attenuation compares 
quite well with the measured value of 6.2 dB, also derived from 111(b) as 
discussed previously. However, the calculated propagation constant derived 
from the impedance measured by the two microphone method is significantly 
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Figure 1 1 2 . Compar 1 son of Experiment and Theory for Power Level Attenuation in 
Liner 1 
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different from that calculated using the impedance tube data, and it is 
obvious that the measured propagation constant lies somewhere in between the 
two computed values. If the measured attenuation of 6.2 dB and the measured 
propagation constant of 0.72 is also entered on figure 88, it is seen that 
the measured value of complex propagation constant lies on a direct line be- 
tween the two computed values. At this point a very critical observation (or 
perhaps hypothesis) can be made. If the measured complex propagation constant 
is taken as accurate, then by an inverse process an impedance value can be 
inferred. This impedance value would not only assume that the duct propaga- 
tion theory is correct, but also that the measurement of complex propagation 
constant is accurate. The value inferred for the effective wall impedance for 
the liner is Z = 1,47 - 0.87i. 

Thus, three values of impedance are inherent in the measurements and it 
becomes quite obvious that an accurate measure of the impedance is crucial. 

The impedance tube data could be in error, since only one sample was tested, 
and quite obviously there is no way of comparing the impedance tube data with 
a mean impedance, since the impedance of the one sample in relationship to the 
mean impedance of the liner is unknown. 

In measurements by Helling (ref. 51), it was noted that, on average, the 
two microphone method results in a lower resistance and a higher mass reac- 
tance measurement than that derived from the impedance tube. Thus, it is 
quite possible that the value of impedance indicated for the two microphone 
(or one microphone traverse) method is in error by a small percentage. 

Melling was not able to conclude if the difference was due entirely to errors 
in measurement (primarily in phase) or i f a real difference exists between the 
two techniques. Dean (ref. 52) ran a controlled test in an impedance tube on 
a 1 1 % perforate in which the standing wave method, the two microphone tech- 
nique and the one microphone traverse technique were compared. Figure 113 
shows the results. The two microphone and impedance tube produce very nearly 
the same results up to 3000 Hertz, but the one microphone traverse method 
shows a very distinct difference. Dean attributed this difference to errors 
in estimating the sound pressure level and phase at the front and back surface 
of the cavity, as discussed in more detail in another section of this report. 


VI. 3. 1.2. Radial characteristics of the 'plane 1 mode . Following the above 
detailed discussion of characteristics of higher order angular mode complex 
propagation constant, it is also interesting to briefly examine the propaga- 
tion of the 'plane' (0,0) mode. The measured amplitude and phase change down 
the duct for an (0,0) mode simulation is shown in figure 114, Although the 
inner and outer wall sound pressure levels should be identical in the hard 
wall source section (Station 37), it is noted that the inner wall is 1.5 
decibels higher than at the outer wall. This measurement is indicative of 
the fact that more than one angular mode exists in the hard wall section. The 
characteristic radial shape of the hard wall angular modes has a lower level 
at the inner than at the outer wall, which is in opposition to the observation 
of a higher inner wall level. However, wi^h more than one angular mode 
present, the zero order radial mode of the various angular modes can add up to 
give the observed behavior. Thus, it must be concluded that a reasonable 
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simulation of the (0,0) mode did not exist for this test. Figure 115 is a 
sample of the angular distribution of sound pressure and phase in the hard 
wall section of the duct. A complex Fourier analysis of this spatial distri- 
bution reveals the angular mode distributions listed in table VII. 


TABLE VII 

ANGULAR MODE DISTRIBUTION FOR 
PLANE (0,0) MODE SIMULATION 


Angu 1 a r 
Mode No. 

Amp 1 i tude 
Decibels 

Phase 
Deg rees 

Phase 

Shift 

Degrees 

0 

118.6 

Hi .9 


1 

103.4 

190.4 

92.7 

2 

86.8 

150.8 

129.4 

3 

100. 1 

83.2 

134.2 


Thus, the (1,0) angular mode is approximately 15 decibels below the (0,0) mode 
or a linear amplitude of 0.17 when normalized to a plane mode of unity value; 
and the 3*0 mode, at 18.5 dB down, is 0,12 relative to the plane mode ampli- 
tude. Since these modes propagate at different speeds and have different 
attenuation rates (i.e., for the (3,0) mode 3 = .68 and the attenuation is 6 
dB per duct width, compared to 3 * 1.08 and an attenuation of 3.5 to 4.0 for 
the plane (0,0) mode), it would be expected that the amplitude difference 
between the inner and outer wall would change down the duct until the effect 
of the higher order modes ceased to be important . However, at 3 duct lengths 
the (1,0) mode would have attenuated, at most, only 2 decibels more than the 
(0,0) mode, thus even at station 24, the effect of the higher order angular 
modes is still present and the composite decay rate should be somewhat higher 
than that expected for the (0,0) mode. 

In spite of this higher order angular mode “noise" the (0,0) mode con- 
trols the propagation constant and compares quite well with the calculated 
value (viz. 3 - 1.08 theoretical vs. 1.08 experimental), whereas the measured 
attenuation rate seems to be about 1 dB higher than the theoretical value (3.8 
dB theoretical vs. 4.8 dB experimental). Theoretical values are derived from 
effective wall impedance and experimental values are averages of inner and 
outer wall values. It is interesting to note that in this particular case, 
the mode is a “slow mode." the propagation constant being greater than 1.0 (or 
conversely, the phase speed, c/6, is less than the speed of sound, c, from 
which the 'slow mode' terminology results). 


VI. 3*1. 3. Fourie r-Besse 1 analysis technique for radial modes . For the case 
of multiple mode propagation, the computerized two-dimensional complex 
1 Fou rier-Bessel ' analysis developed to determine the complete modal content at 
each axial plane, was used. As a means of determining the validity of this 
analysis the angular spatial patterns associated with the (3,0) mode at 1000 
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TEST SECTION ANGLE IN DEGREES 


Figure 115. Example of (0,0) Mode Angular Distribution in Hard Wall Sect ion Upst ream of Liner 


PHASE ANGLE IN DEGREES 



ANGULAR MODE COMPLEX AMPLITUDE AND PHASE 



O 

1 

2 


3 



4 


AXIAL 

STATION 

RADIAL 

STATION 

AMP 

dB 

PHASE 

DEG 

AMP 

dB 

PHASE 

DEG 

PHASE 

SHIFT 

DEG 

AMP 

dB 

PHASE 

DEG 

PHASE 

SHIFT 

DEG 

AMP 

dB 

PHASE 

DEG 

PHASE 

SHIFT 

DEG 

AMP 

dB 

PHASE 

DEG 

PHASE 

SHIFT 

DEG 


1 

09.7 

121.1 

93-4 

287.7 

158.0 

98.7 

32.4 

31.2 

120.3 

329.5 

104.8 

79.1 

76.9 

95-7 


2 

66.9 

114.9 

92.0 

271.0 

162.2 

96.I 

66.6 

43.O 

120.3 

325.6 

104.0 

87.1 

119.7 

157.5 


3 

89.2 

116.6 

93-1 

277.6 

161.1 

98.6 

69.5 

34.4 

120.3 

319.9 

103.2 

82.5 

66.7 

107.2 


4 

67.9 

113.4 

91.3 

262.1 

160.9 

96.2 

53.7 

39.6 

120.1 

315.6 

101.2 

65.5 

110.6 

150.7 

36 

5 

89.7 

109.4 

92.9 

266.0 

156.7 

96.I 

54.2 

33.5 

119.6 

311.6 

100.1 

83.I 

91.1 

112.3 


6 

68.8 

115.8 

93-7 

260.1 

157.2 

98.5 

52.9- 

32.1 

119.3 

309.9 

99.2 

79-3 

322.8 

7-3 


7 

91.9 

78.1 

93.1 

254.5 

158.8 

97.8 

54.5 

39-5 

118.8 

311.0 

103.4 

84.1 

128.8 

179.1 


8 

92.8 

83.6 

94.9 

262.3 

158,8 

98.3 

55.4 

39-4 

118.8 

306.8 

105.3 

79-9 

53-9 

112.5 


9 

91.8 

83.7 

94.5 

264.8 

156.8 

97.6 

53-7 

41.0 

118.8 

306.5 

109.1 

79.7 

342.9 

66.2 


mm 

87.3 

272.1 

87.1 

39-6 

154.1 

94.0 

169.5 

44.4 

114.7 

45.1 

109.1 

80.7 

198.3 

173.6 


I 

87.1 

265.5 

85.6 

23.6 

155.7 

93 .9 

161.5 

44.6 

114.3 

39.9 

109.1 

80.9 

199.7 

172.3 


I 

86.9 

256.1 

86.2 

19.6 

152.2 

93-2 

158.0 

43.1 

113.7 

32.0 

108.4 

78.3 

186.8 

179-8 


1 

85.1 

248.7 

85.5 

15.1 

154.0 

93.4 

150.2 

39-4 

113.3 

26.3 

105.4 

77.6 

179-9 

179.4 

32 

5 

84.7 

234.9 

85.6 

7.7 

153-5 

92.8 

146.5 

38.I 

112.8 

19.5 

105.2 

75.0 

161.0 

I65.O 


6 

84.7 

220.5 

86.3 

355-6 

154.5 

92.5 

132.5 

48.1 

112.1 

12.1 

107.9 

82.6 

2.7 

30.6 


■9 

83.8 

211.7 

85.9 

342.8 

158.8 

92.2 

126.2 

50.1 

111.6 

6.9 

108.1 

83.I 

358.0 

31.8 


H 

83.2 

202.9 

85.8 

336.6 

160.5 

92.2 

123.2 

51.4 

111.2 

3-2 

107.4 

83.5 

356.5 

35-9 


28 


84.7 

202.5 

86.6 

347.2 

83.3 

352.0 

84.9 

128.1 

82.1 

326.4 

85.6 

101.4 

82.6 

314.6 

86.1 

100.0 

82.7 

308.3 

86.2 

99-5 

81.8 

306.2 

86.2 

95.3 

82.0 

305.4 

86.0 

96.5 

80.6 

306.8 

83.4 

95-8 

80.2 

302.8 

82.8 

91.5 

80.1 

302.6 

82.8 

88.8 

79.5 

118.9 

81.8 

264.2 

78.3 

109.4 

80.6 

246.1 

73.7 

82.9 

79.1 

232.4 

73.7 

81.2 

78.0 

231.3 

74.1 

69.9 

78.7 

220.3 

73-4 

69.4 

77.0 

215.5 

73.1 

67.7 

76.4 

211.1 

73-1 

68.2 

76.5 

210.6 

73.7 

67.3 

76.2 

211.1 


166.0 

92.2 

126.3 

53-3 

111.1 

2.0 

109.0 

160.1 

90.4 

256.7 

'48.2 

109.3 

117.8 

109-2 

167.0 

89.6 

246.8 

45.2 

108.0 

103.4 

115.2 

162.9 

89.8 

234.6 

46.0 

108.3 

93.9 

114.7 

116.3 

90.1 

239-6 

46.7 

108.2 

94.6 

114.7 

166.4 

89.2 

235.5 

45.8 

107.3 

88.6 

115.0 

163.7 

88.1 

220.7 

64.4 

106.8 

84.5 

116.8 

159.5 

87.8 

217.3 

69.3 

106.3 

81.3 

116.6 

158.4 

87.6 

213.4 

71.2 

105.8 

78.2 

116.5 

161.4 

88.0 

215.9 

70.6 

105.8 

78.7 

117.7 

160.7 

85.I 

25.7 

57.0 

102.0 

232.5 

119.7 

163.6 

84.0 

15.0 

54.2 

100.7 

224.2 

119.5 

166.3 

83.8 

2.7 

55.6 

99-5 

211.1 

120.5 

I63.8 

82.7 

353.2 

58.7 

98.2 

201.8 

120.7 

168.5 

82.3 

351.2 

57.4 

97.3 

194.7 

123.5 

167.2 

81.7 

341.7 

59.7 

96.4 

187.1 

122.6 

I69.I 

81.1 

339.6 

59.7 

95.7 

182.6 

122.5 

166.9 

80.6 

337.8 

61.8 

95.1 

179.1 

123.5 

164.4 

80.3 

337.2 

62.9 

94.9 

178.3 

123.0 


81.0 

348.2 

41.3 

75.8 

214.4 

133.4 

68.4 

181.9 

117.6 

73.0 

157.1 

98.0 

70.8 

181.3 

86.9 

73.4 

220.7 

76.1 

76.0 

108.8 

68.7 

77.2 

82.0 

31.1 

78.6 

93-1 

37.9 

77.6 

98.O 

46.1 

70.2 

8.9 

93-1 . 

70.0 

18.7 

132.6 

68.0 

328.3 

106.5 

68.4 

340.4 

121.5 

67.9 

314.5 

95-5 

66.7 

315.8 

94.7 

67.1 

315.5 

117.1 

68.1 

321.3 

127.0 

69.0 

326.2 

131.6 

































Hertz were d igi t ized at the nine radial stations indicated in figure 1 1 6 at 
axial stations 36, 32, 28, and 24. As previously discussed, the angular 
patterns were first Fourier analyzed and then the lined duct radial modes were 
used to compute the radial mode contribution for an angular set. 

Table VIII contains the angular mode amplitude at each radial and axial 
station surveyed, for the(3,0)mode at 1000 Hertz. The pressure distribution 
which resulted from the survey will be shown in the next section. The columns 
in the tabulation give the modal amplitude, the phase of the mode referenced 
to the oscillator and the angular phase shift of the first minimum of the 
angular mode in question referred to zero degrees of azimuth. This is tabu- 
lated for the first five angular modes at the nine radial stations for each of 
the four axial locations. The 3rd angular mode, which corresponds to the mode 
set (3,0) , (3,1), ... (3,N) , which was simulated in the source section, is of 
the order of 20 decibels above the other modes which were spuriously excited. 
Closer examination of the table, however, reveals that even the unintentionally 
excited modes are present in a consistent manner in most cases. For instance, 
at station 36 , the amplitude variation across the duct in the first order 
angular mode set ((1,0), (l,l) ...) is uniform to within 3 decibels and the 
phase is uniform to within 20 degrees. Further, if that mode, which is only a 
spurious excitation from the (3,0) mode, is followed down the duct, it is seen 
that the location of the first angular minimum (indicated by the phase shift 
column) is consistent to within about 15 degrees, indicating that even the 
modes which were unintentionally excited, are responding and do propagate, and 
also that the instrumentation system is capable of measuring over a SPL range 
of greater than kO decibels. 

In order to determine the radial mode content in the 3rd angular set, the 
data listed in table VIII for each radial station were Fourier-Besse 1 analyzed, 
using the computed radial modes for the particular frequency and wall impedance 
boundary conditions. Table IX gives the result of this analysis for the first 
two radial modes at each of the four axial stations. 


TABLE IX 

RADIAL MODE ANALYSIS FOR 3RD ANGULAR MODE SET 


Configuration: U" Deep Liner with kO Rayl Felted Metal 
Frequency: 1000 Hertz 

Impedance Used in Computing Radial Modes: 1.22 - 0.65i 


p ^ 



RADIAL MODE NUMBER 



Axial 
Stat ion 


0 



1 


AMP 

Phase 

Phase 

AMP 

Phase 

Phase 

1 nches 

dB 

Deg. 

Shift Deg . 

dB 

Deg. 

Sh i ft Deg . 

36 

122.7 

118.1 

102.9 

113.6 

247.6 

103.4 

32 

226. 1 

186.8 

107.7 

108.9 

303.9 

108.1 

28 

1 10.6 

255.9 

114.3 

102.6 

17.1 

114.2 

2k 

101.6 

6.9 

120.9 

95.2 

111.7 

121 .0 
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Examination of this data points up an inherent problem in the analysis 
system which may have been avoided by more radial data points. Since only 
nine radial points were used, what would visually appear to be minor errors, 
ultimately results In 'noise' of the order of 10 decibels down from the pre- 
dominant radial mode. Thus, any radial mode which is more than 3-6 decibels 
below the predominant mode will be obscured by 'noise| generated from minor 
errors in the predominant mode measurement and analysis. In the case under 
consideration, the (3,1) mode is determined to be on the order of 8 decibels 
below the (3,0) mode at all stations, based on the experimental analysis tech 
nique. At the entrance to the lined section, a significant contribution of 
the (3,1) and even higher radial modes is expected, in order to match the hard 
wall (3,0) mode (which theoretically is 2.3 decibels lower at the inner wall 
than at the outer wall). However, since the (3,1) mode is, for al 1 practical 
purposes, cut off at 1000 Hertz (the computed attenuation rate is 23 decibels 
per duct width), it should be insignificant at station 32. Obviously, this is 
not the case, as the data shows. 


Now let us examine the data for the (3,0) mode and compare it with the 
measurements made previously as now shown in figure 117. The attenuation rate 
determined from this analysis for the (3,0) lined duct mode is 6.1 to 6.8 dB 
and the propagation constant is between 0.66 and 0.79 according to the weight 
given the phase data at station 24. The mean square values are 6.5 dB atten- 
uation, with 6 = 0.75. These values of complex propagation constant are also 
shown on figure 88 for comparison with the values determined previously. The 
median value is very close to the effective impedance value and, in effect, ( 
justifies use of the analysis technique for h i gher orde r modes , if the noise 
problem resulting from insufficient radial data points did not exist. 

As a further test of the data analysis system, a complex angular modal 
pattern was set up at 1000 Hertz with liner 1 installed in the test section. 
The modal pattern consisted of an equal summation of the first four angular 
modes at the transducer plane; that is, the pressure at the source was 


p = 1 + cosG + cos26 + cos 39 * 


( 81 ) 


After Fourier analysis of the angular patterns at the 9 radial stations, at 
each of the four axial stations, the radial Four ier-Besse 1 analysis was con- 
ducted. The results of this analysis for the first radial mode for each of 
the four angular modes is listed in table X. 

The most consistent method of determining the average modal attenuation 
rate and propagation constant from the data listed was to fit a least-square 
regression line through the four data points for SPL and phase for each mode. 
The least-square modal attenuations and 3 are listed in Table XI with the 
computed values for each mode using three impedance values, viz. that 
measured by the impedance tube, the modified two microphone method, and 
the effective impedance as described in this section. 
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TABLE X 

FIRST RADIAL MODE FOR FOUR EQUAL AMPLITUDE 
ANGULAR MODES 

FREQUENCY; 1000 Hz - LINER 1 “ NO FLOW 


- [ 

Axi al 
Stat i on 
Inches 

ANGULAR MODE NUMBER 

0 

1 

2 

3 

__ 

Amp 

dB 

Phase 

Deg. 

Amp 

dB 

Phase 

Deg. 

Phase 
Sh i ft 
Deg. 

Amp 

dB 

Phase 
Deg . 

Phase 
Sh i ft 
Deg. 

Amp 

dB 

Phase 
Deg . 

Phase 
Sh i ft 
Deg. 

36 

32 

28 

24 

1 12.6 
109.7 
104.0 
102.3 

I80.*t 

290.3 

29.9 

15*1.7 

114.3 

109.4 

106.4 
101 . 1 

106.7 
196.2 
301 .9 
97. *t 

128.4 

131 .0 

129.8 

125.2 

116,0 
1 10.0 
104.5 
101.9 

150.0 

248.2 

348.6 

92.4 

82.3 

88.8 

98.2 

87.5 

1 16.6 
109.3 
103.0 
97.1 

117.5 

187.5 
275.0 
3*+7 . 7 

26.1 

32.5 

40.0 

35.4 


TABLE XI 

COMPARISON OF EXPERIMENTAL AND THEORETICAL ATTENUATION 
RATES AND PROPAGATION CONSTANTS FOR FIRST FOUR ANGULAR MODES 

FREQUENCY: 1000 HERTZ - LINER 1 - NO-FLOW 


MODE 

NO. 

EXPERI Mi 
VALUE! 

ENTAL 

5 

IMPEDAh 

Two Mike 
Z= 1 .25 - 0.931 

ICE VALUE FOR CAL Cl 

Effective IMP 
Z=1 .48 - 0 . 88 i 

LATI0NS 

Impedance Tube 

Z= 1.65 - 0 .80 i 


Amp dB 

3 

Amp dB 

3 

Amp dB 

3 

Amp dB 

3 

1.0 

2.0 

3.0 

4.0 

3.77 

4.39 

*i.93 

6.68 

1 .034 
1.091 
0.832 
0.733 

3.85 

4.10 

4.78 

6.45 

1.100 

1.055 

0.958 

0.780 

3.60 

3.94 

4.78 

6.20 

1 .070 
1 .020 
0.914 
0.720 

3.40 

3.70 

4.54 

6.00 

1 .047 
0.998 
0.891 
0.680 


First, it is obvious that the trends exhibited by the 
match the theory in almost all respects. The attenuation 
number increases and the wavelength constant, 3, decreases 
mode number. The effective impedance was chosen as being 
measured data for the (3,0) mode in isolation. Once again 
impedance produces theoretical results which compare well 
data over all the modes. 


experimental data 
increases as mode 
with increasing 
the best fit to the 
, this value of 
with the measured 


Now summarizing the results discussed 
dimensional data acquisition Fou r ie r-Besse 1 
conclusions can be reached. These are (1) 
previously and the 4 angular mode situation 


to this point concerning the two- 
analysis system, two important 
the data on the (3,0) mode shown 
just discussed, further strengthen 
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the earlier conclusion that the eigenvalue theory of sound propagation in an 
annular duct is adequate, and (2) it is shown that mode superposition is a 
valid concept there is no evidence of mode cross-talk for these zero order 
radial modes. 


VI. 3. 1.4. Higher order radial modes . At this stage of the discussion, it is 
appropriate to examine the utility and accuracy of the two-dimensional 'measure- 
ment and analysis system for studying higher order radial modes. Liner 1 was 
tested at 1808 Hertz with an angular mode set-up for the (2,0), (2,1) ... etc. 
modes. The results from the angular Fourier analysis (i.e., the data used as 
input for the radial Fourier-Bessel analysis) are graphically shown in figure 
118. It would appear from visual inspection that the nine radial data points 
follow the expected trend. At the entrance to the lined duct section, it is 
evident that the (2 , 1 ) hard wall mode dominates. As the mode progresses down 
the lined section, the radial shape rapidly assumes a new form, as does the 
phase distribution. 

The results of the radial analysis of the nine data points are listed in 
table XII for experimental values of impedance determined from two methods. 


TABLE XI I 

RADIAL MODE ANALYSIS FOR COMBINED (2,0) & (2 , 1 ) MODES 
FREQUENCY : 1808 HERTZ - LINER 1 - NO FLOW 






ANALYSIS 

IMPEDANCE 



AXIAL 
STATION 
1 NCHES 

Impedance Tube: 

Z = 1.40 + 0.22 i 

Two- Mike: Z - 

1.22 + 0 .22 i 

RADIAL MODE ORDER 

RADIAL MODE ORDER 


0 

1 

0 

' 1 


Amp 

dB 

Phase 

Deg. 

Amp 

dB 

Phase 

Deg. 

Amp 

dB 

Phase 
Deg . 

Amp 

dB 

Phase 
Deg . 

36 

103.3 

158. 4 

122.7 

338.7 

104.0 

146.8 

122,6 

334.3 

32 

110.7 

211.9 

113.8 

10.1 

110.0 

205.8 

1 13 . a 

3.6 

28 

108.2 

4.8 

109.3 

135.6 

107. ^ 

358.9 

109.1 

127.5 

24 

IO6.3 

162.5 

104.9 

297.0 

105.5 

157.0 

104.6 

287.5 


It is immediately obvious that some problem, not previously discussed, exists 
^ the data * ^ e 9 a rd less of the impedance data chosen, the amplitude of the 
U,0; mode increases from station 36 to 32. If the data could be assumed to 
be absolutely accurate, this would be indicative of energy transfer between 
modes. It is emphasized that the analysis incorporates the same set of lined 
duct modes at each axial analysis point. There is no possibility of energy 
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conversion due to switching from hard wall modes to soft wall modes (in the 
anal ys i s , that is). 

The computed mode attenuations range from 2.38 to 2.48 dB for the (2,0) 
mode and 15.31 to 16.93 dB for the (2,1) mode. If only the last three data 
points for the (2,0) mode are considered, an attenuation on the order of 2,2 
dB results. However, the (2,1) mode attenuation resulting from these three 
data points is 4.5 dB; considerably less than the computed value. It is 
assumed that the problem discussed earlier, that of 'noise'* is indeed a 
reality. The only position where it would appear that the(2,1)mode is decay- 
ing much more rapidly than the (2,0) mode is near the liner entrance. This 
rapid decay then merges into the analysis ' noi se* before reaching station 32. 

Thus, it must be concluded, as suspected earlier, that the radial mode 
decomposition analysis method is in doubt, due to error in the measured data as 
a result of too few radial measurement points. Therefore, no further attempts 
to measure higher order radial mode decay will be presented. It is felt that 
the method is valid, as shown^ for the lower order modes, but that either a 
more accurate measurement/analysis scheme must be incorporated or else more 
radial data points must be acquired. If the number of radial data points was 
increased without decreasing the number of angular points, the analysis would 
become prohibitively time consuming. However, it is probable that a better 
balance between angular and radial data points could be achieved. As it now 
stands, the angular analysis is probably over-accurate with 128 points. 

Perhaps if the 128 x 9 matrix were divided such that the same number of data 
points were allocated for each half-cycle of the mode, whether angular or 
radial, the accuracy would be optimized. With this criterion, the number of 
data points for each half-cycle would be 


128 x 9 

2 x m x (n+1) 


(82) 


For a 2,1 mode, the number of angular data points would more appropriately be 
48 and the number of radial points would be 24 (arrived at by allocating 12 
data points to each half-cycle). However, no such system was possible at the 
time this meas urement /ana 1 ys i s problem was discovered. 


VI. 3. 1.5. Complex propagation constant measurements for non-local l y reacting 
.liners . F i nally, further useful data was provided by measurements of atten'ua- 
tion and propagation constant with the two foam liners in place. Figure 119 
is a plot of the sound pressure amplitude at the outer wall. This plot alone 
is a rather meaningless measure of the performance of a liner in terms of 
modal attenuation. It does, however, reveal rather interesting behavior. 


Noise can be defined *i n this context as being inherent variations of liner 
impedance plus digitizing errors and computational approximations; the last 
two components being very small in the no flow case. 
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SOUND PRESSURE IN DECIBELS 



The sound field appears to have a rather strong reflected wave component, 
since what appear to be damped standing waves are in evidence all the way back 
to the test section entrance. However, as a result of the impedance measure- 
ments on the liner, which displayed a considerable variability in impedance as 
a function of position (of the order of ±50%), it was surmised that the 
reflections were the result of local changes in liner impedance. 

The slope of the sound level decay curve is approximately 4 dB/duct width 
(for this one liner, the duct width was 4.75“). The computed value of attenu- 
ation (based on a locally reacting assumption) for the (0,0) mode ranged from 
8 to 14 dB/duct width according to the combination of values selected for 
average impedance of the inner and outer liners. 

Here, it must be concluded that the fact that the material is non- 
locally reacting dominates the attenuation character i st i cs of the duct, since 
the measurements were checked several times and it is believed that the re- 
sults presented are accurate. 

A further test was conducted using a commercially available acoustic foam 
on the outer wall only. The inner and outer SPL and phase are shown over an 
axial distance of 12 inches for 3 frequencies in figure 120. Detailed calcu- 
lations for this case were not made, however, the generalized calculations of 
figures 87 and 88 apply for a locally reacting liner. The impedance at 1000 
Hertz is, Z = 0.35 “ 0.60i. The impedance was uniform within ±20% over the 
five samples measured. The calculated modal attenuation from figure 88 is of 
the order of 10.5 dB and 3 is of the order of 1.45. The measured attenuation, 
although not a linear decay, is approximately 4-5 dB/duct width and 3 is 
approximately 0.67. Thus, the measured data for this liner are also not in 
agreement with locally reacting calculations and it must be concluded that the 
non-local ly reacting impedance effect is very significant for this particular 
geomet ry/f requency combination. Further, the effect is to significantly 
reduce the attenuation from that predicted by the theory. 

It must further be concluded that if non-local ly reacting liners are to 
be used, an analysis incorporating Scott's (ref. 53) assumptions and solution 
techniques must be used. 

This concludes the presentation on propagation constant and attenuation 
in soft-walled annular ducts without flow. Next, the discussion will center 
on mode shapes and radial pressure distribution in the lined duct. 


VI. 3. 2. Mode shapes and pressure distribution . 

The previous display of comparisons between experiment and theory as 
pertains to the complex eigenvalues for lined annular ducts has, in effect, 
provided the basic substantiation for the mode theory of sound propagation in 
ducts. However, one critical detail of the mode theory which has been alluded 
to in the discussion, especially concerning the calculation of power attenua- 
tion, is the spatial distribution of sound pressure and phase for a particular 
mode. In this sub-section, a number of mode shape and pressure distribution 
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measurements for each of the liners will be presented. Also, for many of 
these experiments, valid comparisons with theory are possible. 


VI. 3. 2.1. Liner 1 mode shapes . The liner 1 configuration provided a signifi- 
cant portion of the useful data on pressure distribution and mode shapes. The 
angular mode shapes are not of particular concern, since the presence of an 
attenuating liner, provided it has a uniform circumferential distribution of 
impedance, does not alter the angular pressure solutions. Because of the 
previous concentration on data for liner 1 at 1000 Hertz, the radial distribu- 
tion of sound pressure at that frequency will be presented first. 

Figure 121 is a plot of the measured radial sound pressure and phase at 
1000 Hertz. The measurements were made at one inch intervals beginning one 
inch upstream of the liner entrance and extending twelve inches into the lined 
section. A 1.5 dB reduction occurs between station 37 and 36 indicating a 
reflection at the interface between the hard wall and lined sections. The 
radial pressure distribution at the entrance to the lined section is that of a 
hard wall (3*0) mode, which theoretically has a 2.3 dB lower level at the inner 
wall as compared to the outer wall. As the pressure wave progresses down the 
lined duct, the shape gradually changes, both in pressure and phase. At 
station 24, the pressure distribution is 6 dB lower at the inner wall with a 
55° phase difference between the inner and outer walls. Examination of figure 
91, which is for the same set of physical conditions, would imply a signifi- 
cantly different wall impedance, from those presented before, in order to 
match the mode shape. At the 'effective' impedance value, the mode shape is 
calculated to have a 4.5 dB lower level at the inner wall and a 33 degree phase 
change between walls. Thus, some differences do exist between measurement and 
theory. Figure 122 shows the comparison between the theoretical and experi- 
mental pressure distribution at the entrance and three axial stations spaced 4 
inches apart, based on an assumed plane radial distribution at the liner 
entrance. The calculations are made at an impedance of 2 = 1.4 - 0.8i, which 
lies approximately half-way between the two-microphone and the impedance tube 
experimental impedance values. While generally the comparison could be con- 
sidered excellent, there are differences of up to 2 dB at station 24 in the 
absolute sound pressure level and up to 18 degrees absolute phase difference. 

It is emphasized that all the calculations are referenced to the entrance 
plane, thus effects of modal attenuation and wavelength constant differences, 
as well as differences in mode shape are involved. At station 24, the pressure 
distribution has decayed to the shape of the (3,0) mode, since all the higher 
order soft wall modes have attenuated to much lower levels. 

Figures 123 and 124 illustrate the behavior at two other frequencies, 
consistent with the data presented in figure 112. In figure 123 for 884 Hz, 
the comparison with theory, calculated at an Impedance of Z - 1 .90 - 0.85; , is 
excellent for both amplitude and phase, down to station 24, where the compari- 
son shows a significant disagreement in phase (of the order of 60 degrees 

difference) and approximately 3*5 dB difference at the inner wall. It is not 
known what to attribute this difference to. It is probably an experimental 
error, but it could also be related to a change in liner impedance between the 
liner entrance and station 24. The comparison at 1430 Hz generally follows 



HOLE NUMBER 



S33y030 N I 310NV 3SVHd 


O 




the shape trends both in amplitude and phase. It is significant to note that 
the amplitude is lower at the outer wall for both experiment and theory, in 
opposition to the tests at 884 and 1000 Hz. The cal cul cu 1 at i ons were made at 
a frequency of 1500 Hz, while the experiment was conducted at 1428 Hertz. This 
difference in wavelength would introduce approximately 15 degrees phase error 
between stations 36 and 24. The balance of the error is probably due to 
improper definition of the wall impedance due to spatial inhomogeneities of the 
1 i ne r mate rial. 

One further test of interest is the propagation of the (0,0) mode. Figure 
125 shows the pressure and phase distribution at six axial stations one inch 
apart beginning at station 37 and progressing to station 32. Although the 
pressure distribution is not initially plane at station 36, it rapidly assumes 
the shape of the calculated pressure distribution, which at station 32 , is 
predominantly the (0,0) soft wall mode. It is interesting to observe how the 
pressure distribution changes from a shape which is 1.5 dB lower at the outer 
wall to a shape which is 1.5 lower at the inner wall, in conformance with the 
soft wall (0,0) mode calculations. 


VI. 3. 2. 2. Mode shape and pressure distribution for liner 2 . Another series 
of radial plots were made in the hard wall section and in the lined section 
with liner 2 installed. This liner, as described earlier, was comprised of 
inner and outer cylinders consisting of acoustic material over 1/2 inch deep 
honeycomb with a rigid backing. The acoustic material was a nominal 20 Rayl 
material. Three higher order angular-radial mode combinations were excited 
in the hard wall source section at the response peak cor respond i ng to the 
damped cut-on frequency of the particular mode. The radial amplitude and 
phase plots are shown in figures 126 to 128 for the (2,1), (5,1), and (6,2) 
hard wall radial duct modes. No attempt has been made to provide a theoreti- 
cal comparison directly on the charts. These experimental studies do, 
however, provide a significant bit of data for future analytical studies, such 
as that of Zorumski (ref. 15) on sound propagation in a two section duct. The 
plots are placed onto separate sections of the page for ease in discrimination 
between individual plots, The soft wall section begins at station 35i for 
this particular liner. The parameters defining the test condition are given 
on each figure. 

Theoretical mode shapes in the lined section calculated at frequencies 
near the test frequencies and using average impedance acquired by the two- 
microphone method, are shown in figure 129. The frequency, impedance and com- 
puted complex propagation constant for each mode shown is given in the legend. 
The calculations show that the (2,1) mode should be only slightly attenuated 
at station 32, as is evidenced in the experimental plots (figure 126). The 
(5,1) mode should be attenuated by approximately 3"5 dB and this is also 
confirmed by the experiment (figure 127). The (6,2) mode also seems to follow 
the trend of the predicted attenuation of 4-6 dB (figure 128). 

Comparing figures 129 and 126, it is seen that the pressure distribution 
for the (2,1) condition very closely matches the computed (2,1) mode shape at 
station 32i, except near the minimum. From examination of figures 127 and 129, 
it is seen that the radial shape at station 32 for the (5,1) set seems to be a 
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axial stations 34 - 32 £ 

Figure 126. Radial SPL and Phase Angle Contours in Hard Wall Section and Liner 2 
Test Section for (2,1) Mode at 1800 Hertz 
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Figure .128. Radial SPL and Phase Angle Contours in Hard Wall Section and Liner 
2 Test Section for (6,2) Mode at 377B Hertz 
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a combination of the (5,0) and (5 , 1 ) soft wall modes, but is dominated by the 
(5,1) mode as it should be. From comparison of figures 128 and 129, the 
pressure distribution in the lined section for the (6,2) case does not really 
show any close resemblance to the (6,0), (6,1) or (6,2) soft wall modes pre- 
sented in figure 129(c). Figure 130 is a calculation of the pressure distri- 
bution at one duct width downstream for an assumed radial plane wave at the 
liner entrance for a 6th angular mode variation. As expected, the shape has 
decayed to the calculated shape of the (6,2) mode, since the other two cut-on 
modes have attenuations of 12 dB greater than the (6,2) mode. Thus, the 
experimental behavior at the (6,2) mode cut-on frequency seems to be dominated 
by other modes. 

As was demonstrated in figure 117, the attenuation derived from the two- 
dimensional Fourier-Bessel analysis was quite accurate for the zeroth order 
radial modes. Table VIII listed the Fourier analysis results associated with 
the summary data in figure 117* Figure 13l(a)-(d) shows the angular sound 
pressure distribution in liner 1 at the nine radial stations for each axial 
location. These data are presented as an illustration of the consistency of 
the data acquisition and analysis system for analyzing the zeroth order radial 
modes. Figure 132 is included as an aid in estimating the effect of axial 
distance on the amplitude of the (3,0) mode as it propagates down the liner 1 
conf i gurat i on . 

Because of the exceptionally good correlation with theory shown for the 
case of four equal amplitude angular modes propagating simultaneously, an 
example of the angular distribution of SPL and phase at the four axial sta- 
tions is shown in figure 133. These angular plots were obtained at radial 
station 6 as shown in figure 116,. The angular patterns at each axial station 
provide further insight into the effect of higher order mode attenuation. The 
SPL and phase pattern at station 36 is very complex due to addition of the 
four equal amplitude angular modes. But, at station 2*+, the (1,0) mode has 
decayed 2.5 dB more than the (0,0) mode, the (2,0) mode has decayed 5.0 dB 
more than the (0,0) mode and the (3,0) mode has decayed 11.6 dB more than the 
(0,0) mode. Thus, the amplitude of the (0,0) mode now dominates, as evidenced 
by the reasonably flat angular distribution of SPL. However, it is also 
evident that a strong component of the ( 1 , 0 ) mode still exists by examination 
of the phase plot. 

Finally, in the interest of providing original data for future use, the 
radial mode shapes for liners h and 5 at several axial stations are presented 
in figures 13*+ and 135 (a)- (c). Complete documentation concerning these 
measurements is given on the figures. 


VI. 3. 3. Effect of flow velocity in lined ducts. 


The most effective data for illustrating the effects of velocity were 
obtained in the liner 1 configuration. Two sets of data will be presented and 
in some cases these will be compared with theory. Figure 136 contains the 
inner and outer wall SPL and phase measurements at 1000 Hertz for the (0,0) 
and (3,0) modes at 50 fps . Figure 137 shows the radial pressure distribution 
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Figure 130. Pressure Distribution Across Annulus for Sixth Angular Mode Set at 3750 Hertz 
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Figure 131. Angular SPL and Phase Angle Distributions in Liner 1 at 1000 Hertz without Flow 
for Two-Dimensional ' Fourier-Bessel' Analysis - (a) Axial Position: 36 Inches 
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Figure 132. Comparison of Outer Wall Angular SPL Distribution for Four Axial Stations (see Figure 131) 
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Figure 137 . Radial SPL and Phase Angle Contours for the (0,0) and (3,0) Modes in Liner 1 at 1000 Hertz 
with 50 fps Flow Velocity 


'■j 





at four axial stations spaced 4 inches apart for the two conditions shown in 
figure 136. Figure 1 38 illustrates the attenuation of the (3,0) mode at 150 
fps for 3 frequencies and figure 139 illustrates the radial pressure distri- 
bution at two of those frequencies. 

From figure 136(a), the (0,0) mode attenuation is determined to be 
approximately 4.0 dB, which is essentially that obtained for the no-flow case. 
The propagation constant, 3, has reduced from 1 .08 for the no-flow case to .99 
at 50 fps. Calculations were not made for comparison with the (0,0) mode case 
at 50 fps. Figure 136(b) shows the (3,0) mode attenuation and propagating 

phase slopes. The average attenuation is 5.4 dB and 3 is 0.551. This compares 
with a computed value of attenuation of 5.6 to 5.8 dB. 

The pressure distributions presented in figures 137(a) and (b) are inter- 
esting for illustrating the effect of turbulence on the measurement evidenced 
by small irregular variations in the amplitude and phase. 

In figure 138(b), the (3,0) mode average attenuation is 5.7 dB and the 
average 3 is 0.511. Calculated values are on the order of 4.5 to 5.0 dB 
attenuation and 0.55 to 0.65 for 3. Good values of impedance were not ob- 
tained on liner 1 at 150 fps due to a drastic reduction in driver output with 
increasing flow velocity. It was later discovered that this could be effec- 
tively remedied by approaching an overdriven condition on the source trans- 
ducers, however, there was insufficient time at that stage to rerun the tests. 

The pressure distributions are shown in figure 139(a) and (b). These are 
relatively unexciting, since there appears to be no major effect of flow, as 
predi cted. 

Finally, data complementary to figures 1 3 1 (a) - (d ) are shown in figures 
l40(a)-(d) for 50 fps and I4l(a)-(d) for 150 fps. In this case, the effect of 
turbulence on the measurements is quite obvious at station 24. The pertinent 
data from the two-dimensional ' Fourier-Bessel 1 analysis are given in the 
following table for both flow velocities. 

TABLE XIII 

TWO-DIMENSIONAL ANALYSIS SHOWING EFFECT 
OF AIR FLOW VELOCITY ON (3,0) MODE 

CONFIGURATION: LINER 3 

FREQUENCY: 1000 HERTZ 




FLOW VELOCITY 



AXIAL 

NO 

FLOW . 

50 FPS 

150 

FPS 

STAT 1 ON 

Amp 

Phase 

Amp 

Phase 

Amp 

Phase 

INCHES 

dB 

Deg. 

dB 

Deg. 

dB 

Deg. 

36 

122.7 

118.1 

127.0 

73.5 

126.5 

49.7 

32 

116.1 

186.8 

121.3 

127.6 

120.9 

97. 

28 

110.6 

255.9 

115.3 

215.5 

116.6 

137.0 

24_ 

101.6 

6^.9. 

108.4 

305.6 

109-3 

218.1 
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(a) 884 Hertz 



fb) 1000 Hertz 



Figure 1 38 . Sound Pressure Level Decay and Phase Angle Change for the (3,0) 
Mode in Liner 1 for Various Frequencies at 150 fps Flow Velocity 
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^ -°9 DISTANCE ACROSS DUCT - INCHES 

(b) 1430 Hertz 


Figure 139. Radial SPL and Phase Angle Contours for the (3,0) Mode in Liner 1 
at 884 and 1430 Hertz with 150 fps Flow Velocity 
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Figure U4O. Cont. 
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An analysis of the data in tab 1 e XI 1 1 using a least-squares linear regres 
si on gives the following average values of attenuation and 6 for each of the 
velocity conditions. 


TABLE XIV 


EFFECT OF AIR FLOW VELOCITY ON COMPLEX PROPAGATION CONSTANT 


VELOCITY 

FPS 

USING DATA 
FROM PREVIOUS TABLE 

USING DATA FROM FIGURES 
127 (b) , 152(a), and 15Mb) 

Attenuat i on 
dB 

P ropagat i on 
Constant 

3 

Attenuation 

dB 

P ropagat ion 
Constant 

3 

0 

7. 10 

.768 

6.2 

.719 

50 

6.37 

.735 

5.b 

.551 

150 

5.76 

.508 

5.7 

.511 


Included in the last two columns of this tabulation are the attenuation and 3 
values determined from the average values between inner and outer wall SPL 
levels from figures 111(b), 136(a) and 138(b). The comparison between these 
two sets of data is reasonable, but not good. Both sets of data do, however, 
show a decrease in both attenuation and 3. The 3 decrease is to be expected, 
but the attenuation decrease is dependent on many factors and an increase 
could also occur, provided the physical characteristics such as impedance and 
frequency satisfy the conditions for an increase. 


VII, THE MEASUREMENT OF HIGHER ORDER CIRCUMFERENTIAL MODE 
RADIATION IMPEDANCE FOR AN ANNULAR DUCT WITH FLOW 


V I 1 . 1 . I nt roduct i on 


The theoretical treatment of the generation of sound in, transmission 
through, and radiation from bypass air ducts of the fan-jet engine family can 
be separated into three distinct but partially coupled problems, as discussed. 
The problem of sound radiation from an engine duct has not been considered in 
detail in the literature using this superposition approach, but approximate 
treatments have discussed the radiated field for simplified geometry, without 
the presence of the duct air flow (e.g. Morfey, ref. 8, Sofrin and Tyler, ref. 
7) . If the superposition approach is to be utilized in predicting sound 
radiation from a fan -jet engine, it is necessary to know the acoustic radia- 
tion impedance at the duct termination, since this impedance serves as the 
coupling link between the duct and the radiated field. 
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There have been basic calculations of the radiation impedance (for the 
plane and higher order modes) of a flanged annulus by Morfey (ref. 26), Doak 
(ref. 24), and Haddle (ref, 25), but the effects of axial flow through the 
duct and a realistic termination were not included in these studies. One 
other study by Lansing, Drischler and Pusey (ref, 14) implicitly calculated 
the radiation impedance in determining the radiated field with a realistic 
termination. This work, which extended the Wiener-Hopf type solutions of 
Levine and Schwinger (ref. 11) and Carrier (ref. 12), proceeds directly from a 
source distribution to the radiated field. The radiation impedance can be 
calculated using this procedure. 

There have been few reports on measurements of radiation impedance to 
serve as verification of the theoretical predictions. Plane mode radiation 
impedance was measured by Mechel, Mertens and Schilz (refs, 54 and 55). They 
demonstrated that for frequencies below the first higher order mode cut-on 
frequency the measured plane mode radiation impedance for a circular piston 
compares favorably with calculated values. They also measured the effect of 
axial velocity on the plane mode radiation impedance for a circular aperture. 
The effect of velocity was to reduce radiation resistance up to a frequency of 
ka - 1.7, which was the frequency limit of the tests. The velocity effect on 
reactance is more difficult to characterize, but In general, at very low fre- 
quencies (below ka ■ 1) the reactance was reduced; conversely the peak value 
in the reactance curve increased with flow velocity. 

The major shortcoming with the available data on radiation impedance Is 
the complete lack of measurements for higher order modes. Further, no known 
data exists for annular ducts (or flanged annular pistons) in any frequency 
range. However, utilizing the facility described previously In this report, a 
technique for measuring higher order mode radiation impedance has been devel- 
oped and it has been possible to check the predictions for no-flow radiation 
impedance. The technique requires rather lengthy data reduction methods, 
which are outlined in the following section. Data for higher order mode 
radiation impedance are discussed and compared with available theoretical pre- 
dictions. Also, effects of flow on the radiation Impedance are discussed. 


VI 1.2. Experimental Technique for Determining 
Modal Radiation Impedance 


The measurement of modal radiation impedance necessitated the development 
of several new data analysis methods because of 

o restrictions on microphone degrees-of-f reedom in the test section, 
o complexity of the sound pressure field In the duct, and 
o the added complexity of measurement due to flow. 
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VII. 2.1. Sound field mapping 


The philosophy behind measuring only the axial distribution of sound 
pressure in a plane wave impedance tube is that the sound pressure distribution 
in the plane normal to the tube axis is known. Once the sound field becomes 
complicated enough that the sound field in the plane of the tube or duct is 
unknown, more detailed measurements are required to characterize the wave com- 
ponents. It is certainly valid to assume that radiation impedance of a higher 
crder mode field can be calculated from the axial sound pressure distribution 
of that field, provided enough is known about the wave field to properly 
separate the components of all modes present. The measurements necessary to 
complete this separation are in essence a three-dimensional mapping of the 
duct sound field. Obviously, a continuous volume mapping of the sound field 
is not necessary as the boundary conditions remain uniform along the tube 
axis. Thus, only measurements at properly spaced points are required. 

As discussed earlier, it was not feasible to mount a probe microphone 
projecting into the test section from the end to permit a continuous axial 
survey; therefore, on the outer wall of the test section, microphone mounting 
holes were located on a one-inch spacing down the axis of the duct. With this 
system, it is possible to obtain a fairly complete mapping of sound pressure 
in a plane normal to the axis of the annulus at the one inch spacing locations. 
The system does not, however, provide the usual continuous standing wave plots 
down the duct, with the required resolution for direct measurement of the 
standing wave ratio and location of the first minimum. 


VI 1.2.2. Duct sound field 


A first impression is that it is particularly expedient in the deter- 
mination of radiation impedance to excite a single higher order duct mode in iso 
lation. It has previously been shown that this capability exists for angular 
modes at frequencies below the first higher order radial mode cut-on frequency. 
It is not always possible, however, to achieve the degree of isolation required 
for a particular mode and in some instances mode isolation is less than 10 dB. 

In this instance, it is impossible to determine the radiation impedance of a 
single mode, because of the multiple mode axial standing wave patterns present 
in the duct; thus, the effects of the other modes which are spuriously excited 
must be analytically removed. To accomplish this, a complex Fourier analysis 
technique similar to one described in section III was developed. One of the 
interesting facets of this analysis technique is that several modes could be 
excited simultaneously and the data, after Fourier analysis, could be used to 
determine the radiation impedance of each mode. This method was ultimately 
incorporated, because of the significant time savings in conducting the tests. 

It was verified, however, that this method does degrade the accuracy of the 
results to a minor extent. 
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VII. 2. 3. Measurements in flow 


Two major difficulties were encountered in making measurements in 
flow. The first involved reducing the random noise component. This was 
accomplished with two-stage filtering. First, the microphone signal was 
filtered through a 1/3 octave filter set. The output of this stage was then 
amplified and filtered through a 2 Hertz narrow-band tracking filter. The 
resultant filtered signal was very stable and had a signal to noise ratio of 
greater than 20 dB at 150 fps based on a source sound level of 120 dB. 

The second difficulty was evidenced when trying to maintain a constant 
value of wave number, k = 27Tf/c, during the testing, with the inherent temper- 
ature changes associated with delivering air from the compressor and storage 
tanks. To overcome this difficulty, as briefly discussed previously, an 
analog servo control was developed using temperature as the control signal. 
This servo control maintained a constant value of f req uency-to-sound speed 
ratio for temperature variations of ±20°F. The accuracy was better than four 
decimal places in impressed wave number, k. 


VI 1.3. Data Analysis Methods for Determining Radiation Impedance 


The problems of measuring the necessary data for calculation of 
radiation impedance have been outlined. The technique chosen for determining 
radiation impedance used the axial standing wave pattern in the duct to 
determine the reflection coefficient at the duct termination. The impedance 
was directly determined from the reflection coefficient and the propagation 
constant, as will be subsequently shown. However, to obtain the appropriate 
SPL amplitude and phase for a particular mode at each axial measuring station 
it was necessary to Fourier analyze the data in the circumferential direction. 
Figures 142(a) and 142(b) are typical of the angular (circumferential) SPL 
distributions where single mode setups were attempted. In figure 142(a) for a 

(3.0) mode at 1000 Hertz, from a cursory examination, isolation from the lower 
order cut-on modes has been achieved. But in figure 142(b) for the (1,0) mode, 
the distortion in shape and phase is evident. Use of the peak amplitude of 
the measured angular pressure pattern would not be sufficient to determine the 

(1.0) mode standing wave ratio. Therefore, the data must be reduced to its 
Fourier spatial components. 


VII. 3.1, Fourier analysis of angular pressure distribution 

For a pure standing wave, the phase must change by 180 degrees on 
passing through a nodal point. Any deviation from the 180° phase change is 
indicative of either a spinning wave component of the same mode shape or other 
mode contributions. In either case, the pressure wave must be treated as a 
complex number. 


182 


SOUND PRESSURE LEVEL - DECIBELS 


(a) TYPICAL SOUND PRESSURE DISTRIBUTION IN TEST SECTION FOR (3,0) MODE SIMULATION AT 1000 HERTZ 



(b) TYPICAL SOUND PRESSURE DISTRIBUTION IN TEST SECTION FOR (1,0) MODE SIMULATION AT 1000 HERTZ 



Figure 142. Examples of Angular SPL and Phase Distributions in Hard Wall Test 
Sect i on 
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measured plane mode axial standing wave characteristics. The central idea in 
the higher order impedance technique is to use the peak value of the modal 
pressure pattern at each axial station to determine the axial pressure wave. 
Then using an analogy to the plane mode impedance tube methods, the impedance 
is calculated from the reflection coefficient and an additional factor relat 
ing to the propagation constant, $. Also, if the duct contains soft walls, 
the modal attenuation constant must be included in the duct analysis. 

The modal standing wave in the axial direction is determined from the 
usual assumption of a forward and backward traveling wave set and uses a com- 
plex least-square curve fit to find the amplitude of each component. The 
assumed form of amplitude of the nth mo de pressure wave is 


P n (z) - 


A e 
n 


k <«n 


i 6^)z 

M n 


-k(a 


B n e 


i3 b )z 

n 


( 91 ) 


where 


a 


n 

*n 

b 

i 

n 


is the modal attenuation in the forward direction 

is the modal propagation constant in the forward direction 

is the reflected wave modal attenuation 
is the reflected wave modal propagation constant 
is the amplitude of the incident wave 
is the amplitude of the reflected wave. 


The least square curve fit criteria were applied, minimizing with respect 
to the real and imaginary parts of both A n and B„, as was done In the sP'nnmg/ 
standing wave analysis. In the general case, it would be very difficult to 
determine the correct values of the set (oj, 6 f n , e$, without some P rl ° r 
knowledge of the expected values. The tests to be discussed later were con- 
ducted in a hard walled annular duct and therefore the values of a- and a n 
were, for all practical purposes, negligible. The values of $ n and 3 n are 
reasonably well approximated by assuming plug flow, using the average value^ot 
flow velocity, in the hard wall, mean flow wave equation. The exact expres- 
sion for $ for a plug-flow situation can be found, for instance, from Mungur 
and Plumb lee (ref. 6). That is 


K 




|Ct n - 


[ - M ± yV + (1-M 2 ) (l-k r 2 / (kL) 2 ) 
1-M 2 


( 92 ) 
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and 


„b Q b . b 

K =3 - i a 

n K n n 


-M ± /m 2 + (1-M 2 ) (1-k 2 /(kL) 2 ) 


1 “M 


(93) 


where H is the mean Mach number and k is the eigenvalue. The values of 
propagation constant determined from above, using the theoretical eigenvalue, 
were used in the minimization process for determining the amplitudes of the 
incident and reflected waves. 

After determination of the wave amplitudes, the reflection coefficient is 
determined in the same manner as for the plane mode, i.e., the ratio of the 
reflected wave to the incident wave, 


R 

n 


B /A 
n n 


(94) 


VI 1.3.4. Acoustic impedance calculation 

The calculation of radiation impedance in the presence of axial flow 
requires inclusion of the flow velocity effect on axial acoustic velocity. The 
momentum equation is, after assuming harmonic motion, 


3v _ ? kv _ 1 9p 

3z M p c M 8z (95) 


The solution yields the usual homogeneous and inhomogeneous parts. The 
homogeneous solution is of the form, Ce'kz/M, whIch is a con vected solution in 
the form of a frozen pattern. Because of this and its behavior at zero 
velocity, it is eliminated from consideration. Thus, the solution for acoustic 
velocity in the nth mode, is 


V as 


_L 

pc 


A X 


. f, k K 

1 a n) e 




1 " 


M(B^ - io f ) 
n n 


B ,.b . . b> * k( “n + ! ^n) z 
B n (3 n ia n )e 

1 + M(3^ - ic£) 


1.(96) 


Of course, the nt pressure mode produces velocity components in all 
other modes, but these components have not been considered in the present 
analysis. | n the data presentation which follows, only angular modes coupled 
with the zero order radial mode, have been considered. In a detailed exami- 
nation of the experimental data it was not possible to detect the presence of 
higher order radial modes at the termination, in the frequency range considered 
for these tests (i.e. below 1600 Hertz). 


187 



At the end, z = 0, the normalized acoustic modal radiation impedance is 
expressed as the ratio of modal pressure to modal velocity in the n™ mode, 


1 + R 


rad pcv 


<?n ‘ ia n) R n 


- J«n) 

1 - M(B„ - ia„T 1 + ‘ ia n) 


(97) 


If the duct walls are non-attenuating and the flow velocity is zero, the 
above equation reduces to the usual form, 


- 1 + R n 

rad n B n° " R n 


(98) 


The general expression in equation (98) was programmed for the determina- 
tion of modal termination impedance, utilizing the values of reflection coeffi 

cient. R determined from the least-square minimization process, 
n 


VII. 3. 5. Optimization of propagation constant 

As discussed earlier in this section, the theoretical value of 
propagation constant, 3, was used to determine the axial pressure pattern via 
the least-square curve fit. In some cases, it is evident that a longer or 
shorter wave-length, corresponding to different values of 3 , would yield a 
closer fit to the measured data. As a means of determining the accuracy of 
the curve fit, and the appropriateness of the theoretical propagation con- 
stants, the ’’standard error of the estimate" (the square root of the 'sample 
variance") was computed. This computation used the sum of the square of the 
differences between the measured amplitude and the computed amplitude. Phase 
differences were not included in the estimate, but could easily be incorpo- 
rated. The expression for the error is, 



(99) 


where L (x-) is the measured SPL at the i ^ axial station and L n (x.) is the 
'least-square' SPL at the i ^ axial station. 
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In attempting to determine the best value of 3 to fit the measured data, 
the curve-fit routine was used repetitively and values of 3 ^ and 3 b , corre- 
sponding to values computed from equation 97 for incremental values’of k 
were used. The error, E, was examined for each computation and the pair r of 3 
values which produced the minimum error was taken as the optimum set. 

This procedure appears to be inherently sensitive to experimental error, 
since one bad data point could result in a minimized value of error at a cut- 
on frequency (or wavelength) significantly removed from the actual value. 

After examining the impedance results obtained from the optimization process, 
the impedance data calculated using theoretical cut-on frequency and theoreti- 
cal values of 3 were chosen as being most representative. However, in the 
data presentation which follows, the computed axial standing wave pattern 
corresponding to both the theoretical and minimized values of 3 are shown. 
Also, the values of impedance corresponding to both values of 3 will be given 
in some cases. 


VII. A. Experimental Values of Radiation Impedance 
and Comparison with Theory 


In the experiments conducted, radiation impedance was experimentally de- 
termined at 5 frequencies. This obviously does not yield enough data for 
experimental verification of the detailed shape of the computed radiation im- 
pedance curves, but the data is sufficient for verifying the validity of the 
technique and to confirm the trends predicted by the theory. 


First, it is of interest to examine the axial wave forms determined from 
the measured data. These are shown in figures 143(a) through (143(f). Figure 
143(a) shows data. measured at 1000 Hertz (kR = 4.4205) for 32 measurement 
locations. In this particular case for the (1,0) mode, the 
the estimate, E, is only 0,25 dB, which is considered to be 
accurate to determine the radiation impedance within 1-2%. 

stance, the 3 at which minimum error, E, occurs corresponds __ . v, 

value of 6 of 0.953. In figures 143(b), 143(c) and 143(d) the measured’ data 
and calculated sound pressure patterns in the duct are shown for the plane 
(2,0) and (3,0) modes for 150 fps flow velocity. ? 


standard error of 
sufficiently 
Also in this in- 
to the theoret i cal 


In this series of tests, fewer axial measuring points were used in order 
to reduce. time required for the tests. This series of figures depicts the 
increase in wavelength with increasing mode number at constant frequency 
(equivalent. to a decrease in 3 ). In addition to the axial pressure pattern at 
the theoretical values of 3 , the pressure distribution calculated at the mini- 
mum error 3 is shown in figure 143(d) (dashed line). This difference in 
propagation constant results in a 25% difference in impedance magnitude. No 
firm conclusion was reached on whether it would be more accurate to use the 
impedance determined at minimum error or to use the values calculated at the 
theoretical propagation constant. However, all results presented will corre- 
spond to values at theoretical 3, since the comparison of cut-on frequency, 
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AXIAL PRESSURE DISTRIBUTION IN ANNULAR DUCT FOR (1 ,0) MOOT TEAK AMPLITUDE, 
DRIVED ROM LEAST SQUARE FIT TO MEASLRED DATA 


FREQUENCY PARAMETER: kR * 4.4205, kL « 1 .9194 


EIGENVALUE: 
FLOW VELOCITY: 


Lj * 1 .2922, k L - 0.5411 

IO fp. 


REFLECTION COEFFICIENT: -•( 

RADIATION IMPEDANCE: .« 

PROPAGATION CONSTANTS: 


THEORETICAL B 
-.0970 -.20431 


MINIMUM ERROR 6 
; -.0970 -.2043i 
: .9083 -.35141 



12 I* » 

DISTANCE (INCHES) 


/ U \ AXIAL PRESSURE DISTRIBUTION IN ANNULAR DUCT FOR (0,0) MOOE PEAK AMPLITUDE, 
V D / DERIVED FROM LEAST SQUARE FIT TO MEASURED DATA 


FREQUENCY PARAMETER: kR o - 4.4144, kL - 1 .944* 

EIGENVALUE: - 0.0000, ^ • 0.0000 

FLOW VELOCITY: 150 fp. 

THEORETICAL « MINIMUM ERROR I 

REFLECTION COEFFICIENT; -.0771 -,40«5t ; -.0771 -.40451 

RADIATION IMPEDANCE: .6234 +.41571 ; .6234 +.61571 

PROPAGATION CONSTANTS: .8812 B* - -1 . 154 , » f , - .8*12 4^ - -1.15* 


\ \ •' \ f \ / \ 

\/ V 'V \ 


4 ' 6 12 16 

DISTANa (INCHES) 


( f N AXIAL PRESSURE DISTRIBUTION IN ANNULAR DUCT FOR (2,0) MOOE PEAK AMPLITUDE, 
V v 'p DERIVED FROM LEAST SQUARE FIT TO MEASLRED DATA 


FREQUENCY PARAMETER: kR = 4.4146, kL - 1.9169 


EIGENVALUE: S( ‘ J -57I4, 1 

FLOW VELOCITY: 150 fp» 

THEORETICAL 

REFLECTION COEFFICIENT: -.1012 -.3576 

RADIATION IMPEDANCE: .7999 + .70471 

PROPAGATION CONSTANTS: Bj - .4944 8, ■ - .9690 


minimum error a 

-.1110 -.35381 
.7880 +.68491 

= .6961 8^ = -.9707 



( fj \ AXIAL PRESSURE DISTRIBUTION IN ANNULAR DUCT FOR (3,0) MOOE PEAK AMPLITUDE, 
DERIVED FROM LEAST SQUARE FIT TO MEASL* ED DATA 


FREQUENCY PARAMETER: kR o * 4.146, kl - 1.9149 

EIGENVALUE: Si *3.8268, L l - 1.64)4 

FLOW VELOCITY: 150 Fp* 

THEORETICAL 8 MINIML 

REFLECTION COEFFICIENT: -.3989 -.58401 ; -.4959 

RADIATION IMPEDANCE: .3357 +1.0871 3627 

PROPAGATION CONSTANTS: 6, - .3843 B* ■= -.6589 , 8, = .3951 


MINIMUM ERR (OR 8 
-.4959 -.4822 
.3627 +.81561 

F „„ ,b . .« 



DISTANCE (INCHES) 


DISTANa (INCHES) 


( « \ AXIAL PRESSURE DISTRIBUTION IN ANNULAR DUCT FOR 11,0) MODE PEAK AMPLITUDE, 
' C / DFIIVFD FROM LEAST SQUARE FIT TO MEASURED DATA 


FREQUENCY PARAMETER: kR o = 4.4154, kL = 1 91 72 
EIGENVALUE: Si - I -2W2, X L » 0.5411 


EIGENVALUE: 
,FLOW VELOCITY 


50 fp. 

THEORETICAL 8 


REFLECTION COEFFICIENT: -.1408 -.IB491 

RADIATION IMPEDANCE: .7421 29101 

PROPAGATION CONSTANTS: 8* = 9137 8^ * -.1004 


MINIMUM ERROR 8 


NOT 

CALCULATED 



AXIAL PRESSURE DISTRIBUTION IN ANNULAR DUCT FOR (1,0) MODE PEAK AMPLITUDE, 
' 1 • DERIVED FROM If AST SQUARE FIT TO MEASURED DATA 


FREQUENCY PARAMETER: kR ■= 4.4146, kL - 1 .9169 
EIGENVALUE: Nl o * ' ’ m ‘ l ‘ X L * 0 - M ' ' 

FLOW VELOCITY: 150 Fp. 


REFLECTION COEFFICIENT: -.1015 -.23241 

RADIATION IMPEDANCE: .7779 . 39081 

PROPAGATION CONSTANTS: » f , * .8374 B^ * -.lt!2 


MINIMUM ERROR B 


NOT 

CALCULATED 


K_.r' 


DISTANCE - INCHES 


OISTANW - INCHES 


Figure 143- Examples of Axial Standing Waves in Annular Duct Showing. Effects of 
Mode Number and Velocity 
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which implicitly validates the theoretical 3, was so exceptionally good. It 
must be remembered that the amplitudes shown in the previous series of figures 
at each axial station were derived from the complex Fourier analysis. Thus, 
there are many manipulations performed on the data, before arriving at the 
calculated reflection coefficient and radiation impedance. It would certainly 
not be unusual to expect an error of the order of ±1 dB as a result of all the 
data handling. Thus, the value assigned to 3 as a result of the minimization 
process is possibly dominated by the random, but small, errors introduced by 
the data analysis. This reasoning predominated in selecting the theoretical 3 
as the value for use in determination of radiation impedance. 

The effect of flow on the axial standing wave patterns is shown in 
figures 143(a), 143(e) and 143(f). The first thing to note is the calculated 
effect on the wavelength constants, 3^ and 3 b . Without flow, 3^ = 0.953 and 
3 i 7 °»953, but at 150 fps the valles of ^wavelength constant are 3^ = .837 
and 3 b = -1.112. Next, examination of the standing wave pattern does not re- 
veal any significant differences, but the computed reflection coefficients and 
radiation impedances do reveal a trend. The change in the reactance is most 
significant, changing from -0.35 at zero flow to +0.39 at 150 fps. It is felt 
that the data at 50 fps have been significantly influenced by the data at 
stations 12, 18 and 27. These deviations from the standing wave curve, deter- 
mined by least-squares, result in a sample variance of 1.1 dB, while at no 
flow the variance was 0.24 dB and at 150 fps the variance is 0,84 dB. 

All the data for radiation impedance is summarized in table XV. This 
comprises three test series (i.e., 0, 50 fps and 150 fps flow velocity) for 
five modes, each at five frequencies, where possible. It was not, however, 
possible to determine the radiation impedance of modes below cut-on frequency. 
Thus the number of values decreases with increasing frequency. The data pre- 
sented in the fol lowing section are plotted against a non-dimensional fre- 
quency which is referenced to the mode cut-on frequency in the nomenclature 
suggested by Morfey (ref. 26). Thus, this parameter is included in the tabu- 
lation for each velocity, indicating the minor change in the eigenvalue, k 
with velocity. The most informative data in the tabulation is the sample r * 
variance, which is a good indication of the deviation of the data from the 
computed axial standing wave. This variance varies from a low of 0.25 dB to a 
high of 18.1 dB, but with the majority of data below 3 dB variance and only 
two error values above 5 dB. 

Measured values of radiation resistance and reactance are compared with 
the values calculated from Haddle's (ref. 25) method, for a flanged annulus, 
for the (1,0), (2,0) and (3,0) modes in figure 144. The calculations, using 
Haddle s program were made for a hub-tip ratio of 0,566, corresponding to the 
experimental facility dimensions. 

While the measured resistance values do not duplicate the predicted 
values, the agreement is close and differences are attributed (l) to diffei — 
ences between the physical set-up and the theoretical model and (2) to 
experimental error. The theoretical model assumes a flanged annulus, while 
the facility termination is unflanged. 
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TABLE XV 


RADIATION IMPEDANCE CALCULATED AT THEORETICAL ft BY I^AST-SQUARES METHOD 


APPROXIMATE 

NORM 

FREQ, Hz 

FREQ,k 


346.6 1.9491 

566.7 3.1823 

845.8 4.7468 

1000.0 5.6144 

1350.0 7.5853 

346.6 1.9^91 

566.7 3.1823 

845.8 4.7468 

1004.7 5.5838 

1350.0 7.5853 

845.8 4.7468 

1000.0 5.618 

1350.0 7.5853 

1000.0 5.6441 

1350.0 7.5853 

1350.0 ■ 7.5853 


kR -k R 

O TO 


SAMPLE 

ARIANCE 


± .9390 
± .9530 


P f P° 


kR -k R 
o r o 


-1.047 

-1.0468 

-1.0469 

-1.0469 

-1.047 

- .5940 

- .9052 

- .9859 
-1.0066 
-1.0236 

- .7761. 

- .8599 

- .9507 

- .5463 

- .8177 

- .5781 


.2877 • 

.4881 . 
.5035 • 
.7084 . 

1.200 
- .1218 . 
.^363 . 
.9249 . 

. .7439 • 
.3559 • 
.4505 1. 
1.071 . 


SAMPLE 

VARIANCE 

dB 

P f 

P b 

kR -k R 
0 r 0 

0.28 

.8096 

-1.1562 

1.54 

0.53 

.8817 

-1.1550 

2.51 

1.97 

.8814 

-1.1555 

3.74 

1.06 

.8812 

-1.1558 

4.40 

3-19 

.8811 

-1.1560 

5.94 

1.41 

.4316 

- .7071 

.22 

0.59 

.7404 

-1.0137 

1.19 

1.72 

.8204 

-1.0946 

2.43 

1.11 

.8374 

-1.112 

3.09 

4.75 

.6777 

-1.1326 

4.68 

1.67 

.6114 

- .6655 

1.16 

0.91 

.6944 

- .^690 

1.87 




3.43 

O.96 

.3843 

- .6589 

.63 

3.16 

.7927 

-1.0676 

2.18 

2.58 

.4150 

- .6899 

.96 


SAMPLE 

VARIANCE 

dB 


.4529 0.57 
.5554 1.01 
.6157 1.15 


.3908 0.84 

-1.469 11.8 


. *»c32 18.1 





















Experimental Values from Unflanged Annulus with Hub-Tip Ratio of 0.566 




Figure 144. Comparison of Experimental and Theoretical Values of Radiation 
Impedance for a Flanged Annulus 






The differences between measured and calculated values of reactance are 
much greater and in fact do not follow the calculated values. However, the 
trends predicted appear to be much the same and the, major difference appears 
to be the result of a frequency shift in the peak values of reactance. There 
is enough experimental data to indicate that the shape of the predicted reac- 
tance behavior with frequency is similar to the measured shape. The frequency 
shift in the peak is believed to be due to the difference between the flanged 
theory and the unflanged experiment. 

The measurements shown in the previous two figures were repeated at 50 
and 150 fps. Figures 145(a) to 145(f) show the comparison of the data and 
theory. ( In examining the measured data, it is first noted that the measured 
resistance values for the (1,0) mode with flow are negative. This Is to be 
expected for small values of kR as discussed by Mechel et al (ref. ^4), The 
effect of flow is to increase the reflected pressure when the incident wave is 
in the direction of flow. The reflected energy is not greater than the 
incident energy (i.e., there is no sound source at the duct termination). 
However, the flow changes the phase relationship between sound pressure and 
acoustic velocity, thus to maintain the reflected energy, the pressure ampli- 
tude increases with the phase change. Another way of describing it would be 
to say that the reflected incident pressure wave is "convect ively amplified," 
while the incident wave is "convect i vel y attenuated*" Therefore, the convec- 
tive effects on pressure result In a negative radiation resistance. At higher 
frequencies, the effect of flow seems to result in increased peaks and dips In 
radiation resistance curve. The radiation reactance data are not accurate 
enough to attempt a qualitative assessment of the effects of flow. 


VIII. RADIATED FIELD 


The level and spatial distribution of the radiated field resulting from 
an acoustic field in a duct has been studied relatively little, especially In 
respect to the effects of Jet velocity. There is no theoretical study known 
at present which allows for calculation of the radiated field from a duct mode 
or a complex pressure distribution as it radiates through a jet flow, 

A need exists for experimental data with which to compare against theoret- 
ical predictions, as the prediction techniques emerge In the next few years. 
Also, comparison with the radiation theory predictions without Jet flow Is an 
equally important requirement. Thus, this section presents a series of measur- 
ed radiated field patterns correlated with the in-duct conditions. In all 
cases presented, an attempt has been made to define the condition at the end 
of the duct. In some cases, where modes are excited In relative Isolation, an 
angular SPL plot at one radial location at the duct termination provides 
adequate definition of the in-duct condition. In other cases, where multiple 
modes were excited, a complete tabulat ion of the individual mode amplitudes 
and phase, as well as the angular position of the mode in the duct (In the case 
of higher order stationary modes), relative to the plane of the radiated field 
measurement is given. Also, in this case, the reflection coefficient and the 
radiation impedance is tabulated for future reference. 
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The primary variables are mode order, velocity and wave number. The most 
complete series of tests were run with multiple mode excitation. These were 
restricted to the zero order radial mode, but extend through the (4,0) angular 
mode. Since the modes are primarily stationary angular modes , the level in 
the radiated field is directly related to the radiated field microphone boom 
plane position, relative to the angular position of the stationary mode in the 
duct. The radiated field solution has an e im< J> characteristic shape, where <{> 
is the circumferential angle around the duct axis. 

It is recognized that the data to follow is not complete in scope, but is 
hoped that enough exists to act as a guide in formulating new theoretical con- 
cepts. Certainly, enough data will be presented to serve as a check for 
initial analytical efforts in this field of study. 


VI I 1.1. Description of Experiment 


The noise radiated from the duct termination was measured in an anechoic 
room as shown in figures 1, 3, and 6. The boom on which the microphone was 
mounted traversed a 124 inch radius from the center of the duct termination. 
The configuration is shown in figure 146. The boom was driven by a variable 
speed electric motor and angular position was monitored by a linear potentiom- 
eter. The microphone probe was the same configuration as shown in figure 64, 
in order to minimize the effects of the jet flow velocity on the acoustic 
measurement . 

The duct reference angle for the stationary angular modes is shown in 
figure 147. In the tabular data which is included with each figure, the 
position of the first node in a clockwise direction, from 0=0, is listed. 

The jet flow was a rather odd flow, since it had a hole in the middle in 
the mixing region. This is a result of only having flow in the annular 
section. The mean flow characteristics of the jet were not measured. The 
flow velocity given is that at the exit plane of the annular duct. 

Measurements below 500 Hertz are questionable, since the polyurethane 
foam pads on the walls were effectively anechoic only above this lower fre- 
quency limit. 


VI I 1.2. Experimental Results 


VI I 1.2.1. Effect of mode number and flow velocity at a fixed frequency 

Figures 148-151 are a series of radiated field plots for angular 
stationary modes in relative isolation. The components of other modes in each 
case is given in part (a) of the figure, which illustrates the SPL angular 
distribution in the duct. 
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In observing the effect of mode number on the radiated field, it is seen 
that the plane mode is the most effective radiator; but the main plane mode 
power is on the axis. As the order of the mode is increased, the SPL at 90° 
to the axis increases significantly. 


. P r < ma ry effect of jet flow is to refract the sound away from the jet 

axis. Thus, it would be expected that jet flow would have the greatest influ- 
ence on the plane mode since It is the only mode with a no-flow peak radiation 
angle on the jet axis. This is observed to be correct, as can be seen by care- 
ful examination of figures 148-151. A replot of the amplitude in these figures 
is condensed in figure 152(a)-(d) in the form of polar plots. This vividly 
illustrates the refractive effect of flow, particularly for the plane mode. 


v III* 2.2. Effect of Frequency and Multiple Modes 

Another series of tests was conducted in order to demonstrate the 
effect of both frequency and multiple modes on the radiated field pattern and 
efficiency. These data are shown in figures 153-157 from zero to 200 fps in 
increments of 50 fps. The results are similar to those of the previous sub- 
section. However, the presence of the multiple modes makes the pressure 
patterns non-symmet r i cal about the duct axis. 


The effect of flow is the same as in the previous sub-section at 1000 
Hertz however, in this series of tests, the radiated sound is of less inten- 
sity (the efficiency is lower). Thus, as the flow velocity increases, the 
turbulent pressure in the flow controls the microphone output when it resides 
in the flow at the low frequencies. It is obvious, however, that as frequency 
is increased, the refractive action due to the Jet flow becomes more pro- 
nounced. For the 1356 Hertz test in figure 157, refraction has reduced the 

level by some 14-15 dB at 15 off the iet avic wi*-h a o ao \ . 

• £ . r . Jet axls » with a 2 dB increase observed 

in the region of 50-60 off the axis. 


The general refractive effect of jet flow, as observed in analytical 
les on point source radiation, is to redirect the acoustic energy. This 
action generally observed in the experimental data. 


stud- 
i s the 


° f c ®' Jtion is in order concerning use of the experimental 
radiated field data. The pressure inside the flow region, which at 200 fps is 

' “ 'VT 66 ! a - OU ? the j6t axiS> is not due t0 the sound field radi- 
t ng from the duct. It is instead a pressure due to the turbulence in the 

f ow. Thus, data in this region should be ignored, when studying the refrac- 
tive effects due to flow velocity. 
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IX. CONCLUSIONS 


The primary conclusion reached from this study is that the theory des- 
cribing sound propagation in an acoustically lined annular duct sustaining a 
mean flow, as outlined In the first section of this report, has been verified 
for most effects of practical interest by the experiments reported herein. 

A listing of the major conclusions leading to this primary result as well 
as the subsidiary findings resulting from the experiments are listed below. 


1. A fundamental conclusion is that audio driver units can be effec- 
tively used in the simulation of the primary characteristics of 
fan-generated tone noise. If the internal and mutual impedances of 
these sources are at sufficiently high values, then interaction, 
effects between the sources are negligible. Subsidiary conclusions 
were that placement of the sources on the duct walls and design of 
the transducers to simulate point velocity sources were necessary 
for effective fan noise simulation. 

2. In a series of tests in hard wall annular ducts it was shown con- 
clusively that the fundamental results of the modal theory of sound 
propagation in that type of duct are valid. Mode cut-on frequencies, 
phase speed and radial and angular mode shapes were measured and the 
comparison with theoretical predictions was accurate. Effects of 
flow were also adequately predicted by the theory. 

3. Techniques were refined for measurement of wall impedance with the 
two microphone method and a modification of this method was devel- 
oped using one traversing microphone. These impedance measuring 
techniques were used to quantify the effect of flow on liner 
impedance. 

2*. The effect of flow on the impedance of locally reacting liners was 
to significantly increase the acoustic resistance of the liners at 
low frequencies. Flow effects on reactance were less significant 
but the trend was to change the reactance toward more positive 
values with increasing velocity. 

5. Another series of tests conducted in five different sets of acousti- 
cally lined annular ducts verified many analytical predictions 
concerning sound propagation in lined annular ducts, and at the same 
time revealed some of the problems to be expected in liner design. 
Mode shapes, attenuation and phase speed generally compared qui te 
well between measurement and prediction. However, it was found that 
some of the liners have drastic variations in acoustic properties 
across the liner surface. Variations of ±50% in acoustic resistance, 
were detected in one of the foam liners and in one type of the 
locally reacting liners (porous metal over a honeycomb core). Flow 
predictions were qualitatively verified. 
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6. A method for measuring the radiation impedance of higher order modes 
with flow was developed. This method produced new results and in 
the case of zero flow, verified the theoretical predictions. It was 
shown that, with flow, the radiation resistance could be negative, 
but this result did not imply an energy source at the termination. 
Instead, It resulted from phase changes between sound pressure and 
acoustic velocity caused by the mean flow through the duct. 

7. Radiated field patterns were presented for several modes and combi- 
nations of modes and the effect of Jet flow on refraction of the 
radiated sound field was demonstrated. 
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